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ABSTRACT 


The  dynamics  of  the  reactions  of  N(2d,  2p)  vith  O2  have  been  investigated 
by  observing  the  initial  vibrational  state  population  distributions  of  NOCX^II) 
formed  in  the  interaction  of  oxygen  and  discharge-excited  nitrogen  in  a 
cryogenically  pumped,  low  pressure  reaction  volume.  Infrared  chemilumines¬ 
cence  from  vibrationally  and  rotationally  excited  N0(v=l-14)  was  observed  in 
the  5  to  7  ym  spectral  region,  and  spectral  data  were  analyzed  to  obtain 
NO(v,J)  number  densities.  The  results  show  multimodal  vibrational  and  rota¬ 
tional  distributions,  indicative  of  several  concurrent  processes.  In 
particular,  extensive  rotational  excitation  is  evident  in  several  vibrational 
levels  of  NO,  resulting  in  sharply  peaked  R-branch  band  heads  characteristic 
of  J>80.  Surprisal  analysis  was  applied  to  the  data  to  infer  the  contri¬ 
butions  of  the  various  reaction  pathways.  This  analysis  indicates  that  all 
the  possible  reaction  pathways  for  producing  NOCX^n),  some  of  which  lead  to 
formation  of  metastable  atomic  oxygen,  occur  with  significant  probability. 

The  results  are  discussed  in  terms  of  their  implications  for  the  detailed 
collisional  dynamics  of  the  reactions,  and  their  relevance  to  the  study  of 
upper  atmospheric  auroral  chemistry. 

INTRODUCTION 

The  importance  of  the  reactions 


N(2d,  2p)  +  O2  -»  NO  +  0 


(1) 


as  the  principal  source  of  NO  in  the  Earth's  mesosphere  and  thermosphere  was 
suggested  by  Frederick  Kaufman  and  others  in  1967.^  Following  the  direct  rate 
coefficient  measurements  and  product  yield  determinations  of  Lin  and  Kaufman,^ 
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first  reported  in  preliminary  form  in  1970,^  atmospheric  modeling  studies'^ 
showed  that  the  NC^D)  reaction  accounted  for  the  large  upper  atmospheric 
production  of  NO  observed  in  rocket  experiments.^  A  number  of  subsequent 
laboratory  investigations^"^  have  provided  refinements  of  the  original^*^ 
kinetic  measurements,  and  atmospheric  modeling  studies^*^* have  continued  to 
examine  the  roles  of  these  reactions  in  odd  nitrogen  chemistry  of  the 
quiescent  and  auroral  upper  atmosphere. 

From  both  an  aeronomic  and  reaction  dynamic  point  of  view,  the  potential 
of  the  above  class  of  reactions  to  form  excited  state  products,  such  as  N0(v), 
0(^D,  ^S),  and  N0(a^n,  v),  is  of  great  interest.  In  a  detailed  experimental 
study,  Kennealy  et  al.^2  showed  that  Reaction  (1)  is  an  efficient  source  of 
highly  vibrationally  excited  N0(v)  (up  to  12  quanta  with  a  strongly  non- 
statistical  vibrational  state  distribution),  with  implications  that  0(^D) 
might  be  formed  as  well.  Rocketborne  infrared  spectra  of  a  strong  aurora 
showed  extensive  auroral  excitation  of  N0(v);^3  tliese  results  were  success¬ 
fully  modeled-^^^  in  terms  of  direct  formation  by  the  reaction  of  aurorally 
produced  NC^D)  with  O2  as  studied  in  the  laboratory. 

The  possible  reaction  pathways  for  Reaction  (1)  are  shown  in  Table  1. 

Also  shown  are  the  reaction  exothermicities,  the  maximum  energetically 
accessible  (rotationless)  vibrational  levels  in  the  product  NO  molecules,  and 
the  Cg  symmetry  species  for  adiabatic  conversion  of  reactants  to  products. 

The  symmetry  species  were  determined^*^^  using  the  adiabatic  correlation  rules 
for  atom-diatom  reactions  having  non-linear  intermediate  complexes,  as 
described  by  Shuler. Note  that,  despite  the  correlation  of  both  reactants 
and  products  to  doublet  and  quartet  A'  and  A"  surfaces,  the  reaction  of  N(2p) 
to  form  ground  state  atomic  oxygen  (.Reaction  3a)  has  no  adiabatic  route;  this 
results  from  application  of  the  so-called  "non-crossing  rule"  derived  by 
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As  we  will  discuss 


analogy  to  avoided  crossings  in  atom-atom  collisions. 
later,  this  is  not  sufficient  reason  to  rule  out  the  occurrence  of  this 
reaction  pathway.  Clearly,  if  significant  reaction  exothermicity  goes  into 
product  vibrational  excitation,  it  should  be  possible  to  identify  the  relative 
contributions  of  the  various  NOCX^II)  product  channels  by  examining  the 
infrared  fluorescence  "signatures"  in  the  fundamental  vibration-rotation 
transitions  in  the  5  to  8  um  spectral  region. 

The  objective  of  the  present  study  is  to  examine  the  detailed  vibrational 
and  rotational  energy  partitioning  in  N0(x2n)  formed  from  Reactions  (2) 
and  (3).  We  have  used  the  technique  of  nascent-product  infrared  chemilumines¬ 
cence  spectroscopy  under  near-single  collision  conditions  in  the  COCHISE  (COld 
CHemiluminescence  Infrared  Stimulation  Experiment)  cryogenic  reactor/ 
spectrometer  facility  at  the  Air  Force  Geophysics  Laboratory.  Previous 
investigation^^  of  these  reactions  on  the  COCHISE  apparatus  showed  extensive 
NO  vibrational  excitation,  principally  from  N(2d),  and  inferred  the  possi¬ 
bility  of  0(^D)  formation  via  Reaction  (2b)  from  surprisal  analysis  of  the  NO 
vibrational  distributions.  In  addition,  evidence  of  significant  rotational 
excitation  was  observed  at  the  lowest  pressures  employed.  However,  the 
results  were  affected  by  the  presence  of  unidentifiable  underlying  spectral 
features  at  the  short  and  long  wavelength  ends  of  the  spectra,  and  the  con¬ 
tributions  due  to  N(2p)  could  not  be  unambiguously  assessed.  In  the  present 
work,  we  have  taken  advantage  of  increased  sensitivity  and  dynamic  range  to 
remove  contributions  from  background  radiation  and  to  extend  the  results  to 
higher  vibrational  levels.  We  have  identified  previously  unknown  spectral 
features  as  R-branch  band  heads  resulting  from  extensive  rotational  excitation 
(J-100)  in  each  of  at  least  eight  vibrational  levels,  which  »e  attribute  to 
formation  by  N(2p).  We  have  inferred  detailed  product  channels  via  surprisal 
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analysis  of  the  NO  internal  energy  distributions;  the  results  Indicate  that 
the  N(2d)  and  N(2p)  reactions  with  O2  proceed  by  very  different  molecular 
dynamic  mechanisms. 

EXPERIMENTAL  MEASUREMENTS 

The  design  and  operation  of  the  COCHISE  facility  have  been  described  in 
detail  elsewhere. In  brief,  the  entire  radiative  environment  of  the 
experiment  is  maintained  at  a  base  temperature  of  app’"5ximately  20  K,  which 
effectively  eliminates  background  radiation  within  the  2  to  20  um  operating 
range  of  the  apparatus.  The  detection  system  consists  of  a  scanning  grating 
monochromator  and  a  liquid-helium-cooled  Si:As  detector.  Reagent  gases  are 
introduced  into  the  reaction  cell  through  temperature-controlled  feed  lines  as 
shown  in  Figure  1.  N(2d,  ^p)  are  produced  (together  with  N(^S)  and  excited 

metastable  states  ot  N2)  by  passing  a  flowing  N2/Ar  mixture  at  approximately 
1  torr  through  four  parallel  microwave  discharges  (2450  MHz,  50  W)  prior  to 
expanding  in'‘o  a  low-pressure  (approximately  3  mtorr),  cryogenically  pumped 
interaction  volume.  A  counterflow  gas  (O2  or  Ar)  enters  the  volume  from  the 
opposite  side  to  combine  with  the  discharge-excited  gas  in  a  stagnation  region 
along  the  axis  of  the  cell,  which  coincides  with  the  axis  of  the  cylindrical 
field  of  view  of  the  detection  system.  In  the  case  of  an  unreactive  counter¬ 
flow,  this  interaction  provides  an  arrested  relaxation  condition  for  studies 
of  infrared  fluorescence  arising  from  discharge  excitation  processes. 

However,  in  the  present  experiments,  an  O2  counterflow  reacts  promptly  in  the 
field  of  view  with  the  effluents  of  the  N2/Ar  discharges  to  give  readily 
identifiable  chemiluminescence  from  vibrationally  excited  NO  under  conditions 
where  diffusive  loss  dominates  over  radiative  or  collisional  cascade. 
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For  most  experiments,  reagent  gases  were  introduced  at  controlled 
temperatures  of  80  to  90  K;  no  temperature  dependence  in  the  fluorescence 
intensity  or  general  spectral  distribution  was  observed  at  temperatures  up  to 
120  K.  In  each  discharge  tube,  the  Ar  flow  rate  was  held  at  640  ymole/s,  and 
the  N2  flow  rate  was  systematically  varied  from  5  to  86  ymole/s  (0.8  to 
12  percent  N2).  Mass-balanced  counterflows  of  either  O2  or  Ar  were  used  to 
produce  chemiluminescence  or  determine  background  fluorescence  levels,  respec¬ 
tively.  Under  these  conditions,  the  gas  residence  time  in  the  active  dis¬ 
charge  is  approximately  2  to  5  ms,  the  transit  time  from  the  discharge  exit  to 
the  field  of  view  is  -  0.5  ms,  and  the  residence  time  in  the  field  of  view  is 
~  0.3  ms . 

Infrared  fluorescence  spectra  were  obtained  in  the  4  to  8  ym  region  with 
spectral  resolutions  of  0.007  to  0.040  ym  (2.8  to  16  cm“^  at  5  ym) .  Some  50 
spectral  scans  were  required  to  survey  the  parameter  space  of  N2  mole  frac¬ 
tion,  O2  mole  fraction,  temperature,  discharge  power,  and  spectral  resolution 
(including  foreground  and  background  scans).  The  uncertainty  in  the  observed 
wavelengths  (due  to  a  slight  periodic  fluctuation  in  the  cryogenic  mono¬ 
chromator  scan  rate)  is  ±  0.003  ym;  absolute  wavelengths  were  calibrated  using 
the  well-resolved  band  center  of  the  C0(l-*0)  vibrational  transition  excited  by 
energy  transfer  from  discharge-excited  nitrogen. The  spectral  data  were 
acquired  via  phase  sensitive  detection,  using  a  computer-interfaced  lock-in 
amplifier  and  chopping  the  discharges  with  a  23  Hz  square  wave.  The  observed 
intensities  were  corrected  to  radiometric  units  using  blackbody  calibrations 
of  the  absolute  spectral  responsivity  of  the  detection  system. 
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RESULTS  AND  ANALYSIS 


Spectral  Data 

An  example  fluorescence  spectrum  observed  in  the  interaction  of  dis¬ 
charged  N2/Ar  with  O2  is  shown  in  Figure  2  for  the  highest  usable  spectral 
resolution.  Spectral  scans  in  which  the  O2  in  the  counterflow  was  replaced  by 
Ar  revealed  the  presence  of  weak  underlying  radiation  from  various  sources 
throughout  the  bandpass  of  interest  (4  to  8  ym).  These  consist  of:  1)  Ar(I) 
Rydberg  lines  between  4  and  5  ym;  2)  N2(W^Ay-  B^ng)  bands  near  6.5  and  7.5  ym 
(the  (1,0)  and  (2,1)  bands,  respectively);  and  3)  broad  but  structured 
emission  extending  from  5.2  ym  to  approximately  7.5  ym.  We  have  previously 
investigated  the  first  two  emission  systems  in  some  detail:  the  Ar(I) 
radiation^l  is  scattered  light  from  the  discharges,  and  the  N2(W-B)  emission^^ 
arises  from  discharge-excited  species  which  survive  long  enough  to  radiate  in 
the  field  of  view.  The  third  background  contribution  occurs  with  signal-to- 
noise  ratios  which  are  too  low  to  permit  conclusive  identification;  however  it 
appears  that  this  feature  arises  from  N0(v)  formed  in  the  discharges  due  to 
ppm-level  oxygen  impurity  in  the  discharged  gas.  The  fluorescence  spectral 
distribution  is  "collapsed"  toward  shorter  wavelengths  (lower  vibrational 
levels),  indicative  of  substantial  collisional  deactivation  as  expected  for 
the  -  1  torr  discharges.  The  maximum  intensity  of  this  feature  is  about  one 
order  of  magnitude  less  than  that  of  the  chemiluminescence  signal,  but  the 
contribution  to  the  observed  spectrum  from  the  background  emission  is 
significant  for  wavelengths  below  5.4  ym  and  above  6.4  ym  (v<2  and  v>12, 
respectively).  Accordingly,  the  chemiluminescence  spectra  were  corrected  for 
background  contributions  by  subtracting  companion  background  spectra  obtained 
immediately  before  or  after  each  scan,  with  all  conditions  identical  except 
for  the  substitution  of  an  Ar  counterflow  instead  of  O2.  An  unfortunate 
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result  of  this  subtraction  procedure  was  to  limit  the  dynamic  range  at  long 
wavelengths  and  hence  the  maximum  observable  vibrational  level.  The  bulk  of 
the  remaining  chemiluminescence,  between  5.2  and  6.5  urn,  can  be  readily  iden¬ 
tified  through  its  vibrational  and  rotational  spacings  as  N0(x2n,  v'=l-14). 

The  maximum  spectral  intensity  occurs  for  v'=5-7,  and  the  rotational  distri¬ 
butions  in  this  portion  of  the  spectrum  signify  near-thermal  equilibrium  with 
the  gas  temperature.  Also  present  in  the  background-corrected  spectra  is  a 
series  of  eight  sharp,  red-degraded  features  between  4.95  and  5.55  ym  (cf. 
Figure  2).  These  features,  which  have  a  constant  spacing  of  31±2  cm“^,  are 
reproducible  in  wavelength  and  relative  intensity  for  all  NO  chemiluminescence 
experiments.  The  features  between  4.9  and  5.2  ym  were  noted  in  earlier 
Investigations; however,  poorer  signal-to-noise  ratios  precluded  systematic 
examination,  and  the  usable  spectral  resolution  in  those  experiments  was  too 
low  to  feveal  the  presence  of  the  additional  features  extending  into  the  main 
NO  band.  It  is  clear  that  the  presence  of  these  features  could  significantly 
affect  the  vibrational  analysis  of  the  NO  spectrum.  Furthermore,  the  three 
features  between  4.95  and  5.15  ym  are  coincident  with  similar  spectral 
features  recently  observed  by  a  rocket-borne  Interferometer  in  an  overhead 
aurora, 22  suggesting  that  this  emission  represents  a  previously  unknown 
aeronomic  phenomenon. 

The  observed  spectral  features  do  not  correspond  to  any  of  the  possible 
transitions  that  might  be  excited  in  the  experiments,  e.g.,  vibrationally 
excited  NO"*"  or  N2O,  or  electronic  transitions  of  NO,  N2,  and  O2.  However,  it 
is  straightforward  to  show  that  the  observed  band  positions  correspond  well  to 
those  predicted  for  R-branch  band  heads  of  the  vibrational  transitions  of 
N0(x2n)  itself.  The  presence  of  such  band  heads  in  a  vibrational  transition 
requires  an  extraordinary  degree  of  rotational  excitation.  For  the  (v',v")  = 
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(1,0)  transition,  for  example,  the  vertex  of  the  Fortrat  parabola  occurs  for 
J'-80,  a  rotational  energy  in  excess  of  1.3  eV.  Thus  it  appears  that  che 
vibrational  transitions  for  at  least  v'=l-8  have  strongly  bimodal  rotational 
distributions,  with  Boltzmann  "temperatures"  ranging  from  -  100  to  ~  10000  K. 
The  rotational  states  contributing  to  the  R-branch  band  heads  (J'~  60-100) 
should  also  exhibit  P-branch  structure  extending  to  8  or  9  pm;  unfortunctely , 
due  to  the  low  signal-to-noise  ratios  at  the  longer  wavelengths,  we  cannot 
examine  the  P-branch  distributions  to  determine  the  detailed  rotational 
population  distributions  at  large  J'.  Furthermore,  due  to  extensive  overlap 
by  the  relatively  intense,  rotationally  ' thermal"  component  beyond  5.6  pm,  we 
cannot  determine  whether  there  are  weak  rotationally  excited  components  for 
vibrational  levels  v'>8.  By  comparison,  in  experiments  where  CO/Ar  counter¬ 
flows  are  used,  we  observe  similar  band  head  formation  in  several  vibrational 
levels  of  the  C0(Av=l)  system  from  rotational/vibrational  excitation  by  energy 
transfer  from  active  nitrogen; ^0  these  measurements  have  sufficient  intensity 
to  permit  resolution  of  the  P-branch  structure  at  large  J' ,  and  show  that  a 
two-temperature  bimodal  description  provides  a  reasonable  approximation  to  the 
rotational  state  populations.  In  order  to  examine  trends  in  the  relative 
intensities  of  the  "hot"  and  "cold"  rotational  components  of  NO,  as  well  as  in 
the  vibrational  distributions,  several  experimental  parameters  were  inde¬ 
pendently  varied.  Spectra  were  obtained  over  the  ranges:  discharge  power, 

10  to  50  W;  gas  temperature,  80  to  120  K;  counterflow  O2  mole  fraction, 

0.0  to  1.0  (Ar  diluent);  and  discharge  N2  mole  fraction,  0.008  to  0.12  (Ar 
diluent).  Increasing  the  discharge  power  gave  increased  signal  intensity  with 
no  discernible  effect  on  the  spectral  distributions.  Increasing  the  gas 
temperature  did  not  affect  the  observed  spectral  intensities  or  distributions 
except  for  slight  increases  in  the  "thermal"  rotational  temperatures. 
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Reduction  of  the  counterflow  O2  mole  fraction  gave  a  corresponding  decrease  in 
spectral  intensity,  with  no  significant  effect  on  the  spectral  distributions. 
Reduction  in  the  discharge  N2  mole  fraction  by  a  factor  of  15  gave  a  factor  of 
~  2  decrease  in  overall  signal  intensity,  with  ~  30  percent  apparent  reduction 
in  the  relative  high-J'  contributions.  This  latter  effect  is  significant  for 
the  present  observations,  and  will  be  discussed  further  in  terms  of  the 
derived  vibrational  state  populations  below. 

Spectroscopic  Analysis 

The  spectral  data  were  analyzed  to  determine  detailed  vibrational  state 
populations  using  a  linear  least  squares  spectral  fitting  method  which  we  have 
employed  for  a  number  of  vibrational  excitation  studies. ^2, 13, 16-18  The 
rethod  consists  of  computing  instrument-convolved  basis  functions  for  each 
Av=l  vibration-rotation  transition  using  an  accurate  spectroscopic  formalism, 
and  fitting  those  functions  to  an  observed  spectrum  to  determine  the  indi¬ 
vidual  integrated  intensity  of  each  vibrational  transition.  These  values  can 
then  be  transformed  into  upper  vibrational  state  number  densities  through 
application  of  the  appropriate  band  transition  probabilities. 

The  vibration-rotation  energies  and  transition  frequencies  of  NOCX^Il) 
were  computed  explicitly, including  the  effects  of  spin-orbit  coupling, 
vibrational  anharmonicity,  vibration-rotation  interaction,  and  centrifugal 
distortion;  effects  of  electronic  rotational  interaction  (X-type  doubling)  are 
negligible  for  these  experiments  (<  0.1  cm~^).  The  spectroscopic  constants 
used  for  the  computations  are  shown  in  Table  2.  For  the  most  part,  these  are 
the  same  values  given  by  Goldman  and  Schmidt^^  and  used  by  Rawlins  et  al.^3 
fit  high  resolution  auroral  spectra.  The  single  exception  is  the  value  of  3g, 
the  correction  term  to  the  centrifugal  constant  Dy,  which  we  adjusted  in  order 
to  better  reproduce  the  positions  of  the  R-branch  band  heads.  The  value  of  8g 
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given  in  Table  2  is  about  twice  that  used  by  Rawlins  et  al,l^  who  reversed  an 
apparent  error  of  sign  in  the  value  given  by  Goldman  and  Schmidt. 

Comparisoi  of  the  values  of  Goldman  and  Schmidt  to  subsequent  work  by  Amiot 
and  Verges^^  indicates  a  2  percent  uncertainty  in  Dy  at  room  temperature, 
which  encompasses  the  change  we  require  for  Pg. 

For  the  limited  spectral  resolution  of  these  experiments,  a  complete 
description  of  the  line-by-line  rotational  intensities  is  not  feasible,  nor 
does  adequate  information  exist  for  extending  these  intensities  to  the  high 
rotational  levels  observed  here.  As  an  approximate  treatment,  we  have  used 
band-averaged  Einstein  coefficients  together  with  Honl-London  scaling  of  the 
relative  rotational  strengths.  The  relative  intensity  of  the  i*^^  optically 
thin  rotational  line  emanating  from  the  state  (v',J',S')  is  then  given  by 


I,  =  he 


"i 


S.(P,Q,R) 


A  ,  „  N  , 

y  /  Y**  y ' 


exp(-E(J' ,2' )/kT) 


where  vj  and  vg  are  the  transition  frequencies  of  the  line  and  the  band 
center,  respectively,  Si(P,Q,R)  are  the  appropriate  Honl-London  factors  as 
given  by  Herzberg,23  Ay/yii  is  the  thermally  averaged  Einstein  coefficient  for 
the  band,  Nyf  is  the  population  of  the  emitting  vibrational  level,  and  Q(v' ) 
is  the  spin-orbit/rotational  partition  function  for  the  emitting  vibrational 
level.  Note  that  here  we  have  assumed  that  the  ^^1/2  and  ^113/2  spin-orbit 
states  are  in  thermal  equilibrium,  which  in  fact  seems  to  be  the  case,  at 
least  for  the  rotationally  thermal  component  of  the  spectrum  where  spin  non¬ 
equilibrium  would  be  the  most  apparent. 

The  individual  line-by-line  vibrational  basis  sets  are  convolved  with  the 
instrumental  scan  function  (in  this  case,  an  isosceles  triangle  with  full 
width  at  half  maximum  as  the  spectral  resolution)  and  are  fit  to  each 
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experimental  spectrum  by  linear  least  squares  to  determine  the  products 
Nv'Av»vti.  The  individual  upper  state  number  densities  are  then  determined  by 
application  of  band-averaged  Einstein  coefficients  for  each  transition  (see 
below) . 

In  the  fitting  procedure,  the  bimodal  rotational  distributions  in  each  of 
the  lowest  eight  vibrational  levels  were  approximated  by  two  separate  vibra¬ 
tional  populations  with  Boltzmann  rotational  distributions  corresponding  to 
thermal  (-  100  K)  and  "hot"  (~  10000  K)  temperatures.  These  temperatures  were 
independently  varied  to  optimize  the  fits  to  the  appropriate  rotational  con¬ 
tours,  and  were  found  to  be  typically  120+10  and  10000+4000  K.  The  quoted 
uncertainties  are  based  on  the  variations  required  to  produce  recognizable 
deviations  from  the  observed  bandshapes  near  5.6  to  6.0  um  and  4.9  to  5.2  ym 
for  the  "cold"  and  "hot"  components,  respectively.  Within  these  ranges  of 
temperature  parameters,  the  vibrational  distributions  determined  from  the 
spectral  fits  were  not  significantly  affected.  We  attempted  to  determine 
whether  the  rotational  temperature  for  the  "hot"  component  was  dependent  on 
vibrational  level.  While  there  was  some  indication  that  the  R-branch  band 
head  structures  required  somewhat  lower  rotational  temperatures  for  the  higher 
vibrational  levels,  the  limitations  of  signal-to-noise  ratio,  spectral  reso¬ 
lution,  and  band  overlap  precluded  a  quantitative  determination  of  this  trend. 

The  band-averaged  Einstein  coefficients  for  each  vibrational  transition 
were  chosen  by  logic  analogous  to  that  recommended  by  Rothman  et  al.26  The 
value  for  the  (v',v")  =  (1,0)  band  was  taken  from  the  absorptivity  measurement 
of  Holland  et  al.,^^  and  the  scaling  to  higher  v'  followed  relative  values 
predicted  by  Billingsley^®  from  ab  initio  calculations  of  the  dipole  moment 
function.  This  relative  scaling  is  supported  by  experimental  measurements  of 
the  (Av=2)/(Av=l)  branching  ratios  for  v'<7,29  and  the  resulting  values  of 
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ky/yn  appear  to  be  reasonably  accurate  (approximately  +  30  percent  or  better). 
Recent  COCHISE  measurements  of  ( Av=2)/(Av=l)  ratios  extending  to  higher  v' 
indicate  more  substantial  uncertainty  in  the  theoretical  dipole  moment  func¬ 
tion  for  but  the  relative  values  of  Ay/yu  for  the  fundamental  band  are 

probably  not  severely  affected  by  this  uncertainty.  The  values  of  Ay/yti, 
determined  by  the  above  procedure  and  used  in  the  analysis  of  all  data  dis¬ 
cussed  in  this  paper,  are  listed  in  Table  3.  It  should  be  noted  that, 
although  the  band-averaged  Einstein  coefficients  are  nearly  independent  of 
rotational  temperature  up  to  at  least  600  K,28  there  may  be  effects  at  high  J' 
due  to  rotation-vibration  coupling  and  uncoupling  of  spin  and  orbital  angular 
momentum  (transition  to  Hund's  case  (b)),  as  noted  by  Billingsley^®  for  J'~30. 
To  test  for  this  possibility,  we  performed  more  detailed  computations 
accounting  for  spin  uncoupling  which  indicate  that  the  Case  (a)  Honl-London 
factors  provide  a  good  approximation  up  to  J-100.  More  significantly, 
however,  the  high-J'  states  forming  the  band  heads  are  sufficiently  high  in 
the  NO  manifold  to  probe  the  range  of  internuclear  distance  where  the  calcu¬ 
lations  of  Billingsley  indicate  the  onset  of  downward  curvature  in  the  dipole 
moment  function,  and  where  our  recent  (Av=2)/(av=l)  data®®  indicate 
uncertainty  in  this  calculated  behavior.  Thus  the  computed  relative 
rotational  intensity  scaling  may  not  be  applicable  to  the  R-branch  band  heads; 
the  rotationally  ”hot"  vibrational  populations  determined  in  this  way  are  then 
somewhat  uncertain.  A  more  detailed  determination  of  the  absolute  values  of 
Ayfyit(J'),  as  well  as  experimental  determination  of  the  transition  moment  at 
large  internuclear  separation,  are  currently  in  progress. 

Vibrational  Population  Distributions 

A  typical  comparison  of  observed  and  computed  least-squares  fit  spectra 
is  illustrated  in  Figure  3,  together  with  the  contributions  of  the 
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rotationally  "hot"  and  "cold"  components.  The  values  of  NyfAy/yi*  determined 
from  the  fit  are  plotted  in  Figure  4.  The  error  bars  indicate  the  statistical 
standard  deviations  (±  la)  in  the  solutions.  Application  of  the  Einstein 
coefficients  of  Table  3  gives  the  upper  state  vibrational  number  densities 
with  the  same  relative  standard  deviations.  Uncertainties  in  the  Einstein 
coefficients,  which  cannot  be  evaluated  statistically  and  are  the  same  for  all 
the  data  analyzed,  are  not  shown  in  our  plots  of  the  vibrational  populations. 

In  general,  relative  vibrational  population  distributions  obtained  by 
omitting  the  background  correction  procedure  and  the  "hot"  rotational  con¬ 
tributions  gave  excellent  agreement  with  values  determined  earlier  in  COCHISE 
for  v' =1-1212  vith  values  inferred  for  Reaction  (2),  v'=l-7  from  electron 
beam  excitation  of  air. 31  Application  of  background  correction  reduces  the 
apparent  relative  population  in  v'=l  by  about  10  percent  and  that  for  v'=12  by 
more  than  a  factor  of  two  from  the  previous  determination;  in  addition  the 
present  experiments  provide  information  for  levels  up  to  v'=14.  Inclusion  of 
the  rotationally  "hot"  states  gives  a  further  reduction  of  the  "cold"  v'=l 
population  by  a  factor  of  ~  2,  with  only  minor  effects  on  the  higher  levels. 

As  described  above,  the  relative  distributions  were  not  affected  significantly 
by  variations  in  experimental  parameters.  Accordingly,  the  solutions  for  a 
series  of  five  spectra,  obtained  sequentially  with  spectral  resolution  of 
0.013  ym  and  carefully  optimized  signal-to-noise  ratio,  were  averaged  together 
to  give  a  set  of  absolute  vibrational  state  populations  corresponding  to  the 
condition  T  =  80  K,  discharge  power  =  50  W,  and  discharge  N2  mole  fraction 
=  0.12,  to  be  used  for  subsequent  kinetic  analysis.  This  composite  averaged 
distribution  is  shown  in  Figure  5,  and  is  well  representative  of  relative 
distributions  determined  for  all  experimental  conditions.  The  remainder  of 
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this  paper  deals  with  kinetic  interpretations  of  the  results  in  Figure  5  to 
determine  product  branching  fractions  for  the  reactions  in  Table  1. 


KINETIC  INTERPRETATIONS 

Based  on  the  results  of  previous  studies, ^2, 31  it  is  reasonable  to  expect 
that  the  rotationally  "thermal"  vibrational  distribution  in  Figure  5  arises 
largely  from  the  reaction  of  N(2d)  with  02*  However,  the  rotationally  "hot" 
vibrational  distribution  may  arise  from  another  process  involving  reaction  or 
energy  transfer  of  other  excited  species  formed  in  the  discharge.  The 
principal  energetic^species  emanating  from  the  Ar/N2  microwave  discharges  for 
our  fast  flow  conditions  and  surviving  into  the  field  of  view  are  N(^S,  ^D, 

2p)  and  N2((X^Eg+,v) ,  A3Ey+,  w3A^,  a'^Ey",  The  flux  of  charged 

species  escaping  the  discharge  region  is  negligibly  small;  this  is  verified  by 
the  observation  of  no  excitation  in  CO  and  N2  counterflovs  upon  interaction 
with  the  effluents  of  lOOX  Ar  discharge.  Based  on  a  combination  of  discharge 
modeling  calculations  (using  methods^®  and  electron-impact  cross  sections^^ 
described  elsewhere),  direct  observations  of  discharge  effluents  in  a  flow 
reactor  at  room  temperature, and  observations  of  infrared-radiating  states 
of  N2  in  COCHISE, we  can  assemble  a  set  of  anticipated  discharge  effluent 
number  densities  for  the  present  experiments  at  a  discharge  N2  mole  fraction 
of  ~  0.1,  as  shown  in  Table  4.  The  molecular  metastables  arise  principally 
from  direct  excitation  of  N2  by  electrons  (characteristic  electron  energies 
are  6  to  8  eV  in  dilute  N2/Ar  mixtures)  and  metastable  argon,  and  the  atomic 
metastables  are  produced  primarily  by  electron-impact  excitation  of  N(^S) 
formed  from  N2  dissociation.  Clearly,  N(^S)  is  the  dominant  species,  but  is 
too  unreactive  with  O2  at  these  temperatures  to  produce  significant 
N0(v).12,34  Similarly,  we  do  not  expect  any  of  the  molecular  metastables  to 


640 


react  with  O2  to  form  NO(v),  as  a  four-center  collision  complex  would  be 
required.  An  interesting  possible  route  for  NO(v)  excitation  is  via  energy 
transfer  from  N2(A)  to  that  is  formed  by  N*  reactions  in  the  field  of 

view; 


N2(A3ry+)  +  N0(x2n,v)  -»  N2(Xlrg+)  +  N0(a2i:+.v) 

N0(a2i+,v)  ->N0(x2n,v')  +  hv 

The  NO(v')  distribution  that  is  formed  by  the  rapid  radiative  cascade  step  is 
governed  by  the  Franck-Condon  factors  for  the  transition  and  is  similar  to 
that  observed  for  the  rotationally  "hot"  component  in  Figure  5.  However,  due 
to  the  low  number  densities  of  both  N2(A)  and  NO  in  the  interaction  volume, 
the  NO(v')  number  densities  obtained  from  this  process  are  about  four  orders 
of  magnitude  smaller  than  those  observed  in  the  measurements.  Furthermore, 
although  rotational  excitation  of  N0(A)  has  been  observed  in  this  process, 
there  is  insufficient  exoergicity  to  account  for  the  much  greater  rotational 
excitation  observed  here.  By  similar  kinetic  arguments,  it  is  straightforward 
to  eliminate  other  electronic  energy  transfer  processes,  as  well  as  the 
possibility  of  rotational  excitation  by  energy  transfer  to  NO  from 
translationally  excited  atoms  emanating  from  the  discharge.  Vhile 
translational  excitation  of  atoms  is  known  to  occur  by  energy  transfer  from  Ar 
metastables, » 32  this  energy  is  moderated  within  microseconds  in  the  1  torr 
discharge  tubes,  and  the  number  density  of  "hot"  atoms  reaching  the  field  of 
view  will  be  miniscule.  Finally,  the  possibility  for  collisional  quenching 
from  a  high  vibrational  level  of  NO  into  a  high  rotational  level  of  a  low 
vibrational  level  (V,R  transfer),  analogous  to  processes  observed  for  HF,38 
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can  also  be  dismissed.  To  excite  the  high  J'  levels  we  observe  would  require 
quenching  of  vibrational  levels  5  to  12  or  higher.  Since  we  observe  very  low 
relative  populations  for  v>10  under  conditions  of  20  to  30  collisions  per 
observed  N0(v)  molecule,  this  would  require  essentially  gas-kinetic  quenching 
of  N0(v)  by  Ar.  This  is  highly  unlikely,  since  direct  measurements  of  the 
v-dependent  quenching  of  N0(v)  by  02^^  and  of  C0(v)  by  C02^^  indicate  quench¬ 
ing  rate  coefficients  (by  whatever  mechanism)  which  are  less  then  1  percent 
gas-kinetic  even  for  seven  or  more  vibrational  quanta  of  excitation.  Thus  any 
amount  of  rotational  excitation  produced  by  V,R  transfer  will  be  insignificant 
in  our  experiments.  Similarly,  the  observed  vibrational  distributions  are  not 
perturbed  by  collisional  deactivation  by  0,  O2,  N2,  N,  or  Ar,  as  demonstrated 
in  the  earlier  work  of  Kennealy  et  al.^^  We  conclude  that  the  observed  vibra¬ 
tional  distributions  represent  composites  of  the  initial  product  distributions 
from  the  reactions  of  N(2d,  ^P)  listed  in  Table  1.  Based  on  the  results  of 
Ref.  31,  the  rotationally  "cold"  distribution  appears  to  arise  largely  from 
the  N(2d)  reactions;  by  process  of  elimination,  the  rotationally  "hot" 
distribution  appears  to  result  from  N(2p)  reactions  with  02- 

Due  to  the  composite  nature  of  these  product  state  distributions,  it  is 
difficult  to  assign  the  individual  reaction  channels  unambiguously.  However, 
it  is  instructive  to  examine  the  observed  composite  distributions  with  the  aid 
of  surprisal  theory.  This  approach  postulates  the  existence  of  an  exponential 
gap  law  for  an  individual  metathetic  reaction;^^ 

P(v)  =  P0(v)  exp(-Xfv)/exp(Xo) 

vhere  fy  is  the  fraction  of  the  reaction  exoergicity  appearing  as  vibrational 
energy  in  the  product,  P(v)  is  the  observed  relative  vibrational  population. 
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and  P°(v)  is  the  statistical  "prior”  distribution  obtained  when  all  final 
translational,  rotational,  and  vibrational  states  are  equally  probable.  Thus 
the  ratios  of  the  observed  and  statistical  vibrational  populations  should  be 
exponential  in  the  vibrational  energy  Ey,  with  the  exponential  fall-off  con¬ 
stant  X  quantifying  the  departure  of  the  observed  distribution  from  a  com¬ 
pletely  statistical  product  distribution.  While  this  theory  does  not 
necessarily  hold  for  all  chemical  reactions,  it  does  provide  useful  guidelines 
for  the  interpretation  of  product  state  distributions  in  terms  of  the 
molecular  reaction  dynamics.  The  prior  distributions  were  determined  from  the 
usual  relationship  for  a  vibrating  rotator:^^ 


3/2 


v=0 


3/2 


Computed  prior  distributions  are  shown  in  Figure  6  for  the  reactions  in 
Table  1. 

Application  of  surprisal  analysis  to  the  vibrational  distributions  of 
Figure  5  is  by  no  means  straightforward  as  it  is  impossible  to  obtain  linear 
plots  of  ln(P(v))  versus  fy  for  any  single  assumed  elementary  reaction  from 
Table  1.  Thus  it  appears  that  each  distribution,  thermal  and  rotationally 
excited,  is  a  composite  of  distributions  for  more  than  one  reaction  channel. 
Previous  work,  employing  differential  quenching  of  N(^D)  and  N(2p)  by  N2  in 
electron-irradiated  N2/O2  mixtures,  provides  strong  evidence  that  the 
inverted,  rotationally  thermal  distribution  of  Figure  5a  is  due  to  the  NC^D) 
reaction  channels. We  therefore  believe  that  the  rotationally  excited 
distribution  of  Figure  5b  results  from  reactions  of  N(2p)  with  O2.  The 
rotational  excitation  required  for  band  head  formation,  J  ~  80-120,  is  about 
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2-3  eV.  This  plus  the  number  of  vibrational  levels  observed  (at  least  8) 
gives  a  reaction  exoergicity  of  4-5  eV  that  is  only  satisfied  by 
Reaction  (3a).  A  combination  of  Reactions  (3b)  and  (3c),  with  roughly  equal 
contributions  from  each  in  approximately  statistical  distributions,  can  also 
account  for  the  relatively  large  populations  observed  in  the  v=l,2  states. 
Furthermore,  a  statistical  combination  of  all  three  reaction  paths  (3a,  3b, 
and  3c)  in  equal  contributions  gives  an  excellent  fit  to  the  data,  as  shown  in 
Figure  7.  It  is  also  interesting  to  note  from  Figure  7  that  the  observed 
populations  of  v'=12-]4  from  the  rotationally  "thermal"  distribution  also  fall 
on  the  statistical  curve  for  Reaction  (3).  This  does  not  signify  a  discrepancy 
in  rotational  excitation  effects,  because  the  spectral  contributions  from 
these  states  are  so  weak  that  it  is  not  possible  to  identify  their  rotational 
distributions  with  any  accuracy,  and  we  cannot  determine  whether  their 
production  is  associated  with  the  rotationally  thermal  reaction  channel  (i.e.. 
N(2d))  or  with  the  rotationally  "hot"  reaction  (N(2p)).  As  will  be  shown 
below,  the  data  for  v=13-14  seem  less  consistent  with  the  trends  identified 
for  production  by  N(2d),  so  it  is  possible  that  the  observed  populations  for 
these  states  arise  primarily  from  Reaction  (3).  The  intermediate  states 
v=9-ll  should  also  be  produced  in  this  reaction,  but  at  levels  which  are 
obscured  by  production  of  NO(v)  from  Reaction  (2).  Similarly,  population  of 
levels  v>14  by  Reaction  (3a)  may  occur  but  cannot  be  observed  in  these 
experiments  because  of  signal-to-noise  limitations.  It  should  be  noted  that, 
while  the  contributions  of  Reactions  (3b,  c)  are  implied  by  the  fit  in 
Figure  7,  they  can  by  no  means  be  demonstrated  conclusively  from  these  data, 
as  we  do  not  have  any  positive  identification  of  the  energetic  0*  products  of 
these  reactions.  Nevertheless,  the  inclusion  of  these  reactions  is  consistent 
with  the  data  and  with  the  picture  of  a  completely  statistical  product 
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distribution  for  the  Nf^P)  reaction,  as  will  be  discussed  below. 

Alternatively,  the  entire  rotationally  hot  distribution  may  be  due  to  Reaction 
(3a)  alone,  with  non-linear  surprisal  behavior. 

We  examine  the  "thermal"  component  of  the  observed  NO(v)  distribution  by 
first  noting  that  correction  must  be  made  for  partially  thermalized  contribu¬ 
tions  from  the  reaction  branch  which  produces  rotationally  "hot"  N0(v).  These 
contributions  are  especially  evident  for  levels  v=l-3,  where  the  rotationally 
thermal  populations  closely  parallel  those  for  the  rotationally  excited 
products.  To  illustrate  this  effect,  we  show  in  Figure  8  the  surprisal  plots 
for  Reactions  (2a,  2b)  where  no  correction  is  made  for  thermalized  high-J 
contributions.  This  plot  shows  two  interesting  effects:  1)  the  populations 
for  v'=l,2  deviate  significantly  from  the  trend  for  v'=:3-7;  and  2)  Reaction 
(2a)  appears  to  make  no  contribution  for  v'<8,  as  evidenced  by  the  break  in 
slope  at  v'=7.  This  general  behavior  was  noted  in  the  surprisal  treatment  of 
similar  N0(v)  distributions  in  the  previous  work  of  Kennealy  et  al.^2 

Since  each  initial  NO(v,J)  molecule  is  subjected  to  20  to  30  collisions 
(mostly  with  Ar  and  O2)  during  its  residence  time  in  the  field  of  view,  it  is 
reasonable  to  expect  substantial  rotational  thermalization,  with  deactivation 
of  low  to  middle  rotational  levels  occurring  at  nearly  every  collision  and 
deactivation  of  the  high  rotational  levels  proceeding  more  slowly.  Our  data 
for  v'=l-3  indicate  that,  at  most,  half  of  the  rotationally  excited  N0(v,J) 
initially  formed  is  thermalized  within  30  collisions.  This  corresponds  to  an 
effective  rate  coefficient  on  the  order  of  ~  10“^^  cm^  s~^,  in  agreement  with 
a  ^.imilar  value  obtained  by  Taherian  et  al.^^  for  thermalization  of  highly 
rotationally  excited  N0(v,J)  formed  from  NO2  photofragmentation  in  He.  To 
correct  the  "cold'"  distributions  for  this  contribution  from  the  "hot"  reaction 
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channel,  we  assume  as  an  upper  bound  that  essentially  all  the  v'=l  population 
is  due  to  thermalized  NO(v,J),  scale  the  "hot"  populations  accordingly  to 
approximate  the  contributions  from  the  "hot"  channel,  and  subtract.  The 
maximum  correction  thus  obtained  is  large  for  low  v'  but  is  minor  near  the 
peak  of  the  distribution. 

The  remaining  population  distribution,  plotted  in  Figure  9,  shows  an 
inversion  peaking  at  v'=7.  Noting  that  v'=8  is  the  thermodynamic  limit  for 
Reaction  (2b),  we  show  in  Figure  10  the  combined  surprisal  plots  for 
Reaction  (2b)  (v'=2-7)  and  Reaction  (2a)  (v'=8-ll).  The  composite  population 
distribution  determined  from  the  surprisal  fits  of  Figure  10  is  given  by  the 
dashed  curve  in  Figure  9.  Excellent  linearity  in  the  surprisal  analysis  is 
now  evident,  as  is  the  implication  that  Reaction  (2b)  accounts  for  all  the 
N0(v)  production  for  v'<7.  In  contrast  to  Reaction  (3),  the  reaction  of  N(2d) 
with  O2  appears  to  be  strongly  channeled  into  a  specific  product  config¬ 
uration;  this  will  be  discussed  further  below.  The  corrected  product  state 
distribution  for  Reaction  (2),  inferred  from  Figure  10  and  plotted  in 
Figure  9,  represents  a  lower  bound  for  v'<P;  the  upper  bound  is  that  of 
Figure  5a  (no  "hot"  channel  correction).  The  fact  that  the  lower  bound  gives 
linear  surprisal  for  v'=l-7  and  is  consistent  with  a  reasonable  "hot"  channel 
correction  gives  strong  support  for  the  choice  of  this  distribution  (Figure  9) 
as  the  actual  nascent  vibrational  state  distribution  for  Reaction  (2). 

The  observed  ratio  of  rotationally  thermal  to  rotationally  excited 
N0(v'>0)  is  approximately  3.  Correcting  for  rotational  thermalization  effects 
and  using  surprisal  extrapolations  to  higher  vibrational  levels  (see  below), 
the  apparent  ratio  of  N(2d)  to  N(2p)  contributions  could  be  as  small  as  ~  0.5. 
However,  as  noted  above,  our  determination  of  the  number  densities  of  the 
high-J  states  may  be  uncertain  owing  to  effects  of  curvature  in  the  dipole 
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moment  function  on  the  scaling  of  the  rotational  line  strengths  to  high  J. 

Thus  the  ratio  of  production  rates  for  Reactions  (2)  and  (3)  cannot  be 
precisely  determined  from  the  present  data,  but  may  be  of  order  unity  or  less 
for  temperatures  of  80  to  100  K.  This  is  in  significant  contrast  to  our 
expectations  based  upon  room  temperature  kinetic  data.  As  shown  in  Table  4, 
for  room  temperature  discharges  similar  to  those  used  in  COCHISE,  we  typically 
observe  number  density  ratios  of  [N(2d) ]/[N(2p) ]  ~  3  at  the  discharge  exit 
using  atomic  resonance  absorption  and  fluorescence  techniques. Taken 
together  with  the  rate  coefficient  ratio  of  k2/k3  ~  2  at  room 
temperature, this  would  indicate  a  ratio  of  N0(v)  production  rates  of 
about  6.  There  are  several  possible  implications  from  this  discrepancy: 

1)  our  underestimate  of  the  high-J  line  strengths  may  be  far  more  severe  than 
expected;  2)  there  may  be  a  significant  unreactive  component  (i.e.,  direct 
energy  transfer  to  O2  without  NO  formation)  for  Reaction  2  at  low  tempera¬ 
ture;^^  and  3)  the  rate  coef f ic'ents  for  Reactions  (2)  and  (3)  may  follow 
opposite  temperature  dependencies,  i.e.,  k3  may  follow  a  negative  temperature 
dependence,  such  that  k2/k3  becomes  considerably  smaller  at  low  temperatures. 
All  of  these  possibilities  have  significant  implications  for  auroral 
chemistry,  and  we  are  continuing  to  pursue  them  in  the  laboratory. 

The  final  vibrational  distribution  shown  in  Figure  9  for  Reaction  (2)  is 
in  significant  disaccord  with  the  results  of  Green  et  al.^^  inferred  from 
electron  beam  experiments  at  higher  pressures.  Specifically,  the  corrected 
relative  populations  from  Figure  9  for  v'=l-3  are  much  lower  than  the  initial 
relative  populations  for  these  levels  inferred  from  a  time-resolved  kinetic 
analysis  of  N0(v)  distributions  created  and  partially  relaxed  in  a  pulsed 
electron  beam  experiment  irradiating  N2/O2  mixtures.  Green  et  al.^^  con¬ 
vincingly  eliminated  N(2p)  as  a  precursor  in  their  experiment  by  noting  that 
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the  NO(v)  precursor  was  quenched  by  N2  (with  a  rate  coefficient  appropriate 
for  N(2d)),  whereas  N(2p)  is  not.^  Under  electron  beam  conditions,  much  of 
the  N*  production  is  due  to  dissociative  recombination  of  N2‘*'  and  NO'*’,  which 
would  tend  to  favor  N(2d)  formation,  so  that  the  relatively  small  N(2p)  con¬ 
tribution  would  be  difficult  to  identify.  This  is  in  contrast  to  microwave 
disharge  conditions,  where  there  is  much  less  ionization  and  N(2D,2p)  are 
excited  by  direct  electron  impact.  Thus,  Green  et  al's  large  populations  at 
low  v'  are  probably  not  due  to  Reaction  (3),  unless  there  is  some  unforeseen 
excitation  or  deexcitation  mechanism  which  would  cause  N(2p)  to  follow  an 
unexpected  kinetic  behavior  under  their  experimental  conditions.  However, 
their  kinetic  analysis  relied  heavily  on  the  assumption  of  single  quantum 
collisional  deactivation  of  NO(v).  Vhile  this  appears  to  be  a  reasonable 
approximation  for  describing  the  loss  rate  of  a  given  vibrational  level, it 
may  lead  to  a  significant  underestimate  of  the  production  rate  for  the  lowest 
levels  by  multiquantum  deactivation  of  higher  levels  by  atomic  oxygen.  In 
classical  trajectory  studies  of  0  +  02(v)  collisions,  we  have  observed  a  simi¬ 
lar  effect,  with  efficient  multiquantum  (Av  =  2-4)  deactivation  of  02(v) 
resulting  in  rapid  population  of  lower  vibrational  levels  from  initially 
excited  higher  states. Thus  the  discrepancy  between  the  present  results  for 
Reaction  (2)  and  those  of  Green  et  al.^^  for  the  lower  vibrational  levels  of 
NO(v)  may  be  due  to  the  effects  of  multiquantum  0  +  NO(v)  deactivation  in 
their  experiments.  In  addition,  the  rapid  quenching  of  N(2p)  by  0^^  may 
compete  very  effectively  with  Reaction  (3)  for  their  electron  bombardment 
conditions.  Alternatively,  our  upper-bound  correction  for  rotational  therma- 
llzation  may  be  too  large  (i.e.,  rotational  thermalization  is  slower  than  one 
might  expect  from  the  data  of  Ref.  42),  and  the  product  distribution  of 
Reaction  (2b)  simply  does  not  follow  linear  surprisal  behavior. 
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In  an  attempt  to  discriminate  between  N(2d)  and  N(2p)  effects  on  the 
observed  NO(v)  distributions,  we  obtained  data  from  five  spectral  scans  over  a 
range  of  a  factor  of  ~  15  in  initial  N2  mole  fraction  in  the  discharge  gas 
mixture,  as  described  above.  For  this  range  of  conditions,  the  relative 
production  rates  for  NC^D)  and  N(2p)  in  the  discharges  remain  constant  to 
within  lOX,  but  the  N2  quenching  contribution  to  the  loss  rate  of  N(2d)  in  the 
sidearm  will  provide  some  variation  in  the  N(2D)/N(2p)  ratio  entering  the 
reaction  zone.  Unfortunately,  due  to  the  short  residence  time  in  the  dis¬ 
charge  sidearms,  this  amount  of  variation  in  [N2]  will  provide  only  about  a 
30  percent  variation  in  this  ratio,  with  the  ratio  decreasing  with  increasing 
N2  flow  rate.  In  the  experiments,  we  observe  the  ratio  of  rotationally  thermal 
to  rotationally  "hot"  (N0(v))  (summed  over  all  observed  v  levels)  to  decrease 
systematically  by  a  factor  of  1.3  with  increasing  N2  flow  rate,  as  expected 
for  the  production  of  thermal  N0(v)  by  Reaction  (2)  and  of  rotationally  "hot" 
N0(v)  by  Reaction  (3).  This  result  is  then  consistent  with  the  conclusions 
from  the  analysis  of  the  population  distributions.  However,  there  is  suffi¬ 
cient  scatter  in  the  data  that  an  average  ratio  for  all  the  scans  gives  a  ±  2a 
uncertainty  of  26  percent,  so  the  observed  trend  lies  only  slightly  outside 
the  scatter  of  the  data. 

DISCUSSION 

The  surprisal  parameters  determined  from  the  above  data  analysis  provide 
a  basis  for  extrapolation  of  the  results  to  other  vibrational  levels  for  which 
direct  observations  cannot  be  obtained.  For  the  N(2d)  reaction,  we  find  that 
the  branch  for  0(^D)  production  (Reaction  (2b))  can  account  for  all  the  pro¬ 
duction  of  N0(v'al-7),  with  a  surprisal  parameter  X  =  -6.25±0.14  (99X  confi¬ 
dence  level),  while  the  remaining  production  of  N0(v'>8)  is  due  to  the  0(3p) 
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branch  (Reaction  (2a)),  with  X  »  9.16±1.19.  The  relative  initial  populations 
of  NO(v)  determined  from  these  surprlsal  fits  are  shown  in  Table  5.  Summing 
the  initial  populations  for  each  branch  over  all  v  gives  the  branching  ratio 
k.2b/k2  =  0.76  at  -  100  K.  Similarly,  for  the  N(2p)  reaction  we  infer  (at 
most)  equal  branching  ratios  for  formation  of  0(^S),  0(^D),  and  0(3p),  with 
essentially  statistical  product  state  distributions  (see  Table  5). 

The  aeronomic  implications  of  the  present  conclusions  for  the  auroral 
upper  atmosphere  are  quite  significant.  The  production  of  0(^D,  ^S)  from 
these  reactions  has  an  impact  on  interpretations  and  modeling  of  ultraviolet 
and  visible  airglow  and  auroral  emissions  from  these  states  anct  from  other 
species,  such  as  N2(a3Eu‘'’)  and  metastable  states  of  O2,  which  are  believed  to 
be  linked  to  their  production. Our  interpretation  that  v=l-2  are 
inefficiently  populated  in  Reaction  (2)  is  a  significant  revision  of  previous 
conclusions, 12,31  ^^d  will  have  an  important  bearing  on  the  interpretation  of 
observed  auroral  N0(v)  distributions,  which  often  show  hitherto  unexplained 
"inversions"  in  the  population  of  v=2  (v=l  is  known  to  be  strongly  populated 
by  non-chemiluminescent  processes) These  results  also  impact  the 
current  views  on  the  roles  of  0-atom  deactivation  of  higher  vibrational  levels 
and  of  Reaction  (3)  as  sources  of  the  observed  populations  in  lower  vibra¬ 
tional  levels  under  certain  conditions,  suggesting  that  these  processes  are 
more  important  than  previously  thought.  Furthermore,  the  high  degree  of 
rotational  excitation  provided  by  the  N(2p)  reaction  will  have  a  dramatic 
effect  on  the  high-altitude  infrared  spectrum  of  auroral  N0(v)  emission,  as 
evidenced  by  the  R-branch  band  heads  for  v=l-3  observed  in  a  recent 
rocketborne  interferometer  flight. 22  The  presence  of  these  R-branch  band 
heads  to  the  blue  of  the  main  band  means  that  the  corresponding  P-branch 
distributions  must  appear  to  the  red.  Our  calculations  indicate  that  the 
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unrelaxed  P-branch  contributions  from  the  observed  eight  rotationally  excited 
vibrational  states  can  provide  more  intensity  in  the  6  to  8  ym  range  than  the 
emission  from  high  vibrational  levels  (for  any  rotational  distribution) 
produced  from  Reactions  (2a)  or  (3a).  Our  observations  further  indicate  that 
the  very  high  J  levels  responsible  for  these  contributions  are  deactivated 
relatively  slowly  in  Ar/02/N2  mixtures.  The  extent  of  collisional  coupling  in 
the  COCHISE  experiments  is  similar  to  that  expected  for  NO(v)  at  altitudes  of 
110  to  120  km  (assuming  the  primary  loss  for  N0(v)  is  by  radiative  cascade). 
Thus  sensitive,  high-resolution  measurements  of  the  N0(v,J)  spectrum  above 
~  100  km  in  a  strong  aurora  could  ultimately  provide  a  means  of  estimating  the 
N(2d)/N(2p)  ratio  and  hence  the  energetics  of  the  auroral  event. 

Our  approach  to  describing  the  rotationally  "hot"  distributions  deserves 
some  comment.  Due  to  limitations  of  spectral  resolution  and  sensitivity,  we 
are  unable  to  directly  observe  any  part  of  the  rotational  distributions  other 
than  the  distinctive  R-branch  band  heads.  Thus  there  is  no  clear  indication 
that  the  low  to  intermediate  J-levels  are  populated,  either  in  the  initial 
reaction  or  on  the  time  scale  of  the  measurements.  It  is  conceivable  that  the 
initial  rotational  populations  are  statistical,  or  that  the  initial  population 
of  high-J  is  rapidly  randomized  by  collisions  to  give  an  apparently 
statistical  distribution.  Such  a  distribution  would  have  the  behavior: 

PvO(J)  «  (2J  +  1){(E-Ev)  -  Er}l/2 

where  E^,  Ey,  and  E  are  the  rotational,  vibrational,  and  total  reaction 
energies,  respectively.  Thus  the  rotational  "temperature"  would  depend  on  the 
ribrational  energy  of  the  emitting  state,  with  higher  vibrational  levels 
exhibiting  "colder"  rotational  distributions.  Ve  cannot  clearly  discern  this 
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effect  in  the  NO(v,J)  spectra,  except  for  the  disappearance  of  obvious  band 
heads  above  v'=8.  However,  in  similar  studies  of  CO(v,J)  excitation  by  energy 
transfer  from  active  nitrogen, ^0  we  are  able  to  fit  the  observed  rotational 
distributions  better  with  the  statistical  formula  than  with  a  Boltzmann 
(v-independent)  description.  Planned  improvements  in  the  measurement 
sensitivity  for  N0(v,J)  may  permit  us  to  address  this  issue  more  accurately  in 
the  future. 

The  obser’'ed  product  state  distributions  imply  strongly  contrasting 
behavior  in  the  reaction  dynamics  of  the  N(2d)  and  N(2p)  reactions  with  02- 
For  Reaction  (2b),  the  high  degree  of  selectivity  for  the  highest  accessible 
vibrational  levels  is  characteristic  of  a  direct,  abstraction  reaction  with  a 
repulsive  approach  configuration  and  a  very  short-lived  collision  complex.  It 
would  appear  that  successful  reactive  encounters  require  maximum  vibrational 
motion  along  the  product  N — 0  coordinate,  perhaps  as  a  result  of  a  barrier  in 
the  N — O2  entrance  channel.  The  behavior  for  Reaction  (2a)  is  somewhat 
surprising,  since  the  reaction  appears  to  occur  only  when  Reaction  (2b)  is 
thermochemically  forbidden.  This  behavior  suggests  that  all  reactive  N — O2 
..ncounters  proceed  initially  along  the  surface  leading  to  0(^D)  formation, 
but  some,  with  too  much  N — O2  motion  in  the  product  channel  (insufficient 
motion  along  the  NO — 0  coordinate)  to  form  0(^D)  undergo  enough  motion  in  the 
collision  complex  to  dissociate  to  highly  vibrationally  excited  NO  +  0(^P)  via 
transition  to  one  of  the  other  ^A'  or  ^A"  surfaces  that  are  available 
(cf.  Table  1).  Thus  there  may  be  no  direct  path,  or  perhaps  a  large  energy 
barrier,  to  oO?)  formation  in  Reaction  (2).  If  the  occurrence  of  this 
reaction  branch  indeed  depends  upon  the  lifetime  of  the  collision  complex, 
then  its  contribution  may  diminish  substantially  at  elevated  collision 
energies  (or  higher  temperatures).  The  direct  reaction.  Reaction  (2b),  ought 


652 


to  follow  at  least  a  temperature  dependence,  or  possibly  one  somewhat 

stronger  due  to  the  presence  of  an  energy  barrier  in  the  entrance  channel. 

For  comparison,  Slanger  et  al.^®  observed  essentially  no  activation  energy  for 
quenching  of  Nf^D)  by  O2  between  240  and  370  K;  however,  the  temperature 
dependence  for  NO  formation  has  never  been  determined.  Clearly,  the  effects 
of  temperature  dependence  on  this  reaction  need  to  be  studied  in  more  detail, 
as  the  present  data  base  lends  somewhat  ambiguous  support  to  the  possibility 
of  an  unreactive  energy  transfer  channel  at  lower  temperatures.  This 
possibility  has  been  suggested  in  interpretations  of  aeronomic  data  on 
emission  from  excited  states  of  02,^^  albeit  on  controversial  grounds. 

In  contrast  to  Nf^D),  N(2p)  appears  to  react  with  O2  through  a  long-lived 
collision  complex,  resulting  in  complete  randomization  of  the  reaction  energy. 
The  appearance  of  high  rotational  excitation  is  characteristic  of  a  highly 
attractive  approach  configuration.  The  collisions  must  satisfy  the  principles 
of  conservation  of  angular  momentum  and  energy: 

(2yEi)l/2b.  =  (2MEf)l/2b^  +  jh/2ii 

Ej  +  AE  =  Ef  +  Ey  +  Bj2 

where  u  is  the  reduced  mass  of  N — O2,  Ej  and  Ef  are  the  initial  and  final 
center-of-mass  collision  energies,  bj  and  bf  are  the  initial  and  final  impact 
parameters,  AE  is  the  exoergicity  of  the  reaction,  Ey  is  the  product 
vibrational  energy,  and  Bj2  approximates  the  product  rotational  energy  after 
the  collision.  It  is  straightforward  to  show^^  that  large  values  of  J  can  be 
achieved  through  highly  attractive  collisions  in  which  large  initial  impact 
parameters  are  converted  to  small  final  impact  parameters,  with  the  concurrent 
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result  that  the  final  center-of-mass  kinetic  energy  is  small.  By  analogy  to 
the  well-studied  reaction  dynamics  for  0(^D)  with  and  HCl,^^  this  type  of 

collision  can  occur  through  a  C2v  Insertion  approach,  leading  to  a  highly 
attractive  ONO  collision  complex.  Such  a  complex  would  probably  undergo 
enough  vibrational  periods  to  completely  randomize  the  internal  energy,  giving 
rise  to  statistical  product  state  distributions.  With  so  much  Lissajous 
motion  in  the  complex,  the  adiabatic  correlation  rules  given  by  Shuler^^  would 
have  little  relevance, so  there  is  no  basis  for  the  previously 
suggested^’ prohibition  of  Reaction  (3a).  The  reaction  rate  coefficient  may 
increase, substantially  with  decreasing  temperature,  owing  to  both  the  dynamics 
of  the  initial  insertion  approach  (favored  for  less  O2  rotational  motion)  and 
the  Increased  lifetime  of  the  collision  complex  at  lower  temperatures.  This 
trend,  together  with  those  suggested  above  for  Reaction  (2),  is  consistent 
with  our  observation  of  an  unexpectedly  large  k.3/k.2  ratio  at  ~  100  K. 

It  is  interesting  to  compare  our  results  for  Reaction  (3)  with  those  of 
Taherian  et  al.^^  for  N0(v,J)  observed  in  the  photolysis  of  NO2.  Their 
experiments,  conducted  at  157.6  nm,  are  nearly  isoenergetic  with  the 
N(2p)  +  O2  reaction,  and  may  probe  portions  of  the  same  ONO  surface.  They 
observed  excitation  of  high  rotational  states  comparable  to  those  reported 
here,  together  with  a  partial  vibrational  distribution  which  is  not 
Inconsistent  with  our  statistical  picture.  Most  interestingly,  they  have 
identified  the  excitation  of  v=21,  confirming  that  0(3p)  is  a  primary  product 
of  thr  photodissociation  process.  However,  the  authors  report  upper  bounds 
for  0(^0,  ^S)  formation  which  indicate  these  to  be  minor  product  channels  at 
best.  Such  a  result  for  Reaction  (3)  would  not  give  a  good  surprisal  fit  to 
our  data,  but  certainly  cannot  be  ruled  out  on  that  basis  alone.  Alterna¬ 
tively,  this  disparity  between  our  results  and  theirs  may  arise  from  the  much 
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larger  phase  space  sampled  in  reactive  collisions.  It  is  clear  that  more 
definitive  product  state  identifications  are  required  to  confirm  our 
preliminary  assignment  of  the  product  branching  ratios  of  Reaction  (3). 


CONCLUSIONS 

We  have  used  the  technique  of  low-pressure  infrared  chemiluminescence  to 
investigate  the  vibrational  and  rotational  excitation  of  NOCX^II)  formed  in  the 
reactions  of  N(2d,  ^p)  with  O2  near  100  K.  The  initial  vibrational  distri¬ 
butions  obtained  from  the  experiments  are  composites  from  several  reaction 
channels.  The  results  of  surprisal  analysis  indicate  that  N(^D)  reacts  with 
O2  via  a  direct  abstraction  mechanism,  probably  on  a  repulsive  potential 
surface,  to  give  a  product  vibrational  distribution  which  is  sharply  peaked  to 
high  vibrational  levels.  This  interpretation  also  implies  that  reaction  to 
form  0(^D)  comprises  76  percent  of  the  total  N0(v)  production  at  these  tem¬ 
peratures.  In  contrast,  the  results  indicate  that  N(2p)  reacts  with  O2 
through  a  long-lived  complex,  perhaps  via  an  insertion  mechanism  on  an  attrac¬ 
tive  potential  surface,  to  give  essentially  statistical  product  state  distri¬ 
butions  with  extensive  rotational  excitation  of  the  molecular  products.  The 
observed  N0(v)  distributions  are  consistent  with  approximately  equal  contri¬ 
butions  from  the  reaction  channels  forming  0(^5),  0(^0),  and  0(3p)  as  products. 

The  results  of  this  work  amplify  and  extend  the  information  obtained  from 
previous  studies  of  these  reactions,  and  provide  additional  insights  into  the 
kinetics  and  spectroscopy  of  vibrationally  excited  NO  in  the  auroral  upper 
atmosphere.  However,  the  results  of  the  present  investigations  raise  a  number 
of  issues  requiring  further  study,  including  the  temperature  dependencies  of 
the  various  reaction  pathways,  the  possibility  for  quenching  collisions  which 
do  not  produce  N0(x2n),  confirmation  of  0(1d,  ^S)  production,  and  quantitative 
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determination  of  the  rotational  line  intensities  at  high  J.  We  are  continuing 
to  address  these  issues  through  experiments  on  the  COCHISE  facility  as  well  as 
with  discharge-flow  reactors  using  resonance  fluorescence  and  laser-based 
diagnostics.  In  addition,  these  reaction  systems  lend  themselves  well  to 
theoretical  investigations  of  quantum  mechanical  potential  energy  surfaces  and 
semiclassical  collision  dynamics. 
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Table  1.  Pathways  for  NO  Formation  from  N(2d,  ^P)  +  O2 


Reaction 

Adiabatic 

Surfaces 

— 

^°298  K 
kcal/mole 

— 

Maximum 

V  for  NO 

2  N(2Du)  +  02(x3i:g-) 

a  ^  N0(x2n)  +  0(3pg) 

22a '  +  22a'»  ^A'  +  ^A" 

-  86.80 

18 

b  N0(x2n)  +  O(lDg) 

2a' 

-  41.43 

8 

3  N(2Pu)  +  02(x3Eg-) 

a  N0(x2n)  +  0<3pg) 

None 

-114.28 

26 

b  -»  N0(x2n)  +  O(lDg) 

2a» 

-  68.91 

14 

c  -►  N0(x2n)  +  O(lSg) 

2a« 

-  17.67 

3 

d  -»  NO(a^n)  +  0(3pg) 


2a'  +  ^A'  +  2^A" 


Table  2.  Spectroscopic  Constants  for  NOCX^Il)®*^ 


1 

1903.937 

«eXe 

13.970 

-  3.7  X  10-3 

(OfeZe 

5.0  X  10-5 

Be 

1.70488 

“e 

0.017554 

De 

5.36  X  10-6 

0e 

5.0  X  10-8 

Ae 

123:28 

Xe 

0.256 

^Values  are  from  Ref. 

24  except  for  pg  (see  text). 

^Units  are  appropriate  for  calculation  of  transition 

frequencies  in  cm"^. 

Table  3.  Band-Integrated  Einstein  Coefficients 
Emission^ 


V' 

Ay  ^  y  **  f  S  ^ 

1 

13.4 

2 

25.3 

3 

36.2 

4 

45.8 

5 

54.2 

6 

60.1 

7 

66.3 

8 

71.2 

9 

75.0 

10 

77.7 

11 

79.4 

12 

80.0 

13 

79.5 

14 

78.2 

15 

76.4 

16 

74.1 

17 

71.3 

18 

68.1 

19 

64.3 

20 

60.0 

^Determined  from  data  of  Ref.  27  and  scaled  according  to 
predictions  of  Ref.  28  (see  text). 
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Table  4.  Estimated  Number  Densities  of  Energetic 
Discharge  Effluents^ 


Species 

Number  Density, 

molecules  cm“^ 

Discharge  Exit 

Reaction  Zone 

N(^S) 

2  X  IQI^ 

7  X  loll 

N(2d) 

(1-4)  X  10^2 

0.3-1  X  lOlO 

N(2p) 

(0.5-1)  X  10^2 

1-3  X  109 

N2(xllg+,  v=l-7) 

3  X  10^5 

1  X  10l3 

N2(a3V) 

1  X  10^2 

3  X  109 

N2(v\,v1-2) 

3  X  lOlO 

3  X  108  b 

N2<a'lV) 

1  X  I0I2 

3  X  109 

N2(v1Au,v=0-2) 

6  X  10^0 

1  X  108  b 

^Microwave  discharge  conditions  given  in  the  text, 
approximately  10%  N2  in  Ar. 

^’Based  on  direct  observations  reported  in  Ref.  19. 
Higher  vibrational  levels  are  affected  by 
radiative  losses  during  transit  from  discharge 
exit  to  reaction  zone. 
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Table  5.  Surprisal  Predictions  of  Relative  Initial  Populations 
of  NO(v)  in  the  Reactions  of  N(2d,  2p)  vith  O2 


®k2a/lc2b  »  0.285  *’k3a''*^3b  “  •^3c/*'^3b  =  ^ 

«^l.36(-02)  signifies  1.36  x  10-2 
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Figure  1 


Figure  2 


Figure  3 


Figure  4. 


Figure  5. 


Diagram  of  the  COCHISE  reaction  chamber.  The  physical  dimensions 
of  the  cell  are  0.6  m  in  length  and  0.4  m  in  diameter. 

Spectrum  of  fluorescence  observed  in  the  interaction  of  O2  with 
discharged  N2/Ar  (12%  N2).  Spectral  resolution  is  0.007  pm  FVHM. 
Gas  temperature  is  80  K. 

Comparison  of  computed  and  observed  N0(v)  spectra.  Discharge: 

12%  N2  in  Ar;  Counterflov:  O2;  Gas  temperature:  80  K;  spectral 
resolution:  0.013  pm  FVHM;  corrected  for  background  emission 
from  the  discharges,  (a)  comparison  of  calculated  (heavy  line) 
and  observed  (light  line)  spectra,  (b)  Contribution  to  computed 
spectrum  from  rotationally  "cold”  component,  Trqt  =  ^^0  K, 
v=»l-14.  (c)  Contribution  to  computed  spectrum  from  rotationally 

"hot"  component,  Trqt  *  10,000  K,  v*l-8.  The  sharp  features 
between  4.9  and  5.6  pm  are  the  R-branch  band  heads  for  V3I-8. 

Vibrational  band  intensities  determined  from  spectral  fit 
in  Figure  3 

Average  vibrational  population  distributions  for 
(a)  rotationally  thermal  and  (b)  rotationally  "hot" 
portions  of  the  NO  chemiluminescence 
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Figure  6. 


Figure 


Figure 


Figure 


Figure 


Computed  prior  distributions  for  NO(v)  formed  in  N*  +  O2 
reactions,  (a)  N(2d)  +  02:®»  0(^P)  product;H, 

0(^0)  product;  (b)  N(2p)  +  02:0,  0(3p)  product;D, 

0(^0)  product;  7,  0(^S)  product 

7.  Comparison  of  rotationally  "hot"  populations  for  v=l-8  and 
rotationally  "thermal"  populations  for  v=12-14  to  composite 
prior  distribution  for  the  N(2p)  +  O2  reaction  (curve) 

8.  Surprisal  plots  for  N0(v)  formed  from  two  branches  of  the 
N(2d)  +  O2  reaction,  using  the  population  data  from  Figure  5a. 
The  lines  correspond  to  surprisal  parameters  X  of  -6.1  and  7.4 
for  the  0(1d)  and  0(3p)  branches,  respectively. 

9.  Rotationally  thermal  populations,  corrected  for  maximum 
thermalized  high-J  component.  The  dashed  curve  represents 
populations  determined  from  the  surprisal  analysis  of  Figure  10. 

10.  Surprisal  plots  and  weighted  least  squares  fits  for  N0(v)  formed 
from  two  branches  of  the  N(^D)  +  O2  reaction:  data  corrected 
for  maximum  thermalized  high-J  component.  Error  bars  have  been 
omitted  from  some  of  the  data  points  for  clarity.  The  least 
squares  lines  correspond  to  surprisal  parameters  X  of  -6.25  and 
9.16  for  the  0(^D)  and  0{3p)  branches,  respectively. 
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ABSTRACT 


CO  fundamental  vibration-rotation  band  emission  resulting  from  the  inter¬ 
action  of  discharged  nitrogen  with  carbon  monoxide  has  been  studied  at  low 
pressure  ( ~3  mt)  in  a  cryogenic  apparatus.  The  spectra  exhibit  bimodal  rota¬ 
tional  distributions;  we  have  identified  fourteen  vibrational  levels  of  a 
thermalized  component  and  eight  vibrational  levels  from  a  rotationally  excited 
component.  The  rotationally  excited  emission  is  adequately  reproduced  by  a 
statistical  distribution,  E  =  3.7  eV,  which  provides  sufficient  Dopulation  in 
the  region  of  the  Fortrat  reversal  (J  ~  90)  to  account  for  the  observed 
R-braiiCh  band  head  formation.  The  excitation  reactions  are: 


N2(v)  +  CO  ^  N2(v-1)  +  C0(v=l)  (a) 

N2(v-2)  +  C0(v=2)  (b) 

N2(a'li:u-)  +  CO  ->  N2(X,v)  +  CO(Aln)  C0(v<9)  +  hv  (c) 

-♦  N2(X,v)  +  C0(v<l4,  J)  (d) 


Single  and  two-quantum  transfer  from  N2(v),  reactions  (a)  and  (b),  contributes 
to  the  thermalized  C0(v=l,2)  populations.  Radiative  cascade  from  C0(A)  pro¬ 
duced  by  quenching  of  N2(a'),  reaction  (c),  also  contributes  to  the  thermalized 
populations.  Based  on  kinetic  and  energetic  arguments  we  have  determined 
another  branch  of  N2(a')  quenching  to  be  responsible  for  the  rotationally 
excited  component. 

Surprisal  analysis  of  the  rotational  excited  component  vibrational  distri¬ 
bution  indicates  two  dynamic  mechanisms.  We  have  modeled  this  distribution 
with  equal  contributions  from  a  direct  (v<4)  and  a  statistical  (all  v)  process. 
The  vibrational  distributions  of  the  rotationally  thermal  component  are  repro¬ 
duced  by  CO(A^n)  radiative  cascade  and  a  contribution  from  reaction  (d). 


INTRODUCTION 


The  dynamics  of  energy  disposal  in  quenching/reaction  of  excited  state 
species  may  be  determined  from  the  nascent  electronic-vibrational-rotational- 
translational  distributions  of  the  product  molecules.  The  observed  product 
state  distributions  reflect  the  details  of  the  interaction  mechanisms  and 
provide  insight  into  the  nature  of  the  potential  energy  surfaces.  Ve  have 
previously  reported  observations  of  infrared  chemiluminescence  from  vibra- 
tionally  and  rotationally  excited  NO  generated  from  the  reaction  of  metastable 
atomic  nitrogen  with  molecular  oxygen.^  Surprisal  analysis  of  the  data  sup¬ 
ports  two  mechanisms:  a  direct  abstraction  reaction  (N^D)  -i-  O2  producing 
rotationally  thermal  but  vibrationally  excited  NO  and  an  insertion  reaction 
N(2p)  +  O2  giving  essentially  statistical  vibrational  and  rotational  product 
state  distributions. 

This  papei  reports  observations  of  highly  rovibrationally  excited  CO 
formed  by  energy  transfer  between  carbon  monoxide  and  nitrogen  metastables. 

The  observed  CO  fundamental  spectra  contain  rotationally  thermalized  emission 
for  v<14  and  rotationally  excited  bands  for  v^.  The  rotationally  excited 
component  appears  as  R-branch  band  heads  that  require  significant  populations 
in  high  rotational  levels  near  the  vertex  of  the  Fortrat  parabola 
(Jvertex  ~  energetics  of  both  the  thermal  and  rotationally  excited 

components  are  attributable  to  an  energy  transfer  process  which  deposits  up  to 
-3.7  eV  into  CO  rovibrational  states. 

Modest  energy  depositions  (typically  <1  eV)  in  CO  internal  states  have 
been  observed  in  photochemical  and  abstraction  reactions.  In  general,  large 
fractions  of  the  energy  above  threshold  are  manifested  in  product  internal 
states.  For  example,  photolysis  of  OCS  at  157  nm^  of  acetone^  at  193  nm,  and 
H2C0^»5  produces  vibrationally  excited  CO  with  rotational  excitations  up  to 
0.9  eV.  Studies  of  energy  partitioning  in  CO  from  the  reaction  of  hot  H  atoms 
with  CO2  indicates  nearly  1.0  eV  in  CO  internal  states^  with  the  rotational 
distribution  in  the  v>0  level  following  a  statistical  model. 
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Greater  energy  depositions  are  observed  for  energy  transfer  reactions 
since  these  interactions  have  little  or  no  threshold  energy.  Quenching  by 
0(21d),7  Na(32p),8,9  i(52p^^2),10  and  Br(42p3^/2) , 10  and  HgCe^P^  and  b^PQ)!! 
have  been  observed  to  produce  vibrationally  excited  CO.  Rotational  excitations 
of  <0.8  eV  have  been  reported  for  the  Na(32p)  +  CO  quenching  reaction^  which 
accounts  for  a  large  fraction  of  the  total  2.1  eV  exoergicity. 

The  excitation  process  observed  here  is  energy  transfer  from  a  metastable 
state  of  nitrogen.  Energy  transfer  to  CO  from  discharged  nitrogen  has  been 
examined  closely,  particularly  with  respect  to  the  N2-CO  laser. 12, 13  High  CO 
vibrational  excitation  is  produced  from  near-resonant  N2(v)  energy  transfer  and 
subsequent  C0(v)  collisional  up-pumping, 1^» 1^  reactions  (la  and  b). 

N2(v)  +  CO(v')  -♦  N2(v-1)  +  C0(v'+1)  (la) 

C0(v)  +  CO(v')  ->C0(v'-l)  +  C0(v'+1)  (lb) 

Although  high  CO  vibrational  levels  are  produced  in  this  manner,  no  rotational 
excitation  is  observed,  which  is  consistent  with  the  small  energy  defect  for 
the  near-resonant  process.  The  energy  defect  from  N2(aln)  quenching, 18 
reaction  (2)  has  been  determined 

N2(aln)  +  CO  N2(X,v)  +  C0(Aln,  v' ,  J')  (2) 

to  be  largely  manifested  in  rotation  but  the  total  energy  defect  is  small 
(<1000  cm“l),  much  less  than  the  rotational  excitations  observed  here.  The 
quenching  of  N2(A8Zy+)  by  C0l2  produces  C0(a8n)  with  no  reported  observations 
of  a  C0(v,J)  product  channel. 

The  energy  transfer  process  reported  here  represents  a  previously 
unobserved  quenching  reaction  of  CO  with  metastable  nitrogen.  We  will  show 
that  consideration  of  kinetic  and  energetic  constraints  positively  identifies 
N2(a'lEy“)  as  the  responsible  agent.  We  will  also  present  and  discuss  the 
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results  of  surprisal  analysis  of  the  observed  vibrational  distributions  which 
indicates  the  rotationally  excited  component  to  be  formed  by  two  distinct 
mechanisms. 


EXPERIMENTS 

These  experiments  were  performed  in  the  COCHISE  (COld  CHemiexcitation 
Infrared  Stimulation  Experiment)  cryogenic  discharge  afterglow  apparatus  which 
is  described  in  detail  elsewhere.^®  Excitation  of  nitrogen  is  achieved  with 
four  parallel  microwave  discharges  (2450  MHz,  50  W)  of  flowing  N2/Ar  mixtures 
at  -1  torr  total  pressure.  A  diagram  of  the  reaction  chamber  is  shown  in 
Figure  1.  After  exiting  the  discharge  tubes  the  gas  expands  into  a  low  pres¬ 
sure  (~3  mt),  cryogenlcally  pumped  chamber  ( ~20  K)  where  the  molecules  enter 
the  collimated  field  of  view  of  a  scanning  monochromator/infrared  detector 
assembly.  Residence  times  in  the  discharge  tubes  are  on  the  order  of  3  to 
5  ms;  an  average  time  of  flight  of  0.5  ±  0.1  ms  is  required  for  the  gases  to 
exit  the  discharge  tubes  and  enter  the  field  of  view.  Opposing  flows  of  argon/ 
carbon  monoxide  are  used  to  create  a  quasi-static  interaction  region  along  the 
centerline  of  the  field  of  view,  resulting  in  partial  rethermalization  of  the 
expansion  cooled  rotational  distributions.  The  gas  residence  time  in  the  field 
of  view  is  -0.3  ms.  Gaseous  helium  refrigerant  maintains  all  internal  tempera¬ 
tures  at  20  K,  excepting  the  gas  lines  and  optics  which  are  held  at  80  and 
40  K,  respectively. 

The  infrared  emissions  are  observed  by  a  cryogenic  0.5  m  Czerny-Turner 
monochromator  equipped  with  a  liquid-helium-cooled  arsenic-doped  silicon 
detector  and  a  grating  blazed  at  3  ym.  A  chopper  located  in  front  of  the 
monochromator  entrance  slit  modulates  the  signal  at  23  Hz.  Data  collection  is 
performed  with  a  computer-interfaced  lock-in  amplifier.  The  absolute  uncer¬ 
tainty  in  the  wavelengths  (due  to  monochromator  drive  error)  of  the  data  is 
±0.003  ym.  The  data  were  corrected  for  instrument  responsivity  using  blackbody 
calibration  spectra  taken  in  the  300  to  370  K  range.  The  uncertainty  in  the 
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accuracy  of  the  blackbody  temperature  is  +3  K,  which  results  in  a  relative 
error  of  +14  percent  for  4. 0/6.0  pm  intensity  ratios.  Spectra  were  taken  for 
N2/Ar  mixtures  with  N2  mole  fractions  of  0.005  to  0.12,  with  a  mass-balanced 
counterflow  of  CO/Ar  with  CO  mole  fractions  of  0.018  to  0.35.  The  data  were 
taken  typically  at  a  spectral  resolution  of  0.013  pm  (FWHM). 


RESULTS 

In  all,  24  emission  spectra  of  the  CO  fundamental  region  were  obtained  at 
various  N2  and  CO  mole  fraction  conditions.  The  conditions  are  similar  to 
those  in  which  nitric  oxide  fundamental  emission  was  examined.^  The  general 
features  of  the  emissions  are  relatively  invariant  with  CO  mole  fraction  but 
the  intensity  of  the  CO(l-O)  emission  exhibits  a  strong  mole  fraction 
dependence  as  shown  in  Figure  2. 

The  envelope  degrading  to  the  red  of  the  CO(l-O)  band  center  at  4.666  pm 
is  due  to  the  Av=l  progression  from  v<14.  The  five  sharp  red-degraded  features 
to  the  blue  of  the  CO(l-O)  band  center,  which  have  a  spacing  of  29.2  ±1.4  cm“^, 
cannot  be  attributed  to  CO  thermalized  vibrational  emission.  These  bands  have 
not  been  previously  observed  in  published  spectra  of  CO  fundamental  emission 
taken  at  higher  pressures  (>1  torr).^^"^!  low  pressure  with  the  COCHISE 
apparatus,  these  features  are  observed  under  all  conditions  that  produce  the  CO 
fundamental  emission.  They  have  been  observed  to  remain  unstructured  even  at 
the  highest  resolution  employed  (0.0067  pm).  Under  low  nitrogen  mole  fraction 
conditions,  in  which  the  CO(l-O)  emission  intensity  is  greatly  reduced,  eight 
bands  are  observed  with  three  progressing  into  the  thermalized  C0(Av=l)  enve¬ 
lope.  Considering  the  overlap  of  the  CO  thermalized  emission,  identification 
of  eight  such  features  represents  a  lower  limit. 

These  emissions  are  not  observed  in  the  absence  of  nitrogen  in  the  dis¬ 
charge  mixture.  Thus,  Argon  metastables  and  residual  ions  do  not  contribute. 


6R9 


We  have  determined  these  features  are  not  due  to  CN(v),  electronic  N2  transi¬ 
tions,  or  NCO.  The  spectral  shape  of  these  features,  sharply  peaked  and 
degraded  to  the  red,  is  similar  to  the  nitric  oxide  R-branch  band  heads 
identified  in  spectra  of  chemiluminescence  produced  from  the  reaction  of 
discharged  nitrogen  with  oxygen.^  CO  forms  band  heads  similarly,  at  suffi¬ 
ciently  high  rotational  excitation  band  heads  will  form  in  the  R-branches  and 
the  P-branches  will  extend  to  the  red  and  do  not  form  band  heads.  We  have  used 
a  spectral  generation  technique  to  predict  the  band  shapes  and  positions  of  the 
CO  R-branch  band  heads.  The  methodology  and  result  will  be  presented  in  detail 
in  the  following  section.  Using  this  technique  the  sharp  red-degraded  features 
have  been  positively  identified  as  CO  R-branch  band  heads.  These  features  are 
sufficiently  intense  that  we  have  been  able  to  determine  an  optimum  rotational 
distribution.  Additionally,  the  P-branches  of  the  rotationally  excited 
component  associated  with  the  R-branch  band  heads  are  identifiable  which 
verifies  both  the  spectral  assignment  and  confirms  the  appropriate  choice  of 
rotational  distribution. 

SPECTRAL  ANALYSIS 

The  data  have  been  analyzed  using  a  spectral  generation  linear  least 
squares  fitting  technique. A  computed  infinite  resolution  spectrum  is 
convolved  with  the  instrument  scan  function  (in  this  case,  a  symmetric  triangle 
with  full  width  at  half  maximum  as  the  spectral  resolution)  to  create  a  basis 
set  for  each  vibrational  transition.  The  basis  sets  are  then  fit  to  each 
experimental  spectrum  using  a  linear  least  squares  method  yielding  a  determina¬ 
tion  of  the  product  of  the  upper  state  density  and  the  spontaneous  emission 
coefficient  of  the  transition,  Ny>  Ayf^yn. 

The  spectroscopic  data  used  in  this  study  are  from  Huber  and  Herzberg.23 
The  data  were  sufficient  to  adequately  reproduce  thy  line  positions  of  both  the 
rotationally  thermalized  emission  features  and  the  R-branch  band  heads. 


690 


The  rotational  temperature  used  to  reproduce  the  CO  vibrational  progres¬ 
sion  to  the  red  of  4.6  pm  was  determined  empirically  by  reproduction  of  the 
observed  branch  structure.  The  optimum  temperature  was  determined  to  be  80  K 
indicating  this  emission  system  is  rotationally  thermalized.  The  rotational 
distribution  was  treated  by  a  simple  Boltzmann  expression,  so  that  band- 
integrated  vibrational  number  densities  and  transition  probabilities  were  used. 
The  transition  probabilities  were  calculated  using  the  dipole  moment  function 
of  Chackerian  et  al.^^  Fourteen  vibrational  levels  of  thermalized  CO  emission 
have  been  positively  identified  from  the  data.  This  corresponds  to  vibrational 
excitation  of  3.425  eV. 

The  band-averaged  transition  probabilities  used  for  the  rotationally 
excited  component  are  the  same  as  those  for  the  thermalized  CO  emission.  We 
believe  these  values  to  be  accurate  since  the  CO  ground  state  is  and 
therefore  not  subject  to  spin-uncoupling  at  high  J'  as  discussed  for  NO.  In 
addition,  the  dipole  moment  function  of  Chackerian  et  al.^^  extends  to  v'=:40, 
well  above  the  internuclear  separation  of  the  highest  V'  and  J'  levels 
encountered  in  this  study.  Thus,  Honl-London  scaling  of  the  band-averaged 
transition  probabilities  is  sufficient  to  accurately  determine  the  populations 
of  the  rotationally  excited  component. 

The  5  cm“l  resolution  (at  5.0  pm)  of  the  spectral  data  is  insufficient  for 
rotational  resolution  of  the  R-branch  band  heads  so  the  chosen  rotational  dis¬ 
tributions  are  those  which  best  reproduce  the  band  shapes  and  peak  positions. 

To  adequately  fit  the  R-branch  band  head  features  by  a  Boltzmann  distribution 
requires  temperature  of  approximately  20,000  K.  Although  the  fits  to  the  lower 
vibrational  levels  are  adequate,  higher  vibrational  levels  are  poorly  fit  and 
are  better  described  by  lower  Boltzmann  rotational  temperatures.  This  is 
evidence  for  an  anti-correlation  between  rotational  and  vibrational  excitation; 
i-e.,  the  lowest  vibrational  levels  contain  the  most  rotational  excitation. 

Such  anticorrelations  have  been  observed  in  photolysis  experiments,  abstraction 
reactions  and  energy  transfer  processes. In  these  instances,  the  observed 
rotational  distributions  are  sometimes  best  described  by  a  statistical  model. ^ 
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Such  a  model  distributes  the  population  statistically  over  the  accessible 
states.  The  model  which  incorporates  the  observed  anticorrelation  between 
vibrational  and  rotational  excitation  is  given  by, 

PvO(J)  a  (2J  +  1){(Et  -  Ev)  -  Ej)1/2  (3) 

where  E-]-  is  the  total  energy  available  for  product  states,  E^  is  the  vibra¬ 
tional  energy,  and  Ej  is  the  energy  of  the  particular  rotational  level. 

Figure  3  contrasts  the  relative  population  distributions  predicted  by  Boltzmann 
and  statistical  models.  The  statistical  model  contains  greater  relative 
population  in  the  higher  rotational  levels  at  the  expense  of  the  lower.  Thus, 
rotational  band  head  formation,  which  requires  significant  population  of 
rotational  levels  in  the  region  of  the  Fortrat  parabola  vertex  (Jvertex  =  ^2 
for  v=l),  is  readily  facilitated  by  a  statistical  distribution. 


Ve  have  examined  several  values  of  Ej  to  determine  which  produces  the  best 
fit  to  the  data,  i.e.,  reproduction  of  the  positions  and  spectral  shapes  of  the 
observed  eight  band  heads.  Values  of  E-p  below  3.0  eV  reproduce  the  lower 
vibrational  levels  well  but  the  predicted  band  heads  at  higher  vibrational 
levels  are  broadened  and  red-shifted.  The  responsible  mechanism  may  be  seen 
from  Figure  3.  At  higher  Ey,  and  lower  Ej-,  the  relative  population  of  the 
rotational  levels  near  the  region  of  the  reversal  decreases  which  causes  the 
observed  effects.  Values  of  Ep-  greater  than  4.0  eV  predict  band  head  formation 
for  vibrational  levels  ten  or  greater.  We  prefer  a  value  of  E-p  between  these 
two  bounds.  Comparing  fits  using  several  values  within  this  range,  we  have 
determined  E-p  =  3.7  eV  to  be  optimum.  Values  within  +0.2  eV  of  the  optimum 
value  produce  adequate  fits  with  only  small  differences. 


Table  1  shows  the  energetics  of  CO  band  head  formation  for  Ep  of  3.7  eV. 
Shown  are  the  values  for  Ej.,  the  rotational  level  corresponding  to  the  vertex 
of  the  Fortrat  parabola,  the  rotational  energy  corresponding  to  this  value  of 
J,  and  the  maximum  rotational  level  allowable  from  Ej.. 
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Figure  4  shows  a  typical  fit  to  the  data  using  14  vibrational  levels  of 
rotationally  thermalized  emission  and  eight  vibrational  levels  with  a  statis¬ 
tical  rotational  distribution  corresponding  to  Ef  =  3.7  eV.  All  of  the 
principal  spectral  features  are  accurately  reproduced.  The  only  features  not 
entirely  reproduced  fail  within  the  5.7  to  6.5  ym  region.  Figure  5  shows  an 
enlarged  view  of  this  region  from  Figure  4.  The  fit  is  shown  by  the  heavy 
line.  The  spacing  of  the  bands  shown  in  Figure  4  is  ~6  cm”^.  These  features 
have  been  observed  in  all  spectra  containing  sufficient  intensity  in  this 
wavelength  region.  The  fit  shows  excellent  reproduction  of  the  spectral  shapes 
but  does  not  match  the  absolute  intensity.  Figure  4  shows  that  the  thermalized 
emission  does  not  contribute  to  this  spectral  region;  the  discrete  features  are 
reproduced  by  the  P-branches  of  the  rotationally  excited  component.  Inclusion 
of  v>8  of  the  rotationally  excited  component  improves  the  fits  in  this  wave¬ 
length  region.  Due  to  overlap  with  the  thermalized  CO  envelope,  however, 
unique  determination  of  the  populations  is  not  possible,  so  we  have  generally 
excluded  these  bands  from  the  fits.  We  describe  in  the  Discussion  Section  that 
unique  fits  may  be  obtained  by  fitting  the  rotational  excited  component  as  a 
single  variable  with  the  relative  vibrational  levels  fixed  in  ratios  from  a 
modeled  distribution.  By  comparing  fits  with  different  vibrational  contribu¬ 
tions  to  the  rotationally  excited  component,  we  have  determined  that  the 
majority  of  the  intensity  in  the  5.7  to  6.5  urn  region  arises  from  the  higher 
vibrational  levels  (v'=6-8).  The  relative  line  spacings  may  be  calculated 
(omitting  D^  correction)  from 


Av  =  P(J)  -  P(J=1)  =  (Bv^  +  Bv-O  -  (Bv'  -  Bv..)(2J+l)  .  (4) 

For  v'=6-8  this  wavelength  region  contains  P-branch  line  spacings  of  ~6  cm~^ 
for  rotational  levels  50-80.  The  N2(w3Au,  v'=1  B^Ilg,  v"=0)  emission  in  the 
6.3  to  6.6  urn  region,  which  we  have  previously  examined  and  reported, and 
signal  to  noise  considerations  impede  determination  of  the  full  extent  of  the 
P-branch  structure.  The  reproduction  of  these  features,  however,  by  fits  to 
the  rotationally  excited  component  confirms  the  identification  of  the  CO 
R-branch  band  heads. 
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All  24  spectra  have  been  fit  using  14  thermalized  CO  emission  bands  and 
eight  vibrational  levels  of  the  rotationally  excited  component,  as  shown  in 
Figure  4.  All  the  fits  show  similar  reproduction  of  the  principal  features 
although  many  do  not  accurately  reproduce  the  absolute  emission  intensity  of 
the  P-branch  envelope.  The  discrepancy  shown  in  Figure  4  is  typical.  Fits  to 
this  spectral  region  are  improved  by  inclusion  of  higher  vibrational  levels  of 
the  rotationally  excited  component.  Accurate  determination  of  any  contribution 
from  the  higher  levels  cannot  be  determined  due  to  the  overlap  of  these  band 
systems  with  the  thermalized  envelope.  Bounds  for  their  relative  populations 
will  be  discussed  in  the  next  section. 

The  statistical  rotational  distributions  employed  in  these  fits  provide  a 
better  reproduction  of  the  data,  both  in  the  R-branch  spectral  band  shapes  and 
the  absolute  intensity  of  the  P-branches,  than  do  Boltzmann  rotational  distri¬ 
butions.  Ve  consider  the  statistical  model  employed  here  to  be  a  more  accurate 
representation  of  the  true  rotational  distributions  but  it  cannot  be  considered 
to  be  a  unique  determination.  Signal  to  noise  considerations  and  band  overlap 
do  not  permit  unambiguous  identification  of  the  vibrational-level  dependent 
rotational  distributions.  The  true  distributions  may  deviate  from  a  purely 
statistical  model,  containing  gaussian  character  as  has  been  observed  in  other 
systems.  However,  the  model  employed  here  has  successfully  demonstrated  an 
anticorrelation  between  vibrational  and  rotational  excitation  and  permitted  a 
bound  of  ~3.7  eV  to  be  determined  for  the  rotational  excited  component.  Deter¬ 
mination  of  this  bound  and  its  similarity  to  the  vibrational  excitation  of  the 
thermalized  component,  3.425  eV,  suggests  that  these  two  components  arise  from 
the  same  excitation  process. 


KINETIC  INTERPRETATIONS 

The  kinetics  of  processes  occurring  in  the  COCHISE  reaction  chamber  have 
been  previously  described.^  Owing  to  the  low  number  densities  in  the  reaction 
zone,  the  short  residence  time  in  the  field  of  view  (0.3  ms),  and  the  long 
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radiative  lifetimes  for  the  infrared  chemiluminescent  processes  that  we 
observe,  reactions  and  quenching  (excluding  rotational)  of  the  excited  species 
created  in  the  reaction  zone  can  be  neglected.  Thus,  vibrational  quenching  and 
CO(v)  up-pumping,  reaction  (lb),  can  be  ignored.  The  kinetics  of  [CO(v,J)]  are 
therefore  in  steady-state  according  to  the  production  rate  and  the  lifetime  of 
the  excited  species  in  the  field  of  view, 

d[CO(v,J)]/dt  =  k[M*][C0]  -  Trl[C0(v,J)]  =  0  (5) 

where  M*  denotes  the  responsible  excited  species  created  in  the  microwave 
discharges,  k  is  the  excitation  rate  constant,  and  T  is  the  residence  time  in 
the  field  of  view,  0.3  ms. 

Figure  6  shows  the  absolute  populations  for  the  fit  in  Figure  4.  The 
bimodal  nature  of  the  rotational  distributions  is  apparent.  The  populations 
for  the  rotationally  excited  component  exceed  those  of  the  thermal  component  of 
all  vibrational  levels  except  val,2.  As  noted  in  Figure  2,  the  lowest  thermal- 
ized  CO  vibrational  levels  exhibit  an  N2  mole  fraction  dependence.  We  have 
determined  previously  that  only  N2(v)  exhibits  a  strong  variation  in  its 
discharge  production  rate  with  N2  mole  fraction  over  the  range  used  here. 

Thus,  relative  increases  in  low  CO  vibrational  population  at  higher  nitrogen 
mole  fraction  must  be  due  to  the  near-resonant  energy  transfer  from  N2(v), 
reaction  (la). 

As  illustrated  by  Figure  2  the  lowest  CO  vibrational  populations  exhibit  a 
marked  dependence  with  N2  mole  fraction.  Examination  of  the  determined  popula¬ 
tions  for  all  the  spectra  indicates  that  only  C0(v=l,2)  are  affected  by 
nitrogen  mole  fraction.  Since  C0(v)  up-pumping  cannot  occur  under  our  experi¬ 
mental  conditions,  the  excitation  process  must  be  single  and  two  quantum 
transfer. 


N2(v)  +  CO  ^  N2(v-1)  +  C0(v=l) 


(6a) 
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N2(v)  +  CO  N2(v-2)  +  CO(v=2) 


(6b) 


To  determine  the  relative  contribution  of  the  N2(v)  energy  transfer  process  to 
the  C0(v»l,2)  populations  we  have  plotted  the  ratios  of  these  populations  to 
[C0(v»3)]  as  a  function  of  nitrogen  mole  fraction.  Figure  7  shows  one  of  these 
plots.  Extrapolation  of  these  curves  to  zero  determines  the  relative  C0(vsl,2) 
populations  which  arise  from  sources  other  than  N2(v=l,2)  populations  which 
arise  from  sources  other  than  N2(v);  i.e.,  from  quenching  of  nitrogen  meta¬ 
stables.  For  the  data  in  Figure  7,  the  multiplicative  factors  are  3.0  and  1.5 
for  C0(v*l,2),  respectively.  The  multiplicative  factors  for  all  examined  CO 
mole  fractions  have  been  determined  to  be  approximately  the  same. 

Using  these  factors  the  contributions  to  the  thermalized  C0(v3l,2)  popula¬ 
tions  from  N2(v)  transfer  and  E-V  transfer  may  be  separated.  Figure  5  shows 
the  C0(v=l,2)  number  densities  corrected, for  N2(v)  transfer.  The  relative 
distributions  in  both  the  thermalized  and  rotationally  excited  components  are 
now  remarkably  similar.  The  vibrational  populations  of  the  rotationally 
excited  component  exceed  those  of  the  thermalized  component  by  typically 
-30  percent.  This  is  approximately  the  ratio  determined  for  all  the  analyzed 
spectra.  Ratios  of  the  vibrational  populations  for  each  component  exhibit  no 
dependence  on  nitrogen  or  carbon  monoxide  mole  fraction. 

As  shown  in  Figure  6  the  fraction  of  the  C0(v=2)  population  attributable 
to  two  quantum  transfer,  reaction  (6b),  is  small.  An  average  ratio 
[C0(va2) J/lC0(v«l) J  of  0.06  ±  0.01  from  N2(v)  transfer  has  been  determined  for 
all  the  observed  spectra.  This  ratio  is  related  to  the  rate  coefficients  for 
reactions  (6a,  b)  in  the  following  manner. 


[C0(v=2)] 

[CO(v-l)] 


Z  k^(Av=2)lN2(v)I 
Z  ky(Av*l)[N2(v)] 


(7) 
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Equation  (7)  may  be  resolved  for  the  ratio  £  ky(Av=2)/£  k.y(Av=l)  if  the  N2(v) 

V  V 

distribution  can  be  determined.  Discharge-flow  measurements  using  Penning 
ionization  spectroscopy  indicate  the  effluent  of  our  microwave  discharges  may 
be  represented  by  a  -6000  K  "modified  Treanor"  distribution  in  N2(v).25  ^ 

Treanor  distribution  contains  enhanced  populations  at  higher  vibrational 
levels,  relative  to  a  Boltzmann  distribution,  that  are  created  by  collisional 
up-pumping. Using  the  modified  Treanor  distribution,  we  have  determined 
£  kv(Av=2)/£  ky(Av=l)  =  0.1  ±  0.04.  Owing  to  the  known  increase  in  ky  as  a 

V  V 

function  of  v,^2,13  jhe  ratio  determined  here  is  likely  representative  of 
higher  v  (probably  7  to  8). 

Figure  8  shows  the  average  population  distributions  for  the  thermalized 
and  rotationally  excited  components,  normalized  separately  and  shown  with  one 
standard  deviation  error  bars.  The  population  distributions  from  only  eleven 
spectra  were  chosen  for  this  average.  The  spectra  containing  large  contribu¬ 
tions  to  thermalized  C0(v=l,2)  have  been  excluded  since  this  emission  envelope 
overlaps  the  v=5-8  levels  of  the  rotationally  excited  component  interfering 
with  reliable  population  determinations  of  these  levels. 

The  distributions  shown  in  Figure  8  are  similar  but  the  distribution  of 
the  rotationally  excited  component  appears  to  be  relatively  flat  above  v=4.  We 
have  examined  the  data  carefully  and  determined  this  trend  to  be  accurate. 

Fits  to  the  data  using  a  fixed  relative  vibrational  distribution  of  the  rota¬ 
tionally  excited  component  following  the  distribution  of  the  thermal  component 
seriously  underfits  the  data  in  the  region  of  v'=5-8  of  the  R-branch  band  heads 
and  in  the  region  of  the  P-branches. 


EXCITATION  MECHANISM 

Several  processes  which  may  account  for  all  or  a  portion  of  the  observed 
emissions  require  consideration.  One  possible  mechanism  for  rotational 
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excitation  of  CO  is  energy  transfer  from  translationally  excited  atoms  emanat¬ 
ing  from  the  discharge.  Translationally  excited  atoms  may  be  produced  in  the 
discharges  by  energy  transfer  from  Ar  metastables.  Translationally  hot  H  atoms 
have  been  reported  to  excite  V,R  states  of  CO. 26  ve  may  discount  such 
processes  in  our  apparatus  due  to  rapid  energy  accommodation  of  any  "hot"  atoms 
in  the  1  torr  discharge  tubes.  The  number  density  of  such  species  reaching  the 
field  of  view  will  be  miniscule.  Collisional  quenching  of  high  CO  vibrational 
levels  into  high  rotational  levels  of  lower  vibrational  levels  (V,R  transfer), 
analogous  to  processes  observed  for  HF,27»28  also  be  dismissed.  To  excite 
the  high  J'  levels  we  observed  would  require  nearly  gas  kinetic  multiquantum 
(Av^)  quenching  of  CO(v).  This  is  unlikely  since  direct  measurements  of  the 
v-dependent  quenching  of  CO(v)  by  C02^^  indicate  the  total  quenching  rate 
coefficients  to  be  less  than  1  percent  gas  kinetic.  Under  quasi-resonance 
conditions,  collisional  quenching  diatomics  in  low  v,  high  J  levels  into  high 
V,  low  J  levels  (R,V  transfer)  may  have  rate  coefficients  of  10~^1  cm^  s“^.30 
The  resonance  conditions  are  given  by 


(0 


Cd 


(8) 


For  CO,  resonance  is  satisfied  at  J'-120  which  is  populated  only  at  the  lowest 
vibrational  level  for  the  statistical  model.  Therefore,  R,V  transfer  should 
not  contribute  significantly  to  the  higher  vibrational  levels  of  the  observed 
thermalized  emissions. 

With  other  possibilities  excluded,  the  mechanism  for  C0(v,J)  excitation 
must  be  energy  transfer  from  one  or  more  of  the  metastable  nitrogen  species 
created  in  the  discharge.  We  have  previously  determined  that  to  account  for 
the  C0(v,J)  excitation  observed  here  the  energy  transfer  reaction  must  be  at 
least  ~3.7  eV  exoergic.  This  constraint,  combined  with  kinetic  considerations 
from  Eq.  (8),  may  be  used  to  determine  the  identity  of  the  responsible  species. 
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Summing  the  populations  of  the  thermalized  and  rotationally  excited  com¬ 
ponents,  a  product  k[M*]  ~  0.5  s“^  is  required  to  account  for  the  observed 
emissions.  The  measured  quenching  kinetics  for  many  of  the  metastable  nitrogen 
species  is  shown  in  Table  2.  This  table  contrasts  the  number  densities  of  the 
metastable  species  required  in  the  interaction  zone  to  account  for  the  observed 
emissions  with  the  number  densities  determined  from  modeling  studies.  The  num¬ 
ber  density  calculations  have  employed  the  room  temperature  rate  coefficients 
in  the  absence  of  data  at  80  K. 

Quenching  of  the  metastable  nitrogen  atoms  are  shown  in  reactions  (9)  and 
(10).  N(2d)  may  deposit  2.38  eV  into  CO 


N(2d)  +  CO  -»  N(^S)  +  C0(v<9) 

AE  =  2.38  eV 

(9) 

N(2p)  +  CO  -»  N(2d)  +  C0(v<4) 

AE  =  1.196  eV 

(10a) 

N(^S)  +  C0(v<15) 

AE  =  3.576  eV 

(10b) 

rovibrational  states  which  is  well  below  the  -3.7  eV  needed 

the  observed  emissions.  Additionally,  the  kinetics  of  this 

to  account  for 

reaction  indicate 

that  contributions  to  the  spectra  form  this  source  would  be  negligibly  small. 
The  energetics  of  N(2p)  quenching  to  form  N(^S)  matches  the  required  ~3.7  eV 
well.  However,  the  recently  determined  rate  coefficient^^  for  quenching  of 
N(2p)  by  CO  is  several  orders  of  magnitude  too  small  for  reactions  (10a,  b)  to 
be  the  C0(v,J)  excitation  mechanism. 

Excitation  of  C0(v<14)  by  energy  transfer  from  N2(X,  high  v)  requires 
multi-quantum  transfer  from  v<14.  The  determined  contributions  to  the 
thermalized  C0(v=l,2)  populations  from  N2(v)  energy  transfer  are  consistent 
wiih  a  preleiieu  channel  for  single-quantum  exchange  and  a  less  favored  channel 
for  two-quantum  exchange.  We  have  found  no  evidence  for  multi-quantum 
exchange. 
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Although  the 


N2(A^Eu~)  quenches  with  CO  to  form  CO(a3n).33 

N2(A3j:y-)  +  CO  -»  N2(X)  +  COCa^fl)  (11) 

C0(a3n)  +  CO  -»  2CO(X,5<v<13)  (12) 

rate  coefficient  for  the  quenching  process  is  known^^  the  absolute  yield  for 
reaction  (11)  has  not  been  measured.  Reaction  of  CO(a)  with  another  CO 
molecule  produces  vibrationally  excited  C0^^»2^  with  a  rate  coefficient  f^  the 
process  near  gas-kinetic  (ki2  ~  1  x  10“^®  cm^  at  room  temperature).^^'^^ 
Quenching  of  N2(A)  may  produce  rovibrationally  excited  CO  directly  from  a 
branch  in  reaction  (11).  Since  N2(A)  lies  6.17  eV  above  N2(X,v=0),  such  a 
process  would  be  sufficiently  energetic.  The  determined  N2(A)  quenching  rate 
coefficient,  however,  is  too  small  for  this  process  or  any  derived  reactions, 
such  as  reaction  (12),  to  be  responsible.  Additionally,  the  activation  barrier 
determined  for  this  process  by  Slanger^®  indicates  the  rate  constant  to  be 
<10“^3  cm^  s”^  at  80  K. 

The  N2(w3Au,  w^Ay)  states  are  sufficiently  energetic  to  produce  the 
observed  C0(v,J)  emissions;  7.36  eV  and  7.35  eV  above  the  ground  state,  respec¬ 
tively.  The  quenching  reactions  of  these  two  species  have  not  been  reported  in 
the  literature.  We  have  determined  N2(w3Au,  v=1-5)  and  N2(w1Au,  v=0-2)  number 
densities  in  COCHISE  directly  from  their  IR  radiance  over  the  2  to  4  ym  region. 
Their  concentrations  in  the  interaction  region  are  3  x  10®  cm^  and  1  x  10®  cm^, 
respectively.  Even  if  these  species  quench  CO  with  rate  coefficients  near  gas 
kinetic,  they  cannot  account  for  the  observed  C0(v,J)  emissions. 

Quenching  of  N2(aln)  is  rapid,  forming  CO(A^n).®^  Since  the 

N2(aln)  +  CO  N2(X,v)  +  CO(Aln)  (13) 

yield  of  C0(A)  formation  from  this  reaction  has  not  been  determined,  a  channel 
forming  C0(v,J)  is  possible.  N2(a)  may  deposit  any  fraction  of  the  available 
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8.4  eV  into  CO  rovibrat ional  states  from  such  a  process.  However,  the  radia¬ 
tive  lifetime  of  N2(a)  is  short,  56+4  which  makes  the  number  density  of 

this  species  in  the  interaction  zone  to  be  too  small  to  account  for  the 
observed  emissions.  The  upper  limit  for  N2(a)  number  density  shown  in  Table  2 
has  been  determined  for  the  noise  level  of  discharged  Ar/N2  spectra  at  the 
wavelenth  N2(a-a')  features  would  occur.^^ 


The  possible  quench  ng  reactions  of  the  N2(a'^J^~)  with  CO  are 

N2(a'li:,j')  CO  N2(X,v)  +  CO(Aln) 

N2(X,v)  +  CO(X,v,J) 

->  N2(X,v)  +  COd^E-,  D^A) 


(lAa) 

(14b) 

(14c) 


N2(a')  may  deposit  up  to  8.5  eV  into  CO  electronic,  vibrational  and  rotational 

states.  The  branching  ratio  forming  C0(A)  has  been  measured  as  20'^o^  percent. 

+40 

This  value,  however,  was  based  on  an  80_2q  WS  lifetime  for  the  N2(a)  state  which 
has  recently  been  revised  to  be  56+4  ys.^^  The  corrected  C0(A)  branching  ratio 
is  30+8  percent.  The  remaining  fraction  must  be  divided  between  the  other 
energetically  accessible  spin-allowed  channels.  These  channels  are  rovibra- 
tionally  excited  C0(X),  reaction  (14b),  and  the  I^E“  and  D^A  states,  reaction 
(14c).  At  room  temperature  the  total  quenching  rate  constant  for  reaction  (14) 
is  1.1  X  10“^C  cm3  g-l.  Table  2  shows  that  the  kinetics  of  this  reaction  are 
sufficient  to  account  for  the  observed  emissions. 


Since  N2(a')  is  the  only  species  present  in  the  interaction  zone  that 
satisfies  both  the  kinetic  and  energetic  constraints,  we  postulate  that 
reaction  (14b)  is  responsible  for  the  observed  C0(v,J)  excitation.  For  this 
hypothesis  to  be  correct,  the  branching  fraction  for  reaction  (14b)  must  be 
relatively  large,  constraining  that  for  reaction  (14c)  to  be  small  and,  owing 
to  the  short  radiative  lifetime  of  C0(A),  the  thermalized  populations  must 
reflect  a  contribution  from  C0(A)  radiative  cascade. 
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DISCUSSION 


Further  interpretation  of  the  data  may  be  assisted  by  surprisal 
theory. ^2, 43  This  approach  postulates  the  existence  of  an  exponential  gap  lav 
for  an  individual  metathetic  reaction: 

Py  -  pO(v)  exp(-Xfv)/exp(\j)  (15) 

where  fy  is  the  fraction  of  the  reaction  exoergicity  appearing  as  vibrational 
energy  in  the  product,  P(v)  is  the  observed  relative  vibrational  population, 
and  P^(v)  is  the  statistical  "prior”  distribution  obtained  when  all  final 
translational,  rotational,  and  vibrational  states  are  equally  probable.  Thus, 
the  ratios  of  the  observed  and  statistical  vibrational  populations  should  be 
exponential  in  the  vibrational  energy  By,  with  the  exponential  fall-off  con¬ 
stant  X  quantifying  the  departure  of  the  observed  distribution  from  a  com¬ 
pletely  statistical  product  distribution.  While  this  theory  does  not 
necessarily  hold  for  all  chemical  interactions,  it  has  proved  useful  in  the 
analysis  of  rovibrational  product  distributions  from  several  photo  chemical  and 
reactive  interactions. 

Both  the  vibrational  extent  of  the  thermalized  CO  emission  (v<14, 

E  -  3.4  eV)  and  the  rovibrational  excitation  of  the  rotationally  excited 
component  (Et  -  3.7  eV)  can  be  attributed  to  an  energy  transfer  process  of  3.5 
to  3.7  eV.  Using  a  prior  distribution  for  3.7  eV  determined  from  the  usual 
relationship  for  a  vibrating  rotator, 

V* 

pO(v)  -  (l-fy)3/2/  z  (l-fy)3/2  (16) 

v.O 


the  vibrational  surprisal  plot  for  the  vibrational  populations  of  the  rotation- 
ally  excited  component  has  been  plotted  in  Figure  9.  The  change  in  slope 
indicate  two  dynamic  mechanisms  are  responsible  for  the  rotationally  excited 
component:  one  accounting  principally  for  v»l-4  and  another  for  the  higher 
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vibrational  levels.  The  relative  populations  of  the  rotationally  excited 
component  may  be  reproduced  by  a  model  composed  of  two  such  mechanisms  as  shown 
in  Figure  10.  The  model  is  comprised  by  roughly  equal  contributions  from  a  low 
V  excitation  process  and  a  "statistical”  process  which  contributes  to  all 
vibrational  levels.  The  distribution  used  for  the  statistical  process  is  the 
prior  calculated  for  3.7  eV.  This  model  predicts  populations  for  v>8  of  the 
rotationally  excited  component.  Absolute  populations  of  these  levels  cannot  be 
determined  from  the  spectra  due  to  the  overlap  of  these  bands  with  the  thermal- 
ized  CO  emission  features.  However,  we  have  performed  fits  to  the  data  using 
v=l-14  of  the  rotationally  excited  component  following  the  modeled  distribution 
shown  in  Figure  10.  .The  resulting  fits  do  not  exhibit  any  marked  differences 
from  fits  excluding  v>8  of  the  rotationally  excited  component  and  the  match  to 
the  absolute  intensity  in  the  P-branch  region  is  generally  improved.  Under 
these  conditions  the  determined  populations  of  the  thermalized  component  for 
v>4  are  decreased  by  an  average  of  ~40  percent.  Thus,  the  populations  for  v>8 
of  the  thermalized  component  shown  in  Figure  8  are  upper  limits. 

Due  to  the  implications  of  the  model  and  the  improvement  in  the  P-branch 
fits,  we  believe  that  v>8  for  the  rotationally  excited  component  are  present  in 
the  data.  However,  higher  resolution  studies  will  be  required  to  confirm  their 
presence  and  determine  the  true  population  distribution. 

Similar  analysis  of  the  rotationally  thermalized  component  must  include 
the  contribution  of  reaction  (14a)  to  these  populations  from  C0(A)  radiative 
cascade.  Quenching  of  N2(a'Ey“,  v=0)  via  reaction  (14a)  yields  C0(A^I1,  v^).^5 
Owing  to  the  short  radiative  lifetime  of  the  C0(A)  state,  ~10  ns,  all  of  the 
C0(A)  created  by  this  process  will  radiatively  cascade  to  form  C0(v<9)  with  a 
vibrational  distribution  reflecting  the  known  branching  ratios.  Since  the 
maximum  energy  defect  for  this  process  is  ~3000  cm“^,  N2(a',v=0)  +  CO  — > 
N2(X,v=0)  +  C0(A,v=0),  we  do  not  expect  substantial  rotational  excitation  in 
the  initially  formed  C0(A)  state  or,  therefore,  in  the  C0(X)  populations 
created  by  radiative  cascade. 
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Thermalization  of  the  rotationally  excited  component  may  also  contribute 
to  the  observed  thermalized  component  populations.  The  fraction  which  may  be 
thermalized  will  depend  on  the  particular  rotational  distributions  and  the 
quenching  rates  at  80  K.  There  are  no  published  quenching  studies  of  high 
rotational  levels  of  CO  at  low  temperatures.  However,  studies  of  C0(v=0,J) 
created  by  formaldehyde  photodecomposition  at  room  temperature  have  determined 
the  rotational  relaxation  rate  constant  for  transfer  of  population  out  of  J=12 
to  be  nearly  gas  kinetic,  3  x  10“!^  cm^  s“^,  and  that  for  J=36  to  be 
5  X  10“^^  cm^  Rotational  levels  in  the  region  J=80  have  rotational 

spacings  of  >300  cm“^,  greatly  exceeding  the  room  temperature  value  of  kT. 

Thus,  quenching  of  these  levels  is  expected  to  be  significantly  slower.  At 
80  K,  kT  exceeds  the  CO  rotational  spacings  for  only  J<15.  Assuming  this  to  be 
the  maximum  number  of  rotational  levels  that  may  be  thermalized  within  the  20 
to  30  collisions  in  the  field  of  view  limits  the  fractional  contributions  from 
the  rotationally  excited  component  to  be  30  percent  for  high  v  and  6  percent 
for  low  V. 

We  may  determine  the  maximum  fraction  of  the  CO  thermal  populations  aris¬ 
ing  from  thermalization  of  the  rotationally  excited  component  by  fitting  the 
observed  thermal  component  population  distribution  to  a  model  composed  of 
thermalized  and  C0(A)  radiative  cascade  contributions.  The  results  of  this 
analysis  are  shown  in  Figure  11.  Reproduction  of  the  observed  distribution 
requires  roughly  30  percent  of  the  initially  formed  C0(v,  high  J)  to  be 
thermalized  and  a  branching  fraction  for  C0(A)  excitation  of  14  percent. 

The  determined  contribution  from  C0(A)  radiative  cascade  is  about  a  factor 
of  two  lower  than  the  recalculated  branching  ratio  of  Piper. ^3  our  result 
could  be  consistent  if  the  branching  ratio  for  reaction  (17a)  decreases 
somewhat  at  lower  temperatures.  Conversely,  the  requirement  of  30  percent 
thermalization  of  the  initially  formed  rotationally  excited  o'"  '  nent  is  larger 
than  the  established  bound.  Inclusion  of  v>8  of  the  rotationally  excited 
component  in  the  *^its  decreases  the  v>4  thermal  populations  by  -40  percent 
•which  reduces  the  fractional  thermalization  proportionately.  Even  with  this 
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reduction,  we  still  require  rotational  relaxation  processes  to  be  faster  than 
our  estimates  or  the  nascent  distributions  to  be  bimodal.  Bimodal  rotational 
distributions  have  been  observed  in  HF  under  near-nascent  conditions. These 
studies  clearly  indicate  that  determination  of  rotational  relaxation  must  take 
the  nature  of  the  excited  complex  into  account.  The  important  factors  are  the 
rate  of  rotation  and  the  relative  velocities  of  the  fragments.  If  the 
velocities  are  low  and  rate  of  rotation  high,  rapid  rotational  thermalization 
will  occur.  The  similarities  between  the  thermal  and  rotationally  excited 
vibrational  distributions  strongly  suggest  a  direct  correlation.  Whether  the 
observed  bimodal  distributions  arise  form  excitation  followed  by 
multi-collisional  thermalization  or  near-nascent  relaxation  of  the  excited 
complex  cannot  be  (distinguished  at  these  experimental  conditions. 

In  summary,  the  results  of  surprisal  analysis  indicates  two  dynamic 
mechanisms  are  responsible  for  the  rotationally  excited  component,  one 
producing  the  low  vibrational  levels  and  another,  more  statistical  mechanism 
contributing  to  all  vibrational  levels.  CO(A)  radiative  cascade  and  a  fixed 
fraction  of  the  rotationally  excited  component  are  sufficient  to  account  for 
the  thermalized  populations. 

The  dynamics  of  the  N2(a' )  +  CO  quenching  process  are  worthy  of  some 
comment.  The  interaction  must  satisfy  conservation  of  angular  momentum  and 
energy. 


(2uEi)l/2b.  =  (2uEf)l/2bj  +  Jh/2n  (17) 

Ei  +  AE  =  Ef  +  Ey  +  Bj2  (18) 

where  u  is  the  reduced  mass,  E^  and  Ef  are  the  initial  and  final  center-of-mass 
collision  energies,  bj  and  bf  are  the  initial  and  final  impact  parameters  AE 
is  the  exoergicity  of  the  process,  and  BJ^  approximates  the  product  rotational 
energy  after  the  collision.  High  rovibrational  excitation  of  CO  from  N2(a' ) 
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quenching  is  certainly  energetically  permissible.  The  observed  CO  rovibra- 
tional  excitation  of  ~3.7  eV  represents  only  45  percent  of  the  total 
exoergicity,  the  remainder  presumably  will  be  manifested  in  the  N2  internal 
modes.  The  manifestation  of  fractional  energies  of  this  magnitude  in  internal 
modes  is  common  to  many  processes  containing  a  large  energy  defect. 

From  Eq.  (17),  large  values  of  J  can  be  achieved  through  highly  attractive 
collisions  with  large  initial  impact  parameters  and  small  final  impact  param¬ 
eters  with  the  concurrent  requirement  that  the  final  center-of-mass  kinetic 
energy  is  small. Large  initial  impact  parameters  have  been  determined  for 
the  quenching  of  Na(32p)+  CO  which  is  2.1  eV  exoergic.  Two  mechanisms  were 
identified,  a  "direct"  process  involving  colinear  Na — CO  approach  with  an 
impact  parameter  of  16aQ  and  a  "complex"  formation  mechanism  from  colinear 
Na--0C  approach  with  an  impact  parameter  13ao.  Both  mechanisms  produce 
rovibrational  excitation  of  CO  with  up  to  0.8  eV  in  rotation.  Extensive 
molecular  orbital  calculations  of  the  Na*  +  CO  interaction  indicated  the  source 
of  the  rotational  excitation  to  be  a  strong  angular  anisotropy  of  the  potential 
energy  surface.  There  are  similarities  between  the  Na*  +  CO  process  and  the 
N2(a')  +  CO  quenching  process  observed  here.  However,  these  similarities 
cannot  be  considered  more  than  qualitative  until  extensive  potential  energy 
surface  calculations  can  be  performed. 


CONCLUSIONS 

Ve  have  observed  extensive  rovibrational  excitation  in  carbon  monoxide 
from  quenching  of  discharged  nitrogen  in  a  low  pressure  cryogenic  apparatus. 
The  CO  fundamental  emission  spectra  exhibit  bimodal  rotational  distributions; 
14  rotationally  thermalized  vibrational  levels  and  eight  rotationally  excited 
vibrational  levels  have  been  identified.  The  spectral  features  due  to  the 
rotationally  excited  component  are  best  fit  by  a  statistical  distribution, 

E  =  3.7  eV.  This  distribution  provides  sufficient  population  in  the  region  of 
the  Fortrat  reversal,  J  ~  90,  to  account  for  the  observed  R-branch  band  head 
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formation.  An  energy  transfer  process  of  ~3.5  eV  is  sufficient  to  account  for 
the  vibrational  extent  of  the  thermal  component.  These  two  energetic  limits 
are  not  different  in  our  analysis. 

Analysis  of  the  data  has  identified  a  contribution  to  C0(v=l,2)  of  the 
thermalized  component  from 

N2(v)  +  CO  ^  N2(v-1)  +  CO(v=l) 

N2(v)  +  CO  ^  N2(v-2)  +  C0(v=2) 

Based  on  kinetic  and  energetic  arguments  we  have  determined  the  energy  transfer 
process 

N2(a'^Eu~)  +  CO  N2(X,v)  +  CO(Aln) 

N2(X,v)  +  C0(v,J) 

to  be  responsible  for  the  observed  rovibrational  excitation.  The  branching 
ratio  for  CO(A^n)  formation  is  30±8  percent.  In  order  to  account  for  emissions 
observed  here,  CO(v,J)  excitation  must  account  for  the  remainder. 

The  energy  transfer  mechanism  producing  the  rotationally  excited  component 
exhibits  dynamic  similarities  to  that  determined  for  Na(^P)  +  CO;  two 
mechanisms  have  been  identified,  one  responsible  for  low  vibrational  levels  and 
another  that  contributes  to  all  accessible  vibrational  levels.  Rigorous 
potential  energy  surface  calculations  will  be  required,  however,  to  determine 
the  true  dynamics  of  the  N2(a')  +  CO  interaction. 
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Table  1.  Energetics  of  CO  Band  Heat  Formation  E-j-  =  3.70  eV 


Vibrational 

Level 

E^(eV)^ 

J 

vertex 

Ej  (eV) 

vertex 

J 

max 

1 

0.266 

3.434 

92 

1.966 

123 

2 

0.528 

3.172 

91 

1.907 

118 

3 

0.788 

2.912 

90 

1.848 

114 

4 

1.044 

2.656 

90 

1.831 

109 

5 

1.296 

2.404 

89 

1.774 

104 

6 

1.546 

2.154 

88 

1.719 

98 

7 

1.792 

1.908 

87 

1.665 

93 

8 

2.035 

1.665 

87 

1.648 

87 

9 

2.274 

1.426 

86 

1.600 

81 

10 

2.511 

1.189 

85 

1.544 

74 

11 

.744 

0.956 

84 

1.444 

66 

12 

2.974 

0.726 

84 

1.478 

58 

13 

3.201 

0.499 

83 

1.430 

48 

14 

3.425 

0.275 

82 

1.382 

35 

a-  Eco(v=0)  =  0 

b.  Calculated  from  the  value  for  Ej-Ev 


Table  2.  Measured  CO  Quenching  Kinetics  of  N2*,  N* 


Number  Density  in 

COCHISE,  cm-3 

Species 

Kq  (p  K) 
cm"3  s"^ 

Required  for 
C0(v,J)  Excitation 

Estimated  for 
Interaction  Zone 

N(2p) 

<1.5  (-14)31 

>3 

(13) 

1  to  3  (9) 

N(2d) 

N2(a3V). 

1.7  (-12)32 

>3 

(11) 

3  to  10  (9) 

v=0 

1.5  (-12)^33 

>3 

(11) 

1  to  3  (9) 

v=4 

1.9  (-ll)b33 

3 

(10) 

1  to  3  (9) 

N2(alng) 

2.8  (-10)^34 

2 

(9) 

2  (6) 

N2(a'lEu") 

1.1  (-10)a35 

4 

(9) 

3  (9) 

a.  Includes  excitation  of  CO(Aln) 

b.  Includes  excitation  of  C0(a3n) 


FIGURE  CAPTIONS 
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Figure 

Figure 

Figure 

F’gure 

Figure 

Figure 

Figure 

Figure 

Figure 


1.  Diagram  of  the  COCHISE  reaction  chamber.  The  physical  dimensions 
of  the  cell  are  0.6m  in  length  and  0.4m  in  diameter. 

2.  Data  comparison.  The  dark  line  shows  a  typical  spectrum  obtained 
under  conditions  of  12  percent  discharged  N2/Ar  reacting  with  a 
couterflow  of  33  percent  CO/Ar.  The  light  line  shows  data  taken  at 
3.17  percent  N2/Ar  with  identical  counterflow  conditions. 

3.  Comparison  of  rotational  distributions.  Shown  are  the  relative 
rotational  distributions  for  a  10^  K  Boltzmann  and  several 
statistical  distributions. 

4.  Data  (light  line)  and  best  fit  (dark  line)  (a)  to  thermalized 

(80  K)  CO  fundamental  emission  (v'=l-14)  and  CO  rotatinally  excited 
bands  using  a  statistical  distsribution  with  Ej  =  3.7  eV.  The  data 
is  the  same  used  for  the  light  line  in  Figure  2.  The  spectral 
resolution  is  0.013  pm.  The  thermalized  CO  and  rotationally 
excited  basis  sets  which  comprise  the  best  fit  in  (a)  are  shown  in 
tb)  and  (c),  respectively. 

5.  Enlarged  view  of  the  5.7  to  6.5  us  region  of  Figure  4.  The  data  is 

shown  by  light  line,  the  fit  by  a  dark  line. 

6.  Determined  population  distribution  versus  By,  for  the  fit  shown  in 

Figure  4.  The  thermalized  component  is  shown  as  (  ),  with 

contribution  to  thermalized  C0(v=l,2)  as  (  ),  and  the  rotationally 

excited  component  is  represented  as  (  ). 

7.  Plot  of  the  ratio  of  the  [C0(v=l,2))  to  [C0(v=3)]  populations  as  a 

function  of  nitrogen  mole  fraction.  The  CO  mole  fraction  for  these 

data  was  0.33. 

8.  Averaged  relative  population  distributions  for  11  spectra,  both 

thermalized  (  )  CO  and  the  rotationally  excited  component  (  )  are 

shown  but  have  been  normalized  independently.  The  error  bars 
represent  one  standard  deviation. 

9.  Vibrational  surprisal  plot  for  the  rotationally  excited  vibrational 
populations.  The  populations  are  those  from  Figure  8. 

10.  Model  fit  to  the  rotationally  excited  vibrational  populations.  The 
model,  shown  by  solid  line,  has  been  produced  by  equal 

constributions  from  the  direct  (v>4),  ( - ),  mechanism  and  the 

statistical  mechanism  (all  v),  (•••)• 
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Figure  11.  Model  fit  to  the  rotationally  "thermal"  vibrational  populations. 

The  model  ( _ )  has  been  produced  by  a  contribution  of  14  percent 

branching  fraction  from  reaction  (14a),  (- _ -),  and  30  percent  of 

the  rotationally  excited  model,  shown  separately  as  direct  ( — )  and 
statistical  contributors  (•-•)•  The  relative  populations  have  been 
normalized  to  the  sum  of  the  rotationally  excited  vibrational  popu¬ 
lations.  The  rotationally  excited  component  is  also  shown  (*••)• 
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1  .  INTRODUCTION 


The  usual  mechanism  for  the  recombination  of  ozone  from  oxygen  atoms  and 
molecules  involves  a  two-sttp  process  in  which  a)  an  unbound  complex  of  ozone, 
O3* ,  13  formed  and  b)  O3*  is  stabilized  coll isionally  to  form  bound  ozone. 

The  reaction  scheme  is  written  as  the  following; 


0 


3 


k 

0.  *  M  I  0^{\/,J) 
u  3 

where  M  denotes  a  third  body  collision  partner  and  03(v,J)  are  the 
vlbrationally-rotationaily  excited  ozone  molecules.  It  is  generally  assumed 
that  the  O3*  formation  is  an  equilibrium  or  steady  state  process  so  that 

* 

^  “  [o^JCo] 

This  is  the  energy  transfer  recombination  mechanism. 


Then  the  so-called  one  way  flux  of  ozone  molecule  production  is  given  by 

dlo  1 

To  compute  the  observed  reaction  rate  of  ozone  formation  the  redissociation  of 
the  highly  excited  nascent  ozone  must  be  included  according  to  the  following. 


0^(v,J)  +  M  0  +  +  M  . 

This  reaction  reduces  the  rate  of  formation  from  the  one-way  flux  value. 

An  alternative  mechanism  for  ozone  formation  involves  the  formation  of  an 
oxygen  atom-third  body  bound  complex 

K 

M  +  0  =’  MO 
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which  collides  with  oxygen  molecules  to  form  stabilized  (but  probably  highly 
internally  excited)  ozone  molecules 

•^3 

MO  +  >  0^(v,J)  +  M  . 

The  ozone  formation  one-way  flux  is  given  by 
dlO  1 

-  kjlMOHOjJ  .  . 

This  is  the  so-called  chaperon  mechanism.  As  before,  it  is  necessary  to 
include  the  collislonal  redissociation  of  highly  excited  ozone  molecules  into 
the  computation  of  the  observable  reaction  rate. 

In  the  upper  atmosphere,  at  altitudes  of  Interest  for  ozone  formation, 
the  primary  constituents  are  O2,  N2  and  0  and  it  is  usually  assumed  that  the 
radical  energy  transfer  mechanism  is  dominant.  However  under  circumstances  in 
which  the  third  body  collision  partner  is  in  high  concentration  (and  is  not 
oxygen)  the  chaperon  mechanism  may  be  important.  The  COCHISE  experiments  are 
run  under  conditions  in  which  there  is  a  large  excess  of  argon.  In  order  to 
assess  the  potential  contributions  of  the  chaperon  mechanism,  trajectory  cal¬ 
culations  were  carried  out  for  the  reaction  of  ArO  with  O2  to  form  ozone. 

We  discuss  two  aspects  of  these  calculations:  1)  calculation  of  the  ArO 
equilibrium  constant  and  2)  classical  trajectory  calculations  of  ArO-02  colli¬ 
sions  and  redissociation  of  newly  formed  ozone  molecules  by  collision  with 
argon  atoms. 
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2.  ArO  EQUILIBRIUM  CONSTANT 


We  assume  that  the  ArO  interaction  potential  is  a  Lennard-Jones 
SIX- twelve  potential  given  by 

V(r)  =  4el('7/r)^^  -  (a/r)^] 

wnere  e  IS  the  well  depth  and  a  a  molecular  parameter.  This  form  for  the  ArO 
Interaction  has  been  used  in  our  earlier  calculations  on  energy  transfer  in 
Ar-Og  collisions  and  radical  mechanism  ozone  recombination.  Values  for  e  and 
’  were  determined  by  Aquilanti  et  ai.^  using  a  molecular  beam  scattering 
technique.  The  adopted  values  are: 


and 


£  =  0.0047  e.v. 

0  =  2.78A. 


The  expression  for  the  equilibrium  constant  of  the  ArO  bound  complex  is 
given  by2 


K  = 


4n 

(g;  a 


j  dr  r  exp[ 


-V(r) 

kT 


]  7(2’  } 


where  r(3/2)  and  y ( 3/2 , -V( r ) /kT )  are  complete  and  incomplete  gamma  functions 
respectively.  This  expression  was  derived  by  Hill. 2  Stogryn  and 
Hirschf elder-^  showed  this  result  may  be  simplified  when  V(ri  is  given  by  a 
Lennard-Jones  6-12  potential.  Using  the  notation  of  Kim,'^  the  equilibrium 
constant  may  be  rewritten  as 


128  2,  3, 

K  =  — —  a  (ttz  ) 
45 


1/2 


F(z) 


(  1  ) 
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where 


F{z)  -  I 


nl 


n“0  (l)  (!) 


n  „  / ,  ,  7  9  . 

^  “  2^2  ’’'4’4'^ 


■4^ 

n  n 

{a)n  “  a(a+i)  ...  (a+n-1) 
and  z  “  e/kT 

An  asymptotic  form  in  the  limit  z  -  e/kT  ->•  ®  is  given  by4 


(-)  (-) 

N-1  '-4-'  '■4^ 

,  1  Stt  -2  z  r  „  n  n  -n  -N>,  i 

F(z)  -  - z  e  [  E  - j -  z  +  0[z  Jj 

32^2  n*0 


In  this  limit  the  equilibrium  constant  may  be  written  as 

2  2,^  ,  ,  J/2  e/k.T 

K  ~  ~  Tra  (2irkT/e)  e 


(2) 


This  form  is  appropriate  for  low  temperatures.  The  lowest  temperature  we 
shall  consider  is  80  K  for  which 

kT  »  0.16  kcal/mol  »•  0.0069  e.v. 


so  that  we  shall  always  work  in  the  z  <  1  case. 


The  above  equilibrium  constant  includes  contributions  from  ArO  molecules 
below  the  dissociation  limit.  Classically,  molecules  above  the  dissociation 
energy  (i.e.,  well  depth)  but  trapped  within  a  centrifugal  barrier  may  also  be 
included.  Kim  and  Ross^  computed  equllibriiun  constants  including  contribu¬ 
tions  from  bound,  metastable  and  orbiting  molecules.  Their  form  for  the  equi¬ 
librium  constant  is  given  by: 


K 


8  3 

Ts 


(3) 
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where 


=rel  -  ^  ’ntfe)” 

n-0 

and  {q^}  are  numerical  coefficients.  Values  for  the  first  ten  coefficients 
are  listed  in  the  Kim  and  Ross  paper. ^ 

In  the  calculations  presented  here  we  consider  only  the  contributions 
from  bound  ArO  complexes.  Our  temperature  range  of  Interest  Is  between  80  and 
400  K  so  that  Eq.  (1)  will  be  used  to  compute  the  ArO  equilibrium  constant. 

Calculated  values  for  the  ArO  formation  equilibrium  constant  are  given  as 
a  function  of  temperature  In  Table  1.  The  values  of  the  equilibrium  constant 
for  formation  of  bound  ArO  molecules  Is  denoted  as  Kb*  comparison  we  also 

Include  values  of  the  equilibrium  constant  Including  metastables  and  orbiting 
pairs  (Eq.  (3)  and  (4)),  Km.  Two  points  are  noteworthy:  (1)  as  expected  the 
value  of  the  equilibrium  constants  Kb  and  Km  Increase  monotonlcally  with  de¬ 
creasing  temperature  and  (2)  Km  Is  approximately  a  factor  of  two  larger  than 
K|^  at  all  temperatures  considered.  Note  also  that  the  temperature  variation 
of  both  Kb  and  Km  Is  quite  slow,  the  values  varying  only  an  order  of  magnitude 
between  400  and  80  K.  It  Is  clear  that  the  observed  much  steeper  temperature 
dependence  In  the  ozone  recombination  rate  cannot  be  explained  by  the  tempera¬ 
ture  dependence  of  the  equilibrium  constant  alone.  In  order  to  find  agreement 
with  experiment  there  must  be  strong  temperature  dependencies  for  the  associa¬ 
tion  reaction 


MO  +  O2  ->■  O3  +  M 

and  redlssoclatlon  reaction 

Oj  +  M  02  +  O  +  M  . 

Results  will  be  discussed  later. 
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Table  1 .  Computed  Values  of  ArO  Equilibrium  Constant 


Temperature 

{ cm^/molecule ) 
Equation  (1) 

Km 

( cm^/molecula ) 
Equations 
(3)  and  (4) 

80 

7.33x10-23 

1 .40x10-22 

100 

5.05x10-23 

1  .03x10-22 

120 

3.74x10-23 

7.94x10-23 

150 

2.61x10-23 

5.78x10-23 

200 

1 .66x10-23 

3.8  xlO-23 

300 

8.79x10-24 

2.12x10-23 

400 

5.64x10-24 

1.39x10-23 

Estimates  of  bound  ArO  concentrations  are  given  in  Table  2.  Conditions 
were  chosen  to  approximately  simulate  those  of  the  COCHISE  experiment  with 

'^Ar  *  10^^  atoms/cm^ 
and  no  “  10^^  atoms/cm^ 

To  compute  production  rates  of  ozone  molecules  according  to  this  mechan¬ 
ism,  trajectory  calculations  were  beg\m.  They  were  used  to  compute  the  asso¬ 
ciation  reaction  and  redissociation  reactions.  We  now  mention  some  details  of 
the  computation. 
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Table  2 


Computed  Values  of  ArO  Concentrations;  n^^^  =  10^^  atoms/cm^  and 
no  =  10^^  atoms/cm^,  n^^^-Q  ■  K^’^Ar'^O 


T 

^ArO 

(molecules/cm®) 

80 

7.3x10® 

100 

5.0x10® 

120 

3.7x10® 

150 

2.6x10® 

200 

1.7x10® 

300 

8.8x10® 

5.6x10® 


400 


3.  TRAJECTORY  CALCULATIONS 


All  calculations  were  performed  using  the  Murrell-Varandas  ozone  ground 
state  potential  surface  with  pairwise  additive  Lennard-Jones  potentials 
between  oxygen  and  argon  atoms.  This  potential  surface  has  been  discussed  in 
detail  elsewhere.^ 

Trajectories  were  computed  by  integrating  the  classical  equations  of 
motion.  The  trajectory  analysis  involves  a  three  step  procedure:  (1)  selec¬ 
tion  of  initial  conditions,  (2)  integration  of  Hamilton's  equations  of  motion, 
and  (3)  determination  of  asymptotic  conditions,  i.e.,  products.  Each  of  these 
steps  is  discussed  in  detail  in  Refs.  7  and  8.  We  mention  here  only  details 
specific  to  the  present  calculations. 

Initial  conditions  for  the  ArO  molecule  were  chosen  as  if  the  molecule 
were  bound  classically,  i.e.,  in  terms  of  its  total  energy  below  the  Ar-0  well 
depth.  The  vibrational  energy  of  ArO  was  chosen  from  an  appropriate  classical 
Boltzmann  distribution  with  a  cut-off  at  the  well  depth.  For  simplicity  the 
ArO  rotational  energy  was  set  equal  to  zero.  This  assumption  was  checked  by 
performing  trajectories  with  rotating  molecules  and  comparing  with  the  results 
for  non-rotating  cases.  It  was  found  that  initial  rotational  ArO  energy  made 
negligible  difference  in  the  outcome  of  trajectories.  This  was  so  because  the 
allowed  rotational  energy  is  very  small  because  of  the  shallowness  of  the  ArO 
potential  well.  Initial  conditions  for  the  initial  O2  molecules  were  chosen 
from  classical  rotational  and  vibrational  energy  Boltzmann  distributions. 
Relative  translational  energies  were  chosen  from  the  appropriate  classical 
flux  distribution. 

Hamilton's  equations  were  integrated  using  a  fixed  step  Runge-Kutta-Gill 
integrator  for  the  first  five  time  steps  and  a  fixed  stepsize  Adams-Moulton 
integrator  for  the  rest  of  the  trajectory.  Trajectories  were  checked  for 
accuracy  by  observing  conservation  of  energy  and  angular  momentum. 
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The  Initial  ArO-02  separation  was  chosen  so  that  the  interaction  between 
the  molecules  was  negligible.  Due  to  the  relatively  low  temperatures  consid¬ 
ered  this  necessitated  a  fairly  large  initial  separation.  This  large  initial 
separation  and  also  large  exit  channel  distance  cutoff  caused  the  trajectories 
to  be  fairly  time  consuming,  with  many  trajectories  requiring  some  tens  of 
thousands  of  integration  steps.  This  was  particularly  true  at  the  lowest  tem¬ 
perature  considered. 

Simulations  of  the  association  reaction  were  carried  out  at  temperatures 
of  80,  100,  120,  150,  200,  300  and  400  K.  Batches  of  between  6000  and  12000 
trajectories  were  computed  at  each  temperature. 

At  the  end  of  each  tra^iectory  it  was  determined  if  a  bound  ozone  molecule 
was  formed.  If  so,  the  ozone  coordinates  were  recorded  and  stored.  These 
coordinates  were  used  as  initial  conditions  in  the  redissociation  reaction 
computations . 

The  multiple  collision  method  discussed  in  Ref.  9  was  used  to  compute  the 
redissociation  probability  of  newly  formed  ozone  molecules  in  an  argon  thermal 
bath.  Available  computer  time  only  allowed  these  calculations  to  be  carried 
out  at  temperatures  of  400,  300  and  200  K.  At  these  temperatures  the  number 
of  computed  initial  ozone  bound  state  initial  conditions  was  fairly  small. 
Therefore  each  set  of  initial  ozone  coordinates  was  used  several  times  in 
order  to  generate  a  reasonably  large  ensemble  of  nascent  ozone  molecules.  The 
initial  ensemble  sizes  varied  between  600  and  750  initial  ozone  molecules.  At 
300  and  400  K  the  ensembles  were  propagated  through  a  sequence  of  sixty  colli¬ 
sions.  In  the  200  K  case  computer  restrictions  forced  stopping  at  20  colli¬ 
sions.  This  value  seemed  sufficient  to  determine  the  redissociation  probabil¬ 
ity  at  200  K. 
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4.  RESULTS  AND  DISCUSSION 


Computed  values  of  the  cross  section  for  the  association  reaction  of  ArO 
and  O2  to  form  ozone  are  shown  in  Figure  1.  The  cross  section  is  seen  to 
rapidly  increase  with  temperature.  At  the  highest  temperatures,  300  and 
400  K,  the  fraction  of  trajectories  leading  to  stable  ozone  formation  is  ap¬ 
proximately  0.002.  The  fraction,  P,  is  readily  computed  from  the  values  in 
Figure  1  by  the  following 

2 

P  =  c^react/^bmax 

where  bmax  is  the  maximum  impact  parameter  used  in  the  trajectory  calculations 
(equal  to  4A ) .  The  fraction  of  reactive  trajectories  increases  with  decreas¬ 
ing  temperature  to  approximately  0.05  at  80  K.  Even  at  the  lowest  temperature 
considered  the  probability  of  forming  stable  ozone  is  quite  small.  The  nas¬ 
cent  ozone  molecules  were  found  to  be  very  vibrationally-rotationally  excited 
with  total  molecular  energies  very  near  the  dissociation  limit.  Average  ozone 
molecular  energies  for  the  nascent  ozone  molecules  at  200,  300  and  400  K  are 
listed  in  Table  3.  At  these  temperatures  the  mean  newly  formed  ozone  mole¬ 
cule  lies  approximately  one-half  kcal/mol  below  the  dissociation  limit.  It  is 
clearly  anticipated  from  earlier  work  on  Ar-03  dissociation  dynamics  that 
there  will  be  a  large  degree  of  redissociation  as  the  nascent  ozone  molecules 
undergo  collisions  with  the  thermal  argon  bath. 

Table  3.  Average  Energy  Below  Dissociation  Energy  for  Nascent  Ozone 


Energy  Below 

Temperature 

Dissociation 

(kcal/mol) 

200 

0.58 

300 

0.63 

400 

0.49 
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Figure 


The  fraction  of  nascent  ozone  molecules  undergoing  dissociation  as  a 
function  of  bath  temperature  is  shown  In  Figure  2.  At  400  K  almost  80  percent 
of  the  nascent  ozone  molecules  dissociate.  The  dissociation  probability 
decreases  at  the  lover  temperatures  to  0.70  and  0.57  at  300  and  200  K 
respectively.  Since  this  quantity  is  necessary  for  the  computation  of  the 
overall  rate  we  will  use  a  linear  extrapolation  of  the  available  data  to 
predict  dissociation  fractions  at  lower  temperatures.  This  leads  to  the 
following  values. 


^dlss 

Bath  Temperature  (K) 

0.5 

150 

0.35 

120 

0.25 

100 

0.10 

80 

These  values  are  reasonable  from  our  systematic  study  of  ozone  dissocia¬ 
tion  In  argon  baths. In  that  work  ozone  molecules  were  initialized 
25  kcal/mol  above  the  molecular  well  minimum  (or  about  1.5  kcal  below  the  dis¬ 
sociation  limit) .  The  average  translational  energy  In  a  Ar-03  collision  xs 
2  kT.  In  the  case  under  consideration  here,  at  400  K,  the  average  relative 
translational  energy  is  1 .6  kcal/mol  so  that  the  ratio  of  available  energy  for 
dissociation  to  energy  required  for  dissociation  (Table  2)  Is  about  3.  The 
values  of  this  ratio  for  the  300  and  200  K  cases  are  about  2  and  1.6.  We 
shall  assume  that  the  dissociation  probability  Is  only  dependent  on  the  value 
of  the  above  ratio  so  that  we  may  use  our  previous  results  to  Infer  dissocia¬ 
tion  probabilities  at  lower  temperatures. 

To  check  this  assertion  we  note  that  In  the  previous  work  for  a  1000  K 
bath  1.5  kcal/mole  below  the  dissociation  limit,  the  ratio  of  available  to 
needed  energy  for  dissociation  is  1.33  and  the  computed  dissociation  probabil¬ 
ity  was  found  to  be  0.61.  This  compares  favorably  with  the  value  of  0.57  for 
a  ratio  of  1.6  at  200  K.  Clearly  ^is  procedure  requires  more  elaborate  and 


744 


TEMPERATURE  (K) 


745 


Function  of  Bath 


extensive  examination.  However  we  shall  use  It  here  without  further  ado.  We 
also  assume  the  nascent  ensemble  average  energy  to  be  constant  with  tempera¬ 
ture  and  equal  to  0.6  kcal/mol  below  the  well  depth.  The  following  values  for 
the  dissociation  probabilities  are  found. 


^diss 

Bath  Temperature  (K) 

0.45 

150 

0.25 

120 

0.15 

100 

0.02 

_ 

80 

This  procedure  yields  values  sufficiently  close  to  the  simple  linear 
extrapolation  values  to  be  in  agreement.  Note  both  procedures  are  used 
because  of  necessity  and  not  desirability.  The  latter  values  will  be 
adopted. 


We  now  compute  effective  ozone  recombination  rates  via  the  chaperon  mech¬ 
anisms.  According  to  earlier  discussion  the  net  production  rate  of  ozone, 
d[03j/dt,  is  given  by 


d[0  J 
_ 3_ 

dt 


[0][0_][M](1 
3  b  2 


diss 


) 


where  ^3  is  the  rate  constant  for  the  association  reaction  given  by 
^3  »  <va> 


where  a  is  the  cross  section  for  association,  v  the  relative  velocity  and  the 
brackets  denote  a  thermal  average .  The  value  of  k3  is  adequately  approximated 
as 


kg  »  V(T 
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where  v  is  the  average  relative  thermal  velocity.  Therefore  the  overall 
chaperon  mechanism  ozone  formation  rate  constant,  k,  is  given 

k  =  vaKb(l-P)  . 

Computed  values  of  k  are  shown  in  Figure  3  as  a  function  of  temperature. 
Recent  experimental  results  are  shown  in  Figure  4. 

The  trajectory  calculation  values  for  the  ozone  recombina tion  rate  con¬ 
stant  are  clearly  much  smaller  than  measured  experimentally.  This  discrepancy 
may  have  several  sources:  1)  the  chaperon  mechanism  may  be  ineffective  in 
ozone  recombina tioh  and  the  recombination  is  dominated  by  the  energy  transfer 
mechanism,  and  2)  uncertainties  in  the  Ar-Og  potential  surface. 

The  primary  goal  of  this  effort  is  to  assess  the  importance  of  the  chap¬ 
eron  mechanism.  Therefore  the  first  potential  source  is  a  result  we  seek  and 
may  not  be  invoked  to  explain  any  discrepancy  between  theory  and  experiment. 
The  potential  surface  used  in  these  calculations  is  a  state-of-the-art  poten¬ 
tial  surface  empirically  determined  ozone  surface  and  the  best  one  currently 
available.  The  ArO  potentials  are  from  molecular  beam  scattering  data  and  are 
expected  to  be  reasonably  accurate.  In  any  case  the  anticipated  uncertainties 
are  not  expected  to  produce  orders  of  magnitude  variations  in  the  dynamical 
results.  However  this  conclusion  must  be  viewed  as  only  an  assertion. 

The  low  values  of  the  ozone  recombination  rate  constant  are  due  to  the 
magnitude  of  the  equilibrium  constant  for  ArO  complexes  and  the  small  exchange 
reaction  cross  section  for  the  reaction  of  ArO  -t-  03.  The  small  exchange  cross 
section  may  be  seen  as  due  to  the  difficulty  with  which  the  relatively  stiff 
oxygen  molecule  bond  is  stretched  to  reach  the  ozone  equilibrium  configuration 
and  the  difficulty  of  transferring  the  released  bond  energy  to  relative  argon 
translational  energy.  The  small  value  of  the  ArO  equilibrium  constant  is 
a  result  of  the  small  ArO  well  depth.  A  factor  of  two  increase  in  the  well 
depth  only  leads  to  increases  of  4-5  over  the  current  value. 
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Figure  3.  Computed  Ozone  Recombination  Rate  Constant  as  a  Function 
of  Temperature. 


Within  the  framework  of  these  calculations  we  may  conclude  that  the 
chaoeron  mechanism  does  not  play  a  major  role  in  ozone  recombination  in  the 
COCHISE  environment. 

As  a  final  point  we  note  that  the  temperature  dependence  of  the  computed 
rate  constant  is  in  fairly  good  agreement  with  the  experimental  values  at  low 
temperatures.  This  agreement  would  seem  to  make  the  chaperon  mechanism  a  dis~ 
turblngly  interesting  candidate  mechanism.  However  as  noted  above  its  magni¬ 
tude  is  substantially  below  experimental  values.  A  more  complete  comparison 
of  the  chaperon  mechanism  with  the  energy  transfer  mechanism  requires  computa¬ 
tion  of  the  rate  constant  via  the  energy  transfer  mechanism.  Such  calcula¬ 
tions  have  not  yet  been  performed. 
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1 .  INTRODUCTION 


According  to  the  radical  molecule  mechanisms  the  recombination  of  ozone 
occurs  by  the  following  steps: 

>^1  * 

0  +  O2  J  O3  (R1) 

*  1^3 

O3  +  M  +  03(v,J)  +  M  (R2) 

Reaction  ( 1 )  Involves  the  collision  of  an  oxygen  atom  with  an  oxygen  molecule 
to  form  an  unbound  ozone  complex,  O3*.  Since  O3*  Is  formed  above  the  disso¬ 
ciation  limit  of  ozone  It  has  a  finite  lifetime  and  will  dissociate  unless 
stabilized  by  a  collision  with  another  molecule  or  atom  (denoted  as  M  In 
Reaction  (2)).  The  newly  recombined  ozone  molecules,  03(v,J),  will  In  general 
be  very  vlbratlonally-rotatlonally  excited  and  undergo  collisions  with  bath 
molecules  that  will  lead  to  either  relaxation  or  redlssoclatlon: 

k4 

03(v,J)  +  M  +  03(v',J')  +  M  (R3) 

or 

kd 

O  {v,J)  +M+  0+0+M  .  (R4) 

The  net  rate  of  production  (In  the  bath)  of  stable  ozone  molecules  Is 
given  by  the  rate  of  production  of  nascent  ozone  corrected  by  the  fraction  of 
nascent  molecules  that  subsequently  dissociate,  P(j,  l.e., 

d[0  ] 

-JJ-  -  kj[03*][Ml(t  -  P^l  ■  Cl 

The  computation  of  the  dissociated  fraction  requires  a  multiple  collision 
treatment. 

In  this  report  we  discuss  the  results  of  a  classical  trajectory  study  of 
the  formation  of  long-lived  complexes  of  ozone,  their  stabilization  by 
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collisions  with  argon  atoms  and  subsequent  relaxation  of  newly  formed  ozone 
molecules  in  a  200  K  argon  atom  thermal  bath.  The  calculations  were  designed 
to  gain  an  understanding  of  the  molecular  quantities  that  control  each  of  the 
processes  involved  in  the  overall  recombination  reaction  in  order  to  aid  in 
the  interpretation  of  available  experimental  data  from  the  COCHISE  experiments 
at  AFGL. 

The  organization  of  this  report  is  as  follows.  We  begin  with  a  discus¬ 
sion  of  the  classical  trajectory  formulation  for  calculation  of  the  ozone 
recombination  dynamics.  A  brief  discussion  of  the  computational  details  and 
calculations  performed  follows.  Results  and  discussion  are  then  presented  and 
finally  the  implications  of  this  work  to  the  COCHISE  experiments  are  given. 
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2.  CLASSICAL  IRAJECTORV  METHODOLOGY 


A  Steady  state  approximation  for  the  O3*  concentration  yields 


[OJ  - 


kitO] [O^ 
k2  +  k3[M] 


The  value  of  ^2“^  is  the  average  lifetime  of  an  O3*  complex  and  the  term  k3[M] 
Is  the  rate  of  colllslonal  removal  of  O3*.  We  will  assume  that  the  total  gas 
pressure  Is  sufficiently  low  and  the  complex  lifetime  sufficiently  short  so 
that 


k2  >>  k3[M]  .  (4) 

This  Is  surely  valid  at  pressures  up  to  several  atmospheres.  The  ozone 
complex  concentration  becomes 

.  k  [o]to  ] 

[O3]  -  ■  --  -  =  KCDlolto^] 

where  K(T}  Is  the  equilibrium  constant  for  the  formation  of  an  unbound  ozone 
complex.  The  rate  of  production  of  stable  ozone  molecules  Is  now  written 
from  Eq.  ( 1)  as 

-  k ,K(T)[0][0,][MJ(1-P^)  .  (5) 

~dF“  ^  2  d 

The  dynamical  simulation  of  the  ozone  recombination  process  requires 
three  separate  computations.  First,  a  complex  equilibrium  constant  K(T)  must 
be  computed.  Second,  the  single  collision  rate  constant  for  complex  stabili¬ 
zation,  k39,  must  be  determined,  and  finally  the  multicolllslonal  redlssocla- 
tlon  probability  of  newly  formed  ozone  molecules  is  computed  (R4}. 
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We  follow  the  work  of  Bunker^  and  Blint^  for  the  calculation  of  the  ozone 
complex  formation  equilibrium  constant.  The  complex  equilibrium  constant  may 
be  related  to  the  cross  section  for  the  formation  of  an  O3*  complex  by  colli¬ 
sion  of  O2  with  0,  0(EY,Ej^,Et)  at  fixed  O2  vibrational  energy,  Ey,  rotational 
energy,  Er  and  relative  translational  energy,  E^,  by  the  following  relation 


K(T) 


r('l'|1/2  ('iL_')3/2  — L_ 
UirJ  '-kT'*  QrQv 


y  V  Ey)AT 


(6) 


V  -E  At 

[/  (t(E^,Ej^,E^)  e  vdv]] 


where  Qr  is  the  diatom  partition  function,  the  vibrational  partition 

function,  V  the  atom-diatom  relative  velocity  and  y  the  reduced  mass.  In  a 
purely  classical  treatment  the  sums  are  replaced  by  integrals. 


A  relationship  is  needed  between  the  complex  formation  cross  section  and 
classically  computed  quantities.  This  relationship  is 


09 


where  T(Ev>^R>^t^  lifetime  of  the  complex  and  P(E^,Ej^,E^)  is  the 

complex  formation  probability. 

To  compute  the  values  of  T(Ev>^R>£t)  P(Ev,Er,E^),  the  phase  space 

volume  of  an  ozone  collision  complex  must  be  defined.  A  convenient  although 
not  unique  procedure  is  that  a  complex  is  formed  whenever  the  distance  between 
either  of  the  diatom  oxygen  atoms  and  the  incoming  oxygen  atom  is  less  than  a 
given  cutoff  value.  An  intuitively  appealing  cutoff  value  is  a  multiple  of 
the  equilibrium  0-0  distance  in  stable  ozone.  The  value  should  be  chosen 
large  enough  so  that  no  multiple  crossings  of  this  region  occur.  If  there  are 
no  multiple  crossings  the  lifetime  of  the  complex  is  simply  the  time  between 
when  the  incoming  oxygen  atom  first  crosses  the  cutoff  value  and  when  it 
recrosses.  Since  no  further  recrossing  is  assumed  to  occur,  the  colliding 
pair  separate  asymptotically  after  recrossing. 
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The  classical  procedure  for  the  computation  is  therefore  the  following: 
a)  collision  pairs  of  oxygen  atoms  and  oxygen  molecules  are  chosen  from  the 
appropriate  Boltzmann  and  flux  distributions,  b)  for  collisions  in  which  the 
oxygen  atom  separation  decreased  below  the  cutoff  value,  the  value  P(Ey,Ej^,£^) 
is  set  equal  to  unity  and  the  time  between  the  initial  crossing  and  recrossing 
is  recorded;  otherwise  P(  Ey ,  Ej^,  E^. )  equals  zero  and  c)  the  procedure  is  summed 
over  all  impact  parameters. 

The  choice  of  complex  formation  cutoff  parameter  does  not  seem  critical. 
This  i''  because  for  two  values  of  the  cutoff  parameter,  the  time  spent  trav¬ 
ersing  the  distance  between  should  be  small  compared  to  the  lifetime  of  the 
complex.  Some  computer  experimentatir>r«  is  necessary.  As  mentioned  earlier  it 
is  necessary  to  choose  a  cutoff  parameter  such  that  multiple  crossings  do  not 
occur.  If  they  occur  the  computation  can  be  substantially  complicated.  A 
reasonable  remedy  to  that  problem  would  be  an  asymmetric  set  of  cutoff  param¬ 
eters.  For  example,  a  smaller  value  in  the  entrance  channel  for  the  critical 
oxygen  atom  distance  and  a  larger  value  in  the  exit  channel  would  likely  elim¬ 
inate  multiple  crossings. 

Computer  programs  for  the  calculation  of  K(T)  have  been  completed  and 
tested  for  accuracy.  However  calculations  have  not  yet  been  undertaken. 

The  next  step  in  the  calculation  of  the  ozone  recombination  rate  is  the 
determination  of  the  collisional  stabilization  probability  of  ozone  complexes 
by  argon  atoms.  This  calculation  may  be  carried  out  by  fairly  standard  lonte- 
Carlo  classical  trajectory  techniques  except  for  the  initial  ozone  complex 
conditions. 

The  ozone  complex  initial  conditions  may  be  chosen  using  either  of  two 
procedures:  1)  random  sampling  from  appropriate  distributions  or  2)  using 

coordinates  of  complexes  determined  from  the  calculation  of  K(T).  The  first 
method  has  the  advantage  of  being  self-contained,  i.e.,  it  does  not  depend  on 
the  previous  calculations.  However  it  has  the  very  serious  difficulty  of  not 
readily  relating  to  the  phase  space  of  collisionally  formed  complexes.  One 
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may  readily  choose  initial  conditions  from  the  appropriate  distribution  func¬ 
tions  to  match  the  total  average  energies  and  angular  momentum  of  collision- 
ally  formed  complexes  (assuming  them  to  be  known;  such  averages  require  fairly 
extensive  complex  lifetime  calculations).  However  it  seems  impossible  to  show 
(without  extensive  additional  calculation)  that  the  lifetime  and  phase  space 
distributions  are  the  same  for  collisionally  formed  complexes  as  those  chosen 
by  a  Monte-Carlo  phase  space  sampling  procedure.  Therefore  it  is  much  more 
efficient  to  use  the  coordinates  of  complexes  computed  directly  from  O-O2  col¬ 
lisions  as  initial  conditions  for  complex  stabilization  computations. 

The  procedure  to  be  followed  is  as  follows:  a)  random  initial  conditions 
for  O-O2  collisions  are  chosen,  b)  the  trajectory  is  integrated  and  the  forma¬ 
tion  of  a  complex  is  determined  using  the  aforementioned  criteria,  and  c)  if  a 
complex  is  formed  the  coordinates  are  recorded  and  stored  as  initial  condi¬ 
tions  for  the  stabilization  calculation.  In  order  not  to  waste  too  much  com¬ 
putational  time  with  very  short  lived  complexes  that  cannot  contribute  effec¬ 
tively  to  the  ozone  production  rate,  a  cutoff  value  for  the  complex's  lifetime 
must  be  chosen.  An  appropriate  criterion  is  to  use  the  time  for  an  argon  atom 
to  move  from  the  asymptotic  regime  to  its  distance  of  closest  approach  in  the 
stabilization  calculation,  i.e.,  a  lifetime  long  enough  so  that  the  complex 
does  not  fall  apart  before  the  argon  atom  crrives.  Some  computational  experi¬ 
mentation  is  clearly  warranted. 

A  practical  advantage  of  this  procedure  is  that  the  initial  complex  coor¬ 
dinates  can  be  generated  without  significant  additional  calculation  during  the 
calculation  of  K(T).  This  may  be  quite  substantial  since  th;:se  calculations 
will  be  quite  computer  Intensive. 

The  final  step  in  the  simulation  of  the  recombination  process  is  to  com¬ 
pute  the  redissociation  probability  and  relaxation  of  newly  formed  (nascent) 
ozone  molecules.  For  purposes  of  computing  new  reaction  cates  only  the  redis¬ 
sociation  probability  is  necessary.  However  it  is  of  interest  to  examine 
relaxation  of  the  recombined  ozone  molecules  in  order  to  predict  the  radiative 
characteristics  of  such  a  recombining  gas.  The  CCX^HISE  facility  is  deaigned 
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to  neasura  radiation  from  recombining  ozone.  However,  this  is  not  possible  in 
a  collision  free  environment  and  therefore  measured  radiative  levels  have  been 
affected  by  some  degree  of  relaxation.  It  is  therefore  important  to  gain 
insight  into  the  relaxation  processes.  This  information  will  also  be  useful 
in  understanding  processes  involving  ozone  recombination  in  the  upper  atmos¬ 
phere  where  the  rate  of  relaxation  of  the  nascent  distribution  may  differ  sub¬ 
stantially  from  that  in  the  COCHISE  experiments. 

The  appropriate  simulation  of  the  relaxation/redissociation  processes 
employs  the  multicollision  simulation  of  the  Boltzmann  equation.  This  simula¬ 
tion  technique  has  been  discussed  earlier^  so  only  a  brief  outline  will  be 
given. 

The  Boltzmann  equation  for  the  recombining  ozone  in  argon  mixture  may  be 
written  as; 

3[n  (t)  f  (E.T.t)] 

3 _ 3 _  -  -  Ln  (t)  f  (E,T,t)  +  R  (8) 

at  3  3 

where  ng^^t)  is  the  concentration  of  ozone  molecules,  fo2(E»T,t)  is  the  phase 
space  distribution  function  of  ozone  molecules,  L  is  the  Boltzma.tu  collision 
operator,  R  is  a  reactive  sources  term  and  E  and  T  refer  to  the  ozone  Internal 
coordinates  and  bath  temperature  respectively.  The  reactive  source  term  is 
computed  from  the  preceding  calculation  and  is  the  so-called  one  way  flux 
production  term.^ 

Using  the  notation  of  Ref.  3  the  time  dependent  rate  coefficient  may  be 
written  as 


r 


n-1 


I 

j-1 


/  Xj  Ld  dV  dE/  [n^n^J 


(9) 


where 


Xj  -  (1  -  LT)^Xg  ,  (10) 

Xo  -  TR(T,E)nQnQ 
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and  is  the  part  of  the  Boltzmann  collision  operator  leading  to  dissocation. 
The  one-way  flux  distribution  of  recombined  ozone  molecules  is  given  by  Xq  and 
is  computed  from  the  prior  stabilization  calculation.  The  ozone  distribution 
that  results  from  j  bath  collisions  with  is  Xj  and  is  the  primary  quantity 
computed  in  the  multicollisional  calculation. 

After  a  sufficient  number  of  collisions  enough  relaxation  occurs  so  that 
the  redissociation  probability  of  a  recombined  molecule  becomes  zero.  The 
rate  constant  then  achieves  a  steady  state  value.  Previous  calculations  on 
diatomic  systems^  indicate  between  20  and  50  collisions  are  sufficient  for 
achieving  a  steady  state  recombination  rate.  At  that  point  recombined  mole¬ 
cules  are  still  highly  excited  but  far  enough  down  the  potential  so  that  col- 
lislonal  dissociation  does  not  occur.  It  is  anticipated  that  similar  numbers 
of  collisions  will  be  required  to  describe  ozone  redissociation.  Further 
relaxation  of  the  recombining  system  is  followed  by  simply  computing  more 
sequential  collisions. 

The  methodology  described  above  provides  a  systematic  computational 
scheme  for  the  calculation  of  the  ozone  recombination  rate  constant.  Compu¬ 
tational  limitations  have  prevented  such  a  program  from  being  completed.  A 
limited  set  of  calculations  were  completed  on  the  formation  of  long-lived  com¬ 
plexes  and  their  subsequent  dissociation  and  relaxation.  These  calculations 
were  quite  extensive  and  systematically  studied  the  dynamics  of  complex  forma¬ 
tion.  They  are  reported  in  the  following  sections  of  this  report. 
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3.  CALCULATIONS 


Classical  trajectory  calculations  were  performed  for  the  formation  of 
long-lived  complexes  of  ozone  by  collision  of  0-atoms  with  O2  molecules  and 
stabilization  of  formed  complexes  by  argon  atoms.  The  calculations  were 
purely  classical  with  the  Initial  diatom  vibrational  energy  chosen  from  the 
classical  distribution  function.  All  calculations  Involving  O-O2  collisions 
were  performed  using  the  Murrell-Sorbie-Varandas  ozone  potential  energy  sur¬ 
face.^  Calculations  of  Ar-03  collisions  were  performed  using  the  Murrell- 
Sorbie-Varandas  potential  energy  surface  for  the  ozone  portion  and  additive 
Lennard-Jones  6-12  potentials  between  the  argon  atom  and  each  oxygen  atom.^*^ 

The  dependence  of  long-lived  complex  (>1.2  x  10~^^s)  formation  on  rela¬ 
tive  translational  and  diatom  energies  was  computed.  Calculations  were  per¬ 
formed  at  relative  translational  energies  of  0.239,  0.478,  0.717,  0.956, 

1,194,  1.433,  1.672,  1.911,  2.150,  2.389,  and  2.628  kcal/mol  and  diatom  rota¬ 
tional  energies  of  0.0,  0.239,  0.478,  0.717,  and  0.956  kcal/mol  (for  each 
translational  energy).  Higher  values  of  relative  translational  energy  calcu¬ 
lations  (2.867,  3.106,  3.344,  3.583,  and  4.061  kcal/mol)  were  performed  at 
Initial  diatom  rotational  energy  equal  to  zero.  At  translational  energies 
below  or  equal  to  2.150  kcal/mol  batches  of  at  least  750  trajectories  were 
computed  for  each  value  of  Initial  relative  translational-diatom  rotational 
energy.  At  higher  translational  energies  batches  of  500  trajectories  were 
computed.  A  total  of  approximately  42,000  trajectories  were  computed.  At  low 
translational  and  rotational  energies  the  computing  times  became  very  large 
due  to  the  relatively  long  lifetimes  of  the  O3*  collision  complexes. 

Calculations  were  performed  In  order  to  determine  stabilization  proba¬ 
bilities  of  complexes  by  collision  with  argon  atoms.  A  200  K  argon  atom 
thermal  bath  was  employed  In  all  stabilization  calculation.  The  calculations 
were  designed  to  Investigate  variations  In  stabilization  probability  with  com¬ 
plex  energy  and  Initial  formation  conditions.  It  Is  convenient  to  character¬ 
ize  the  complexes  by  the  Initial  relative  translational  energy  and  the  Initial 
diatom  rotational  energy,  (E^,E{^).  Thirteen  cases  were  used  for  stabilization 
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calculations.  These  included:  (0.239,0.0),  (0.478,0.0),  (0.717,0.0), 
(0.956,0.0),  (1.672,0.0),  (2.389,0.0),  (0.239,0.239),  (0.478,0.239), 
(0.717,0.239),  (0.956,0.239),  (1.433,0.239),  (0.478,0.478),  and  (0.956,0.478). 
Six  hundred  or  more  trajectories  were  computed  for  each  stabilization  calcu¬ 
lation.  Since  the  number  of  long-lived  trajectories  determined  from  the  0-02 
trajectories  was  always  less  than  600,  each  of  the  initial  conditions  was  used 
several  times  in  the  stabilization  calculations.  For  the  higher  total  energy 
cases  very  few  long-lived  complexes  were  found  and  the  statistics  of  stabili¬ 
zation  may  be  suspect.  A  total  of  approximately  80,000  trajectories  were  com¬ 
puted  for  these  calculations. 

Multicollisional  calculations  were  carried  out  to  compute  redissociation 
probabilities  and  energy  relaxation  of  nascent  recombined  ozone.  The  cases 
considered  were:  (0.239,0.0),  (0.478,0.0),  (0.478,0.239),  (0.717,0.0), 
(0.478,0.478),  (0.717,0.239),  (0.956,0.0),  (0.956,0.239)  and  (0.956,0.478). 

In  each  case  the  initial  batch  of  stabilized  ozone  molecules  were  subjected  to 
a  sequence  of  50  collisions  with  a  200  K  ozone  thermal  bath.  Initial  batches 
varied  from  120  to  150  ozone  molecules.  These  calculations  required  the  com¬ 
putation  of  more  than  50,000  trajectories. 
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4.  RESULTS  AND  DISCUSSION 


In  this  section  we  discuss  the  results  of  the  calculations  for  the  forma¬ 
tion,  stabilization  and  subsequent  relaxation  of  long-lived  complexes.  Our 
goal  Is  to  Interpret  the  results  of  the  formatlon/stablllzatlon  calculations 
In  terms  of  the  energy  transfer  characteristics  found  In  Ref.  7.  The  results 
of  the  relaxation  calculations  will  be  used  as  an  aid  to  the  Interpretation  of 
the  COCHISE  experimental  results  and  atmospheric  processes. 

4. 1  Long-Lived  Complex  Calculations 

Probabilities  of  long-lived  complex  formation  as  a  function  of  Initial 
rela.tlve  oxygen  atom-oxygen  molecule  translational  energy  Is  shown  In  Figure  1 
for  the  case  of  an  Initially  non-rotating  diatom.  As  Is  readily  seen,  complex 
formation  probabilities  decrease  rapidly  with  Increasing  translational  energy. 


Figure  1.  Probability  of  Long-Lived  Complex  Formation  as  a  Function  of 
Relative  Translational  Energy 
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Even  at  the  lowest  translational  energy  the  probability  of  forming  complexes 
of  lifetime  greater  than  1.2  x  10~^^s  is  quite  small.  At  relative  transla¬ 
tional  energies  than  about  2  kcal/mol  there  was  generally  one  (or  no)  trajec¬ 
tory  that  produced  a  long-lived  complex.  A  complete  list  of  computed  proba¬ 
bilities  is  listed  in  Table  1. 


Table  1.  Computed  Probabilities  for  the  Formation  of  Long-Lived  Complexes 


E^/Et 

0.0 

(kca 1/mol) 

0.239 

(kcal/mol) 

0.478 
(Icca  1/mol) 

0.717 
(kca 1/mol) 

0.956 

(kcal/mol) 

0.239 

0.062 

0.002 

0.001 

0.00 

0.00 

0.478 

0.044 

0.010 

0.008 

0.003 

0.001 

0.717 

6.033 

0.009 

0.007 

0.00 

0.001 

0.956 

0.026 

0.011 

0.008 

0.005 

0.005 

1.194 

0.011 

0.005 

0.005 

0.005 

0.003 

1.433 

0.009 

0.007 

0.005 

0.004 

0.004 

1.672 

0.007 

0.0)1 

0.001 

0.001 

0.001 

1.911 

0.00 

0.003 

0.004 

0.002 

0.001 

2.150 

0.004 

0.001 

0.003 

0.001 

0.003 

2.389 

0.002 

0.006 

0.002 

0.002 

0.00 

2.628 

0.002 

0.00 

0.002 

0.006 

0.00 

2.867 

0.006 

0.00 

0.002 

3.106 

0.002 

3.344 

0.002 

3.583 

0.002 

4.061 

0.00 

4.778 

0.00 

The  values  listed  in  Table  1  clearly  show  that  increasing  Initial  diatom 
rotational  energy  decreases  the  probability  of  long-lived  complex  formation. 
This  is  most  clear  at  low  translational  energies.  At  the  lowest  translational 
energies  and  highest  rotational  energies  the  diatoms  are  relatively  hot  and 
readily  transfer  some  of  their  rotational  energy  to  relative  translational 
energy  thereby  reducing  the  complex  formation  probability.  A  complex  can  only 
be  formed  when  some  of  the  relative  atom-diatom  translational  energy  is  trans¬ 
ferred  into  Internal  motion  of  the  complex.  The  incoming  atom  becomes  trapped 
when  enough  energy  is  lost  so  that  it  cannot  overcome  the  attractive  forces  in 
the  complex.  Excess  rotational  energy  is  easily  transferred  from  the  diatom 
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so  that  a  rotatlonally  excited  diatom  will  have  little  chance  of  forming  a 
long-lived  complex.  Also  at  higher  translational  energies  the  amount  of 
energy  that  must  be  transferred  Into  Internal  complex  modes  Is  greater  and 
that  reduces  the  probability. 

At  Intermediate  translational  energies  (above  about  1  Iccal/mol)  the 
effect  of  rotational  energy  Is  present  but  not  so  powerful.  This  reduced 
effect  may  simply  be  seen  as  due  to  the  fact  that  at  higher  translational 
energies  the  same  amount  of  rotational  energies  produces  less  hot  diatoms. 
Above  approximately  1.5  kcal/mol  relative  translational  energy  very  few  long- 
lived  trajectories  were  computed  and  the  statistical  uncertainties  of  the 
results  listed  In  Table  1  are  large. 

An  In^llcatlon  of  the  rapid  decrease  of  the  formation  of  long-lived 
complexes  Is  that  the  equilibrium  constant,  K(T)  will  show  a  strong  decrease 
with  Increasing  temperature.  This  Is  simply  because  the  anticipated  lifetime 
of  complexes  formed  at  higher  energies  will  be  shorter.  A  systematic  discus¬ 
sion  of  the  temperature  variation  awaits  computation  of  K(T). 

Some  characterization  of  the  long-lived  complexes  Is  useful  for  under¬ 
standing  their  subsequent  colllslonal  behavior.  The  energy  of  the  complex  Is 
of  course  the  sum  of  the  Initial  relative  translational  and  diatom  Internal 
energies.  Therefore  complexes  formed  at  higher  relative  translational  and 
diatom  rotational  energies  (assuming  constant  diatom  vibrational  energy)  will 
have  to  lose  more  energy  by  collision  with  third  bodies  to  form  a  stable  ozone 
molecule  than  those  formed  at  lower  Initial  energies.  A  naive  assessment 
would  predict  a  decreasing  stabilization  probability  with  Increasing  long- 
lived  complex  energy.  However  another  factor  must  be  taken  Into  considera¬ 
tion,  the  total  angular  momentum  of  the  complex.  In  earlier  work^  on  Ar-03 
colllslonal  energy  transfer,  the  trlatomlc  angular  momentum  was  found  to  be 
the  controlling  factor  In  both  the  magnitude  and  the  sign  of  energy  trans¬ 
ferred  between  highly  excited  ozone  and  argon.  In  particular,  molecules  with 
low  angular  momentum  are  found  to  gain  energy  by  collision  despite  very  high 
Internal  energy.  The  significance  of  this  observation  for  stabilization  of 
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complexes  Is  that  energy  relaxation  to  form  stable  molecules  will  be  favored 
by  higher  complex  angular  momentum.  However,  higher  complex  angular  momentum 
Implies  In  general  complexes  formed  at  high  total  energy  requiring  larger 
energy  transfers  for  stabilization. 

The  variation  In  the  average  long-lived  complex  angular  momentum  as  a 
function  of  total  energy  of  the  complex  Is  shown  In  Figure  2  .  Three  cases 
are  shown:  complexes  formed  by  collisions  with  non-rotating  diatoms,  formed 
by  collisions  with  diatoms  having  0.239  kcal/mol  Initial  rotational  energy  and 
for  Initial  diatom  rotational  energies  of  0.478  kcal/mol.  Since  at  higher 
collision  energies  the  number  of  computed  complexes  was  very  small,  the 
behavior  shown  In  Figure  2  Is  not  reliable  from  a  quantitative  point  of  view 
due  to  statistical  uncertainties.  The  Increase  of  complex  angular  momentum 
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Figure  2.  Angular  Momentum  of  Ozone  Complexes  as  a  Function  of  Complex 
Energy.  Curves  denote  Initial  diatom  rotational  energy  as 
follows:  A,  0.0  kcal/mol;  B,  0.239  kcal/mol;  C,  0.478  kcal/mol. 
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with  complex  energy  is  readily  seen  for  the  non-rotating  case.  The  Murrell- 
Sorble-Varandas  ozone  potential  surface  Is  quite  directional  so  that  a  complex 
Is  formed  only  when  the  three  atoms  are  In  a  narrow  range  of  configurations. 
This  fixes  the  allowed  range  of  Intact  parameters  for  complex  formation.  As 
the  relative  translational  energy  Increases  the  orbital  angular  momentum  for 
this  range  of  Impact  parameters  Increases.  The  orbital  angular  momentum  Is 
converted  Into  complex  angular  momentum. 

For  the  rotating  cases  the  situation  Is  more  difficult  since  the  vector 
coupling  of  orbital  and  diatom  angular  momentum  determine  the  total  complex 
angular  momentum.  However  the  computed  values  show  a  clear  Increase  In  the 
complex  angular  momentum  with  Increasing  complex  energy.  In  the  next  section 
the  results  of  stabilization  calculations  will  be  presented. 

4.2  Stabilization  Calculations 


Calculations  were  performed  to  determine  the  stabilization  probability  of 
long-lived  complexes  by  collision  with  argon  atoms.  An  argon  atom  bath  of 
200  K  Is  used  for  all  the  calculations.  Calculations  were  done  using  coordi¬ 
nates  from  a  previously  determined  long-lived  complex  as  the  Initial  trlatom 
coordinates.  Relative  atom-trlatom  coordinates  were  chosen  randomly  from 
appropriate  Boltzmann  distributions. 

Computed  probabilities  of  stabilization  verses  complex  energy  are  shown 
In  Figure  3.  Each  set  of  connected  points  are  for  complexes  formed  from 
diatom  collisions  with  a  given  value  of  Initial  rotational  energy  (A: 0.0; 
B:0.239  hcal/mol  and  C:0.478  kcal/mol,  initial  diatom  rotational  energy).  The 
most  obvious  feature  of  the  computed  stabilization  probabilities  Is  that  there 
is  no  dlscernable  dependence  on  complex  energy.  From  the  earlier  discussion 
this  may  be  attributed  to  the  near  balancing  of  two  opposite  effects.  Higher 
energy  complexes  require  larger  energy  transfers  for  stabilization  and  should 
have  lower  stabilization  probabilities.  On  the  other  hand  higher  energy  com¬ 
plexes  have  higher  angular  momentum  which  facilitates  energy  transfer  from  the 
trlatom.  These  two  effects  seem  to  cancel  one  another  leading  to  no  obvious 
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dependence  for  stabilization  probability  on  complex  energy  (for  the  range  of 
energies  and  angular  momenta  studied) .  Again  these  conclusions  should  be  tem¬ 
pered  by  the  fact  that  the  sample  sizes  were  exceedingly  limited  due  to  the 
low  probability  of  forming  long-lived  complexes  at  higher  energies. 

The  implication  for  the  general  computation  of  recombination  rates  for 
ozone  is  that  the  temperature  dependence  arises  primarily  from  the  dependence 
of  K(T)  on  temperature.  Higher  energy  complexes  would  be  more  prevalent  at 
higher  temperatures.  On  the  basis  of  the  present  calculation  stabilization 
does  not  seem  to  depend  strongly  on  complex  energy.  We  have  not,  however, 
examined  stabilization  probability  as  a  function  of  argon  atom  bath  temper¬ 
ature.  Earlier,  systematic  studies  of  energy  transfer  between  highly 
vlbratlonally-rotationally  excited  ozone  and  argon  indicates  that  the  energy 
transfer  is  most  controlled  by  the  triatom  angular  momentum.  Stabilization 
calculations  at  different  bath  temperatures  are  clearly  warranted.  On  the 
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Figure  3.  Collisional  Stabilization  Probability  of  Ozone  Complexes  by  Argon 
Atoms  as  a  Function  of  Complex  Energy.  Curves  are  labeled  as  in 
Figure  2. 
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basis  of  the  limited  data  from  the  present  calculations,  stabilization  proba¬ 
bilities  will  be  weakly  dependent  on  temperature.  Complex  formation  is  (as 
shown  earlier)  strongly  dependent  on  the  colliders'  energy. 

We  now  consider  the  behavior  of  newly  recombined  ozone  molecules  in  an 
argon  thermal  bath. 

4.3  Multicolllslonal  Calculations 


Calculations  were  performed  to  simulate  the  relaxation  of  newly  recom¬ 
bined  ozone  molecules  in  a  200  K  argon  bath.  The  calculations  were  done  using 
the  methodology  described  in  Ref.  3.  The  evolution  of  each  batch  of  nascent 
ozone  molecules  was  followed  for  50  collisions. 

The  primary  quantities  of  interest  from  these  calculations  are  the  disso¬ 
ciation  probability  of  recombined  ozone  molecules  and  the  rate  of  energy 
relaxation.  First  we  discuss  dissociation  probabilities. 

Computed  redissoclatlon  probabilities  are  listed  in  Table  2.  The  nascent 
ozone  ensembles  are  labeled  by  the  energies  of  the  unbound  complexes  from 
which  they  are  colllsionally  stabilized.  Generally  the  computed  redissocla¬ 
tlon  probabilities  are  quite  large,  varying  from  59  to  15  percent.  In  a  com¬ 
plete  calculation  of  a  thermal  recombination  rate  these  values  (when  properly 
averaged)  would  be  used  in  Eq.  (1)  and  would  constitute  a  significant  correc¬ 
tion  to  the  so-called  one-way  flux  recombination  rate  constant. 

An  obvious  but  curious  trend  in  the  values  listed  in  Table  2  is  that  the 
molecules  formed  from  higher  energy  complexes  have  lower  redissociation  proba¬ 
bilities.  To  understand  this  behavior  we  consider  the  energies  of  the  nascent 
ozone  molecules  formed  by  stabilization  of  the  ozone  complex.  These  values 
are  listed  in  Table  3.  It  is  quite  clear  from  the  listed  values  that  the 
higher  energy  complex  produce  lower  energy  molecules.  This  is  readily  seen  as 
due  to  the  larger  angular  momentum  of  the  higher  energy  complexes.  The  effect 
of  molecular  angular  momentum  on  energy  transfer  was  discussed  earlier. 
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Table  2.  Dissociation  Probabilities  of  Nascent  Ozone  Ensembles  in  a  2C0  K 
Argon  Atom  Bath.  Ensembles  are  labeled  by  the  energies  of  the 
long-lived  complexes  from  which  they  are  derived. 


Et/Ej 

0.00 

(kcal/mol) 

0.239 

(kcal/mol) 

0.478 

(kcal/mol) 

0.239 

0.59 

0.478 

0.37 

0.44 

0.28 

0.717 

0.26 

0.30 

0.956 

0.22 

0.15 

0.20 

Table  3.  Average  Nascent  Ozone  Ensemble  Energy  Below  Dissociation 
Formed  by  Stabilization  of  Complexes  of  Varying  Energy. 
All  energies  in  units  of  kcal/mol. 


0.00 

0.239 

0.478 

Et/Ej 

(kcal/mol) 

(kcal/mol) 

(kcal/mol) 

0.239 

-0.599 

0.478 

-0.861 

-0.812 

-0.975 

0.717 

-0.896 

-0.964 

0.956 

-1.145 

-1.474 

-1.397 

Zt  is  reasonable  that  lower  energy  molecules  will  suffer  fewer  dissoci¬ 
ation  collisions.  Therefore  the  trend  observed  in  Table  2  for  redissociation 
probabilities  is  simply  that  higher  energy  complexes  produce  lower  energy 
molecules  that  have  lower  redissociation  probabilities. 


The  multlcolllsional  energy  relaxation  of  recombined  ozone  molecules  is 
now  considered.  The  average  total  internal  energy  of  a  relaxing  ensemble  of 
recombined  ozone  molecules  is  shown  as  a  function  of  collision  with  the  200  K 
argon  bath  in  Figures  4  through  12.  The  behavior  of  each  case  is  fairly  simi¬ 
lar,  relaxing  to  about  6  kcal/mol  below  the  dissociation  limit  of  the  ozone 
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ENERGY  BELOW  DISSOCIATION  (kcal/noU  T  ENERGY  BELOW  DISSOCIATION  (ltc«1/iM>1) 


-1.00 


NUMBER  OF  COLLISIONS 


A-3821 


Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 
Sequential  Collisions  with  Argon  Atom  Thermal  Bath.  Ozone 
ensemble  Is  labeled  by  the  energies  of  the  long-lived  complexes 
from  which  they  are  derived.  Calculation  Is  for  (0.01,0.0) 

case. 


Figure  5.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.02,0.0)  case. 
See  Figure  4  for  notation. 
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Figure  6.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.03,0.0)  case. 
See  Fig;ire  4  for  notation. 


NUMBER  Of  COLLISIONS 


A-38Z4 


Figure  7.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.04,0.0)  case. 
See  Figure  4  for  notation. 
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ENERGY  BELOU  DISSOCIATION  (kcal/mol) 


Figure  8 


Figure  9 


NUMBER  OF  COELISIONS 


A-38^5 


.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 
Sequential  Collisions  with  Argon  Atom  Bath  for  (0.02,0.01)  case 
See  Figure  4  for  notation. 


NUMBER  OF  COLLISIONS 


A-3826 


Average  Energy  of  Ozone  Molecules  as  a  Function  of  Numoer  of 
Sequential  Collisions  with  Argon  Atom  Bath  for  (0.02,0.02)  case 
See  Figure  4  for  notation. 


Figure  10.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.03,0.01)  case. 
See  Figure  4  for  notation. 


Figure  1 1 .  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.04,0.01)  case. 
See  Figure  4  for  notation. 
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Figure  12.  Average  Energy  of  Ozone  Molecules  as  a  Function  of  Number  of 

Sequential  Collisions  with  Argon  Atom  Bath  for  (0.04,0.02)  case. 
See  Figure  4  for  notation. 

molecule  after  fifty  collisions.  Redlssociatlon  of  newly  formed  molecules 
occurs  only  at  early  times.  Very  little  redlssociatlon  occurs  beyond  the 
twentieth  sequential  collision.  The  duration  of  redlssociatlon  processes  Is  a 
function  of  the  argon  thermal  bath  energy.  Hence  It  Is  not  possible  to  extend 
this  conclusion  to  different  temperature  baths. 

The  energy  relaxation  appears  to  show  two  regimes:  1)  an  early  time 
regime  during  which  redlssociatlon  occurs  and  the  energy  relaxation  la  faster 
and  2)  a  long  time  regime  with  slower  energy  relaxation.  The  slightly  faster 
energy  relaxation  of  the  early  time  regime  may  be  attributed  to  redlssocla- 
tlon.  Redlssociatlon  selectively  removes  higher  energy  molecules  from  the 
ensemble  and  thereby  Increases  the  apparent  rate  of  relaxation.  This  sort  of 
behavior  has  been  observed  earlier  for  the  relaxation  of  recombined  bromine 
molecules  in  an  argon  thermal  bath.^ 
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The  overall  ozone  internal  energy  relaxation  Is  fairly  slow.  The  average 
initial  ensemble  internal  energies  are  in  each  case  between  0.5  and 
1.5  hcal/mol  below  the  dissociation  limit,  l.e.,  they  have  vibrational- 
rotational  energies  of  about  24.5  to  25.5  hcal/mol.  After  50  collisions  the 
molecules  are  still  very  vibrationally-rotatlonally  excited  with  energies  of 
about  19-20  kcal/mol.  This  slow  energy  relaxation  is  consistant  with  earlier 
single  collision  and  ensemble  relaxation  calculations. 

Histograms  of  internal  energy  distributions  are  shown  for  each  case  in 
Figures  13  through  21.  Each  figure  shows  the  initial  internal  distribution 
after  colllslonal  stabilization  of  the  ozone  complex  and  the  internal  energy 
distribution  after  fifty  sequential  collisions  with  the  200  K  argon  bath.  The 
behavior  in  each  case  is  similar  so  that  we  may  discuss  Figure  14  as 
characteristic  of  all  the  cases. 


OZONE  INTERNAL  ENERGY  (kcal/ntol)  A- 3830 


Figure  13.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  is  (0.01,0.0). 
See  Figure  4  for  notation. 
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020NE  INTERNAL  ENERGT  (kcal/mol)  A- 3831 

Figure  14.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  Is  (0.02,0.0) 
See  Figure  4  for  notation. 


020NE  INTERNAL  ENERGY  (kcil/nwll  A- 3832 

Figure  15.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  Is  (0.03,0.0) 
See  Figure  4  for  notation. 
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Figure  16.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  Is  (0.04,0.0). 
See  Figure  4  for  notation. 


Figure  17.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  Is  (0.02,0.01). 
See  Figure  4  for  notation. 
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OZONE  INTERNAL  ENERGY  (kcal/mol)  A-3835 

Figure  18.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  vith  Argon  Atom  Bath.  Case  shown  is  (0.02,0.02) 
See  Figure  4  for  notation. 


Figure  19.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  is  (0.03,0.01) 
See  Figure  4  for  notation. 
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Figure  20.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  is  (0.04,0.01) 
See  Figure  4  for  notation. 


OZONE  INTERNAL  ENERGY  (kcal/mol)  A-3838 


Figure  21.  Distribution  of  Recombined  Ozone  Total  Internal  Energies: 

a)  Nascent  Distribution  and  b)  After  Fifty  Sequential 
Collisions  with  Argon  Atom  Bath.  Case  shown  is  (0.04,0.02) 
See  Figure  4  for  notation. 


Most  of  the  newly  recombined  ozone  molecules  are  within  1  kcal/mol  of  the 
dissociation  limit.  The  lowest  energy  molecules  are  less  than  4  kcal/mol 
below  the  dissociation  limit.  Over  90  percent  of  the  nascent  distribution 
lies  within  2  kcal/mol  of  the  dissociation.  The  narrowness  of  the  nascent 
distribution  is  most  likely  a  function  of  the  bath  temperature  since  the  mag¬ 
nitude  of  collisional  energy  transfer  depends  on  the  relative  translational 
energy.  It  is  anticipated  that  at  higher  temperatures  than  200  K  a  broader 
nascent  distribution  would  be  found. 

After  fifty  sequential  collisions  the  ensemble  internal  distribution  has 
shifted  to  lower  energies  and  broadened  compared  to  the  nascent  distribution. 
The  shift  is  readily  seen  by  noting  that  very  few  molecules  remain  in  the 
range  of  0  to  4  kcal/mol  below  the  dissociation  limit  (Figure  14b).  This  is 
of  course  due  to  energy  transfer  from  the  molecules  to  the  thermal  bath  and 
redissociation  of  molecules  close  to  the  dissociation  limit.  The  distribution 
is  reasonably  symmetric  about  the  peak  range  of  5  to  7  kcal/mol  below  disso¬ 
ciation  and  has  a  width  of  about  9  kcal/mol.  The  lowest  energy  molecules  are 
9  to  10  kcal/mol  below  dissociation. 

The  above  behavior  is  characteristic  of  all  the  cases  studied.  The  ini¬ 
tial  distribution  is  fairly  narrow,  3  to  4  kcal/mol,  with  most  molecules 
within  1  kcal/mol  of  dissociation.  After  fifty  collisions  the  high  energy 
portion  of  the  distribution  is  depleted  with  very  few  molecules  in  the  0  to  3 
or  4  kcal/mol  energy  range.  The  distributions  are  broadened  with  widths 
approximately  8  to  10  kcal/mol  of  the  dissociation  limit.  It  should  be  noted 
that  within  approximately  fifteen  to  twenty  sequential  collisions  relaxation 
of  each  ensemble  is  very  similar  due  to  loss  of  initial  state  memory. 

In  either  the  upper  atmosphere  or  the  COCHISE  experiments  newly  recom¬ 
bined  ozone  molecules  will  undergo  collisions.  Estimates  based  on  the  charac¬ 
teristics  of  the  nozzle  expansion  in  the  COCHISE  apparatus  indicate  perhaps  as 
many  as  several  thousand  collisions  before  the  recombined  ozone  molecules 
enter  the  collision  free  zone.  The  present  calculations  may  be  extrapolated 
to  this  longer  time  regime.  The  procedure  indicates  that  the  recombined 
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molecules  will  have  relaxed  to  about  half  way  down  the  molecular  well  after 
several  hundred  collisions.  Since  these  collisions  occur  In  a  time  Interval 
short  compared  with  the  vibrational  radiative  lifetime,  It  Is  anticipated  that 
radiation  from  the  upper  levels  near  the  dissociation  limit  will  be  collision- 
ally  quenched  and  therefore  not  observed.  This  prediction  Is  consistent  with 
the  eiqjerlmental  observation  that  radiation  Is  seen  from  vibrational  levels 
only  up  to  V  -  5. 
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INTRODUCTION 


In  this  report  we  describe  the  results  of  a  multicollisional  study  of  the 
relaxation  of  vibrationally-rotationally  excited  ozone  molecules  in  a  bath  of 
argon  atoms.  The  primary  objective  of  the  study  is  to  systematically  study 
the  effects  of  the  ozone  initial  state  and  the  thermal  bath  on  the  relaxation 
of  highly  excited  ozone  molecules. 

The  relaxation  of  highly  excited  ozone  molecules  is  of  considerable 
significance  in  the  study  of  the  upper  atmosphere.  Ozone  molecules  are 
thought  to  be  formed  by  the  following  three-body  recombination  reaction; 

O  +  O2  +  M  -►  O2  ( V ,  ij )  +  M 

where  the  initially  recombined  ozone  molecule,  O^CvjJ),  is  formed  in 
vibrational-rotational  states  near  the  dissociation  limit.  The  newly  formed 
ozone  molecule  can  undergo  a  number  of  processes  Including:  1)  collisional 
dissociation  by  energy  transfer  into  the  molecule,  2)  collisional  relaxation 
to  a  lower  internal  energy  state,  and  3)  radiative  emission  of  an  infrared 
photon.  Additional  processes  might  include  reactions  with  various  atmospheric 
species.  These  processes  are  not  included  in  this  study. 

The  relative  rates  of  relaxation/dissociation  versus  the  emission  rate 
will  determine  the  mode  of  relaxation  of  highly  excited  ozone.  Considerable 
Interest  exists  in  understanding  the  ozone  radiative  contribution  and  colli¬ 
sional  rate  of  relaxation.  A  theoretical  understanding  of  these  processes  is 
also  necessary  for  proper  interpretation  of  the  ozone  recombination  and  emis¬ 
sion  experiments  conducted  at  AFGL. 

Calculations 


All  calculations  reported  here  were  performed  using  the  multicollisional 
classical  trajectory  technique. ^  The  computation  of  classical  trajectories 
has  become  standard  and  is  well  reviewed  by  Bunker^  and  Porter  and  Raff.^  The 
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multicollisional  classical  trajectory  procedure  is  discussed  in  Ref.  1.  We 
mention  only  calculational  procedures  specific  to  the  present  calculations. 

Calculations  were  performed  to  systematically  investigate  the  effect  of 
thermal  bath  temperature  on  the  relaxation/dissociation  of  highly  excited 
ozone.  Calculations  were  performed  for  initial  ensembles  of  ozone  molecules 
excited  25  kcal/mol  above  the  well  minimum  (above  1.5  kcal/mol  below  disso¬ 
ciation)  with  zero  molecular  angular  momentum  at  argon  atom  bath  temperatures 
of  200,  300,  400,  500,  750,  and  1000  K.  Initial  ensemble  sizes  varied  between 
150  (at  200,  300,  and  400  K  bath  temperatures)  and  200  (at  500,  750,  and 
1000  K  bath  tempertures ) .  In  the  200  K  bath  temperature  case,  the  evolution 
of  the  ensemble  was  followed  for  a  sequence  of  50  collisions.  In  all  other 
cases,  the  ensemble  was  subjected  to  a  sequence  of  100  bath  collisions. 

Several  calculations  were  performed  at  an  argon  atom  bath  temperature 
of  300  K  to  investigate  the  effect  of  initial  triatom  rotational  energy  on 
relaxation/dissociation  and  the  effect  of  total  triatom  internal  energy. 

Several  other  cases  were  also  considered:  1)  ozone  initially  excited 
26  kcal/mol  above  the  well  minimum  and  zero  rotational  energy,  2)  26  kcal/mol 
above  the  well  minimum  with  an  initial  triatom  rotational  energy  of  2  kcal/mol, 
3)  27  kcal/mol  internal  energy  with  zero  rotational  energy,  and  4)  27  kcal/mol 
internal  energy  with  1  kcal/mol  rotational  energy.  Note  that  the  total  initial 
internal  energy  of  27  kcal/mol  is  above  the  dissociation  limit.  On  the  em¬ 
ployed  potential  surface,  the  dissociation  limit  of  ozone  to  an  oxygen  molecule 
and  atom  is  26.3  kcal/mol.  Sequences  of  50  successive  collisions  were  computed 
for  cases  2-4  and  100  for  case  1. 

Calculations  were  performed  for  collisional  dissociation  of  ozone  mole¬ 
cules  dissociated  by  energy  transfer  from  the  argon  atom  thermal  bath.  The 
criterion  for  dissociation  was  baaed  on  the  distance  between  any  pair  of  oxygen 
atoms  exceeding  a  cutoff  distance.  This  cutoff  distance  was  taken  as  3.0A. 
Hence  in  several  cases  molecules  with  greater  than  the  dissociation  energy  were 
retained  for  further  collisions.  In  fact,  for  the  higher  initial  energy 
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ensembles,  the  average  ensemble  energy  exceeded  the  dissociation  limit.  The 
retention  of  these  unbound  or  quasibound  molecules  amounts  to  assuming  that 
their  lifetimes  are  very  long  and  that  they  will  not  dissociate  before 
undergoing  further  collisions.  In  other  computations  on  ozone  complexes  very 
long-lived  unbound  complexes  have  been  observed  and  some  test  calculations 
indicated  this  behavior  often  occurs  for  the  complexes  considered  here.  This 
assumption  however  implies  that  the  dissociation  probabilities  are  lower 
bounds.  More  than  150,000  Ar-03  trajectories  were  computed  in  this  effort. 

Potential  Energy  Surface 


All  calculations  were  performed  using  the  Murrell-Sorbie-Varandas^ 
potential  energy  surface  for  ozone  and  pairwise  additive  Lennard-Jones  6-12 
potentials  between  the  argon  atom  and  each  oxygen  atom.  The  details  of  the 
potential  energy  surface  are  given  elsewhere. 5 

Results  and  Discussion 


The  behavior  of  total  ozone  internal  molecular  energy  and  angular  momen¬ 
tum  as  a  function  of  number  of  collisions  with  the  argon  atom  thermal  bath  is 
shown  in  Figures  1  to  10  for  all  the  cases  considered.  We  first  consider  the 
results  shown  in  Figures  1  through  6.  These  results  are  for  ensembles  of 
ozone  molecules  vibrationally  excited  25  kcal/mol  above  the  ozone  well  minimum 
with  zero  initial  molecular  angular  momentum.  Each  figure  corresponds  to  a 
different  argon  atom  bath  temperature.  The  behavior  for  argon  atom  bath  tem¬ 
peratures  of  200,  300,  400,  500,  750,  and  1000  K  are  shown  in  Figures  1 
through  6,  respectively.  The  behavior  shown  in  Figure  2  for  a  300  K  argon 
atom  bath  is  qualitatively  similar  to  the  others.  Therefore,  we  discuss  it  in 
more  detail  as  representative  of  all  the  cases  shown. 

Initially  the  ozone  molecular  angular  momentum  in  the  ensemble  is  set 
equal  to  zero.  This  corresponds  to  non-rotating  molecules.  Collisions  with 
the  bath  rapidly  transfer  angular  momentum  into  the  molecule.  This  process 
occurs  within  approximately  5  to  10  molecule-bath  collisions.  After  the 
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Figure  1.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Initial  total  ozone  internal  energy  equals  25  kcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  200  K. 


Figure  2. 


Variation  of  Average  Ozone  Internal  Energy  and  Molecuxar  Angular 
Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Initial  total  ozone  internal  energy  equals  25  kcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  300  K. 
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Figure  3.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentvim  with  Nvunber  of  Sequential  Collisions  with  Argon  Atom  Bath 
Initial  total  ozone  internal  energy  equals  25  hcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  400  K. 


Figure  4.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath 
Initial  total  ozone  internal  energy  equals  25  kcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  500  K. 
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AVERAGE  OZONE  ENTERNAL  ENERGY  u,  AVERAGE  OZONE  ENTERNAL  ENERGY 

(ABOVE  WELL  MINIMUM,  kcal/mol)  •  (ABOVE  WELL  MINIMUM,  kcal/mol) 


Variation  of  Average  Ozone  Internal  E:nergy  and  Molecular  Angular 
Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath 
Initial  total  ozone  internal  energy  equals  25  kcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  750  K. 


Figure  6. 


NUMBER  OF  COLLISIONS  WITH  ARGON  BATH 


A-4008 


Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 
Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath 
Initial  total  ozone  internal  energy  equals  25  kcal/mol;  initial 
ozone  angular  momentum  equals  zero;  bath  temperature  equals  1000  K 
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initial  transient  period,  the  molecular  angular  momentum  varies  only  slowly. 
(Fluctuations  in  angular  momentum  or  energy  are  most  probably  the  result  of 
using  finite  size  ensembles  rather  than  physical  processes.)  As  noted  in 
earlier  work  on  single  collision  Ar-O^  dynamics,  rotational  energy  is  trans¬ 
ferred  readily  to  relative  translational  energy.^  Thus  the  value  of  the 
molecular  angular  momentum  reflects  the  argon  atom  bath  temperature.  The 
rapidness  of  rotational-translational  exchange  is  verified  by  the  few, 

5  to  10,  collisions  required  to  bring  the  ozone  ensemble  angular  momentum  to 
its  near  steady  state  value. 

The  behavior  of  the  ozone  total  Internal  energy  is  slightly  more  compli¬ 
cated,  Two  regimes  are  evident.  During  the  first  few  collisions,  i.e., 

1  to  10,  energy  is  pumped  into  the  molecule  as  the  molecular  angular  momentum 
increases.  This  occurs  despite  the  fact  that  the  molecules  are  very  highly 
excited.  Similar  behavior  has  been  seen  earlier.^  After  the  initial 
transient,  the  molecular  energy  relaxes  monotonically . 

As  noted  earlier  the  ozone  energy  relaxation  is  qualitatively  similar 
over  the  temperature  range  considered.  The  two  most  significant  variations 
with  temperature  are  the  increase  in  energy  of  the  ozone  during  the  initial 
transient,  and  the  steady  state  value  of  the  ozone  molecular  angular  momentum. 

The  amount  of  energy  pumped  into  the  molecule  as  the  molecular  angular 
momentum  is  built  up  to  its  steady  state  value  increases  monotonically  with 
the  argon  thermal  bath  temperature.  At  a  bath  temperature  of  200  K,  the 
Increase  in  ensemble  average  energy  is  about  0.7  kcal/mol.  At  a  bath 
temperature  of  1000  K,  this  value  increases  to  1.9  kcal/mol.  This  variation 
will  become  significant  when  we  consider  dissociation.  In  the  1000  K  bath 
case,  the  average  ozone  molecular  energy  reaches  26.9  kcal/mol.  This  value  is 
above  the  dissociation  limit  of  ozone  on  the  employed  potential  surface 
(26.3  kcal/mol). 

Closely  related  to  transient  energy  increase  is  the  value  of  the  steady 
state  molecular  angular  momentum.  Since  molecular  rotational  energy  and 
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relative  translational  energy  are  exchanged  relatively  rapidly,  the  molecular 
angular  momentum  tends  to  reflect  the  thermal  bath  temperature.  The  steady 
state  value  at  200  K  is  about  20  while  at  1000  K  its  value  is  about  45 
Since  rotational  energy  scales  as  the  square  of  angular  momentum,  these  values 
scale  as  the  increase  in  relative  translational  energy  of  the  argon  atom 
thermal  bath  as  the  temperature  is  increased  from  200  to  1000  K. 

The  rate  of  molecular  internal  energy  relaxation  is  seen  from  Figures  1 
through  6  to  increase  slightly  as  the  bath  temperature  is  increased.  The 
interpretation  of  this  variation  is  complicated  by  the  fact  that  the  ensemble 
is  simultaneously  undergoing  relaxation  and  dissociation.  Dissociation  is 
biased  towards  higher  energy  molecules.  Therefore,  without  further  study,  we 
may  not  infer  a  faster  energy- transfer  relaxation  rate  at  higher  temperatures. 
In  fact  the  rate  probably  slows  down.  The  variation  shown  in  Figures  1 
through  6  show  the  rate  of  energy  relaxation  of  molecules  undergoing  energy 
transfer  and  dissocation. 

We  now  discuss  the  results  shown  in  Figures  7  through  10.  In  each  of 
these  cases  the  argon  atom  bath  is  300  K.  The  long-term  (past  ten  collisions 
with  the  bath)  variation  of  the  ozone  molecular  angular  momentum  is  controlled 
by  the  atom  bath  temperature.  This  is  due  to  relatively  facile  transfer  of 
rotational  and  relative  translational  energy.  This  was  mentioned  earlier  for 
cases  in  Figures  1  through  6.  Comparisons  of  Figures  7  through  10  with 
Figure  2  show  the  same  long-term  behavior  of  the  molecular  angular  momentum  at 
300  K  bath  temperature. 

Figures  7  and  9  show  the  internal  relaxation  of  initially  non-rotating 
ozone  molecules.  For  molecules  excited  26  )ccal/mol  above  the  well  minimvim 
(Figure  7)  the  behavior  is  virtually  identical  to  the  earlier  case  of  molecules 
excited  25  kcal/mol  in  a  300  K  atom  bath  (Figure  2).  For  molecules  excited 
27  kcal/mol  above  the  well  minimum  there  is  some  variation  in  relaxation 
behavior.  As  mentioned  earlier,  these  molecules  are  above  the  ozone  dissoci¬ 
ation  limit.  Therefore,  unless  stabilized  by  collisions  with  the  argon  bath. 
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Figure  7.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Bath  temperature  equals  300  K;  initial  ozone  rotational  energy 
equals  zero;  initial  ozone  internal  energy  equals  26  hcal/mol. 


Figure 


8.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 
Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Bath  temperature  equals  300  K;  initial  ozone  rotational  energy 
equals  2  kcal/mol;  initial  ozone  internal  energy  equals  26  kcal/mol 
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Figure  9.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Bath  temperature  equals  300  K;  Initial  ozone  rotational  energy 
equals  zero;  Initial  ozone  internal  energy  equals  27  kcal/mol. 


Figure  10.  Variation  of  Average  Ozone  Internal  Energy  and  Molecular  Angular 

Momentum  with  Number  of  Sequential  Collisions  with  Argon  Atom  Bath. 
Bath  temperature  equals  300  K;  Initial  ozone  rotational  energy 
equals  1  kcal/mol;  Initial  ozone  Internal  energy  equals  27  kcal/mol 
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they  will  dissociate.  The  short-time  increase  in  internal  energy  due  to 
buildup  of  rotational  energy  is  absent  in  this  case  (Figure  9).  This  may 
simply  be  attributed  to  the  fact  that  molecules  that  do  not  lose  energy  leave 
the  ensemble  by  dissociation.  In  the  earlier  cases,  Figures  1  through  6, 
molecules  can  acquire  total  internal  energies  above  the  dissociation  limit  by 
rotational  energy  transfer  from  the  bath.  In  most  of  those  cases,  however, 
the  molecular  vibrational  energy  is  insufficient  for  dissociation.  Dissoci¬ 
ation  occurs  most  effectively  through  vibrational  motion.  Thus  molecules  with 
energies  above  the  dissociation  limit  but  without  sufficient  vibrational 
energy  alone  for  dissociation  may  survive  for  long  periods.  This  occurs  at 
the  higher  bath  temperatures  25  kcal/mol  initial  excitation  and  for 
26  kcal/mol  initial  excitation  where  the  ensemble  average  interval  energy  goes 
above  the  dissociation  limit.  On  the  other  hand  when  the  initial  vibrational 
energy  exceeds  the  dissociation  limit  (Figure  9)  the  lifetime  of  the  unbound 
molecules  is  too  short  and  many  molecules  that  gain  energy  are  removed.  Thus 
the  average  energy  of  the  ensemble  does  not  show  the  characteristic  increase 
in  energy  at  short  times. 

The  effect  of  initial  molecular  rotation  on  the  relaxation  is  shown  by 
comparison  of  Figure  7  with  Figure  8  and  Figure  9  with  Figure  10.  Tha 
presence  of  initial  rotational  energy  reduces  or  eliminates  the  initial  short 
time  increase  in  the  ensemble  internal  energy.  This  is  clearly  shown  in 
Figure  8  and  10.  The  memory  of  initial  rotational  energy  disappears  after 
about  ten  collisions.  Energy  relaxation  after  that  point  is  the  same  for  all 
cases.  We  now  discuss  dissociation  in  these  relaxing  ensembles. 

The  computed  probability  of  dissociation  for  each  ensemble  is  presented 
in  Table  1.  The  listed  values  are  the  asymptotic  dissociation  probabilities 
(i.e.,  the  fraction  of  the  initial  ensemble  that  is  removed  by  dissociation 
after  a  large  number  of  sequential  collisions.  Several  trends  are  immediately 
evident. 

Dissociation  probabilities  increase  monotonically  as  the  thermal  bath 
temperature  is  increased.  For  initially  non-rotating  molecules,  the  dissoci¬ 
ation  probability  increases  as  the  ozone  internal  energy  is  increased.  The 
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Table  1.  Computed  Dissociation  Probabilities 


Ensemble 

Energy 

(kcal/mol) 

200 

Bath  Temperature  (K) 

750 

1000 

300 

400 

500 

■ 

0.05 

0.23 

0.29 

0.34 

0.55 

0.61 

Ej  -  26 

^  0#0 

0.37 

E-j.  -  26 

Er  =  2.0 

0.13 

Ej  -  27 

*  0*0 

0.76 

Ej  -  27 

*  1  •  0 

0.43 

effect  of  rotation  is  to  decrease  the  dissociation  probability,  i.e.,  for 
molecules  of  the  same  total  internal  energy,  vibration  is  substantially  more 
efficient  in  inducing  dissociation  than  rotation. 

The  primary  result  of  this  effort  is  to  demonstrate  the  relative  effects 
of  rotational  and  vibrational  energy  in  collisional  dissociation  of  excited 
ozone.  Vibrational  energy  is  clearly  more  effective  in  enhancing  dissociation 
than  rotational  energy. 

For  molecules  excited  within  1-2  kcal/mol  of  the  dissociation  limit, 
dissociation  probabilities  can  be  quite  large.  This  can  have  a  substantial 
effect  on  computational  models  for  ozone  formation  in  the  upper  atmosphere. 
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1 .  INTRODUCTION 


In  the  following  report  we  present  a  classical  trajectory  study  of  the 
energy  transfer  processes  occurring  In  collisions  between  highly 
vlbratlonally-rotatlonally  excited  oxygen  molecules  and  oxygen  atoms.  This 
system  has  several  features  that  make  its  study  Interesting. 

Most  theoretical  studies  of  Inelastic  processes  In  atom-diatom  collisions 
have  focussed  on  rare  gas  atom-diatom  collisions.  This  Is  probably  so  because 
reasonably  adequate  potential  energy  surfaces  may  be  constructed.  The  general 
features  of  the  Interactions  for  such  systems  are  characterized  by  a  weak 
long-range  attraction  and  a  strong  short-range  repulsion.  In  such  systems 
facile  energy  transfer  only  occurs  beti^een  rotation  and  translation.  Clas¬ 
sical  trajectory  studies  Indicate  vibration  Is  lost  very  Inefficiently  and 
often  only  through  Intermediate  transfer  to  rotation.  The  energy  transfer  is 
dominated  by  small  energy  transfers. 

A  different  sort  of  behavior  is  anticipated  for  systems  with  attractive 
short  range  forces.  In  systems  such  as  O-O2,  in  which  a  stable  three-body 
configuration  exists,  qualitatively  different  behavior  Is  anticipated.  In 
these  systems  relaxation  can  be  dominated  by  large  transfers  resulting  from 
the  strong  attractive  forces.  Additionally,  energy  transfer  may  be  facil¬ 
itated  by  complex  formation  and  exchange  collisions.  Of  particular  Interest 
here  Is  the  rapid  relaxation  of  excited  recombined  oxygen  molecules  by  oxygen 
atoms  In  the  upper  atmosphere. 

In  the  upper  atmosphere  oxygen  atoms  are  produced  by  a  number  of  chemical 
and  photochemical  processes.  These  atoms  can  then  recombine  by  a  three-body 
mechanism  to  form  oxygen  molecules.  Theoretical  and  experimental  studies  on 
diatom  recombination  Indicate  that  such  recombined  diatoms  will  be  very  highly 
Internally  excited.  The  post-recomblnatlon  relaxation  of  oxygen  molecules  Is 
of  considerable  Interest  to  the  energy  and  chemical  balance  In  the  upper 
atmosphere.  A  primary  purpose  of  this  work  Is  to  elucidate  energy  transfer 
processes  that  might  occur  during  the  recomblnatlon/relaxatlon  process. 
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2.  POTENTIAL  ENERGY  SURFACE 


The  dynamics  of  the  O-O2  system  were  computed  using  the  Varandas-Murrell 
ground  state  potential  surface  for  ozone.  This  surface  was  derived  by  Murrell 
et  al.^  using  available  spectroscopic  and  ab  Initio  data.  It  has  been  used 
earlier  by  Varandas  and  Murrell^ to  study  the  dynamics  of  the  0^02  exchange 
reaction  under  thermally  equilibrated  conditions.  Additionally  this  potential 
surface  has  been  used  as  the  ozone  portion  of  potential  surface  for  computing 
Ar-03  trajectories.^  The  detailed  form  of  the  Varandas-Murrell  ozone  poten¬ 
tial  surface  Is  given  In  Refs.  2  and  3. 
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3.  CALCULATIONS 


The  technology  of  classical  trajectory  calculations  has  become  standard. 
Only  details  relevant  to  the  calculations  presented  here  will  be  discussed. 

For  detailed  descriptions  of  the  classical  trajectory  technique  the  reader  Is 
referred  to  the  review  articles  of  Bunker'*  and  Porter  and  Raff.^  Trajectories 
were  computed  using  a  fixed  step  Adams-Moulton  integrator  with  a  Runge-Kutta 
integrator  for  the  first  five  steps.  The  accuracy  of  the  trajectories  was 
checked  by  observing  conservation  of  energy  and  angular  momentum.  Generally 
both  energy  and  angular  momentum  were  conserved  to  better  than  one  part  In 
10^.  Selected  trajectories  were  back  Integrated  to  Initial  conditions. 

Two  types  of  calculations  were  performed  In  this  study:  1)  single 
collision  studies  to  systematically  Investigate  energy  transfer  In  O-O2 
collisions  as  a  function  of  relative  translational,  vibrational  and  rotational 
energy  and  2)  multiple  colllslonal  calculations  In  which  the  evolution  of  an 
ensemble  of  O2  molecules  is  followed  by  computing  the  effects  of  sequences  of 
0-atom  collisions  on  the  ensemble  properties. 

Single  collision  calculations  were  performed  at  relative  translational 
energies  of  0.24,  0.48,  0.96  and  1.91  kcal/mol,  at  Initial  diatom  vibrational 
energies  of  2.39,  4.78,  7.17,  9.56,  14.33  and  21.50  kcal/mol  and  at  initial 
diatom  rotational  energies  of  0,  0.24,  0.48,  0.96  kcal/mol.  At  least  500  tra¬ 
jectories  were  computed  for  each  Initial  value  of  relative  translational, 
rotational  and  vibrational  energy.  Additionally,  several  calculations  were 
performed  at  fixed  Initial  diatom  Internal  energy  with  varying  partitioning 
between  Initial  rotational  and  vibrational  energies.  Total  Initial  diatom 
Internal  energy  was  fixed  at  21.50  kcal/mol  with  Initial  vibrational  energies 
of  16.72,  11.94,  11.67,  2.39  and  0.0  kcal/mol.  These  calculations  were  per¬ 
formed  at  relative  translational  energies  of  0.48  and  0.96  kcal/mol.  Batches 
of  five  hundred  or  more  trajectories  were  computed  for  each  set  of  Initial 
conditions. 
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Multiple  collision  calculations  were  performed  to  study  the  relaxation  of 
an  ensemble  of  internally  excited  oxygen  molecules  in  a  bath  of  oxygen  atoms. 
Given  the  model  nature  of  this  study  we  assumed  that  only  collisions  occurring 
on  the  ground  state  potential  surface  need  be  considered.  Additionally,  only 
collisions  between  the  excited  oxygen  molecules  and  oxygen  atoms  are  consid¬ 
ered.  In  effect  we  assume  a  dilute  oxygen  molecule  concentration  in  an  oxygen 
atom  thermal  bath.  This  situation  may  well  be  relevant  to  the  initial  stages 
of  relaxation  of  a  recombining  0-atom  gas.  The  atom  gas  was  assumed  at  a 
fixed  temperature.  The  methodology  of  these  multiple  collision  calculations 
has  been  discussed  in  detail  elsewhere.^  Some  brief  comments  are  given 
later. 


Calculations  were  performed  for  initial  diatom  internal  energies  equal  to 
26.28  kcal/mol.  Several  atomic  bath  temperatures  were  considered:  300,  400, 
500,  750  and  1000  K.  In  most  calculations  the  initial  diatom  rotational 
energy  was  set  equal  to  zero  so  that  all  the  initial  diatom  energy  was  present 
as  vibrational  energy.  Some  calculations  were  performed  with  initial  diatom 
rot«tional  energy  equal  to  2.39  kcal/mol.  For  this  case  relaxation  was  fol¬ 
lowed  at  thermal  bath  temperatures  of  300  and  400  K. 

Batches  of  250  oxygen  diatoms  were  used  to  simulate  ensemble  relaxation. 
The  relaxation  was  followed  through  at  least  25  successive  collisions  with  the 
atom  bath.  After  each  collision  tlie  diatom  internal  energy  and  angular  momen¬ 
tum  were  saved  and  used  as  initial  conditions  for  the  next  collision.  Rela¬ 
tive  translational  energies  for  each  collision  were  computed  assuming  a 
Boltzmann  distribution.  With  this  procedure  the  only  memory  of  previous  col¬ 
lisions  is  contained  in  the  diatom  internal  coordinates.  No  memory  is 
retained  in  the  0-atom  bath.  This  is  reasonable  for  a  dilute  concentration  of 
diatoms. 
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4.  SINGLE  COLLISION  CALCULATIONS:  PLSULTS  AND  DISCUSSION 


Average  energy  transferred  per  collision  as  a  function  of  initial  diat'^m 
vibrational  energy  at  several  relative  translational  energies  is  shown  in 
Figure  1.  These  calculations  were  performed  with  the  diatom  initial  rota¬ 
tional  energy  equal  to  zero.  In  these  and  all  other  calculations  presented  in 
this  paper  the  maximum  impact  parameter  was  set  equal  to  4A.  Two  points  are 
noteworthy.  The  magnitude  of  the  diatom  energy  loss  per  collision  (i.e., 
energy  transferred  from  the  diatom  Internal  to  relative  translational  energy) 
increases  as  the  initial  diatom  vibrational  energy  is  increased.  This 
behavior  is  clearly  shown  in  Figure  1.  For  a  given  value  of  initial  diatom 
vibrational  energy  the  diatom  energy  loss  decreases  monotonically  with 
increasing  relative  translational  energy.  Both  these  variations  are  readily 
understood  in  terms  of  the  relative  excess  of  diatom  initial  internal  energy 
compared  to  the  initial  translational  energy.  Molecules  having  more  initial 
vibrational  energy  may  transfer  more  energy  into  translation.  Higher  transla¬ 
tional  energy  reduces  the  relative  excess  of  diatom  internal  energy. 


INITIAL  DIATOM  VIBRATIONAL  ENERGY  (kcal/mol) 

A-3261 

Figure  1.  Energy  Transferred  from  O2  i  \  an  O2-O  Collision  as  a  Function  of 
Initial  O2  Vibrational  Energy.  Results  are  shown  for  several 
values  of  relative  translational  cnorgy,  E^. 
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The  effect  of  initial  diatom  rotational  energy  on  average  energy  trans¬ 
ferred  per  collision  is  more  complicated  and  depends  on  the  relative  transla¬ 
tional  energy.  At  the  lowest  translational  energy  considered  (0.24  kcal/mol) 
initial  diatom  rotational  energy  substantially  reduces  the  average  energy 
transferred  per  collision.  This  is  shown  in  Figure  2.  However,  at  the  high¬ 
est  considered  value  of  relative  translational  energy,  1.91  kcal/mol,  the 
initial  diatom  of  rotational  energy  increases  the  average  energy  transferred 
(see  Figure  3).  At  a  relative  translational  energy  of  0.96  kcal/mol  and 
initial  diatom  vibrational  energy  of  2.3  kcal/mol  an  increase  in  diatom  rota¬ 
tional  energy  has  little  effect  on  the  magnitude  of  the  energy  transferred 
(Figure  3). 


A-3262 

Figure  2.  Energy  Transferred  from  O2  in  an  O2-O  Collision  as  a  Function  of 
initial  O2  Rotational  Energy.  Results  are  shown  for  two  values 
of  initial  relative  translational  energy,  E^,  and  initial 
vibrational  energy,  E^. 
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INrTI*L  DIATOM  ROTATIONAL  ENERGY  (kcal/mol) 

A-3263 

Figure  3.  Energy  Transferred  from  O2  in  an  '^2“^  Collision  as  a  Function  of 
Initial  O2  Rotational  Energy  for  Several  Values  of  Initial  Diatom 
Relative  Translational  Energy,  E^.,  and  Vibrational  Energy,  Ey 

These  calculations  indicate  that  in  this  strongly  interacting  system 
rotational  energy  decreases  the  efficiency  of  energy  transfer  when  the  mole¬ 
cules  are  relatively  vibrationally  hot  compared  to  translation.  However,  at 
lower  relative  vibrational  excitations  {i.e.,  higher  translational  energies) 
the  behavior  is  reversed  and  increasing  rotational  energy  increases  energy 
transfer  from  the  diatom. 


Calculations  of  energy  transfer  from  highly  vibrationally-rotationally 
excited  O2  molecules  were  performed  with  fixed  total  internal  energy  equal  to 
21.5  Iccal/mol  for  several  initial  partitionings  between  rotational  a  id  vibra¬ 
tional  energies.  Since  vibrational-rotational  modes  are  not  stricnly  sepa¬ 
rable  (especially  at  higher  internal  energies)  it  is  more  correct  to  state 
that  we  are  examining  the  variation  of  energy  transfer  with  diatom  angular 
momentum.  The  initial  diatom  rotational  energies  were  computed  assuming  the 
diatom  at  its  equilibrium  internuclear  distance.  The  results  of  these  calcu¬ 
lations  are  presented  in  Figure  4. 
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TOTAL  02  ENERGY  -  21.50  kcal/mol 


INITIAL  DIATOM  ROTATIONAL  ENERGY  (kcal/mol) 

A-3264 

Figure  4.  Energy  Transferred  from  O2  In  an  O2-O  Collision  as  a  Function  of 
Initial  O2  Rotational  Energy  with  Total  O2  Initial  Interval 
Energy  Equal  to  21.50  kcal/mol.  Results  are  shown  for  two 
values  of  relative  translational  energy,  E^. 

At  the  lower  relative  translational  energy  of  0.48  kcal/mol  the  average 
energy  transferred  decreases  rapidly  as  rotational  energy  is  added  to  the 
diatom.  Above  0.5  kcal/mol  Initial  rotational  energy  the  average  energy 
transferred  becomes  essentially  constant  and  Indepedent  of  the  Initial  parti¬ 
tioning  between  vibration  and  rotation.  At  1.91  kcal/mol  relative  transla¬ 
tional  there  Is  similar  behavior  although  the  decrease  within  Initial  rota¬ 
tional  energy  Is  substantially  less.  The  difference  In  rotational  energy 
dependences  at  these  two  translational  energies  has  the  following 
consequences. 

At  zero  or  very  low  Initial  rotational  energies,  energy  transferred 
decreases  with  Increasing  translational  energy.  However,  at  high  Initial 
rotational  energies,  energy  transfer  increases  with  Increasing  translational 
energy. 


830 


Histograms  o£  post  collision  diatom  Internal  energy  transfers  are 

presented  In  Figures  5  and  6.  Note  that  the  total  classical  cross  section, 

2  2 

Jib  ,  equals  50. 3A  .  In  the  cases  shown  In  Figures  5  and  6  the  Initial 
max 

diatom  Internal  energy  Is  equal  to  21.5  kcal/mol  and  Initial  translational 
energy  1.91  kcal/mol.  For  the  calculations  in  Figure  5  tha  initial  rotational 
energy  Is  equal  to  zero  and  for  those  In  Figure  6  initial  rotational  energy 
equals  4.78  kcal/mol.  In  both  cases  there  is  a  very  large  proportion  of 
extremely  large  energy  transfers.  Energy  transfers  of  more  than  10  kcal/mol 
are  quite  frequent.  The  cross  section  for  energy  transfers  greater  than 
10  kcal/mol  is  7A^  for  the  case  in  Figure  5  and  4A^  for  that  in  Figure  6. 

This  behavior  is  very  dif rerent  than  seen  in  rare  gas  atom-diatom  collisions. 
In  rare  gas  atom-diatom  collisions  large  energy  transfers  are  extremely  rare 
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Figure  5.  Distribution  of  O2  Internal  Energy  Transfers  in  a  Collision 
with  O.  Collision  initial  conditions  are:  O2  internal 
energy  equals  Zt.jO  kcal/mol,  O2  rotational  energy  equals 
zero,  and  relative  translational  energy  equals  1.91  kcal/mol. 
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Figure  6.  Distribution  of  O2  Internal  Energy  Transfers  in  a  Collision  with  0. 
Collision  initial  conditions  are:  O2  internal  energy  equals 
21.50  kcal/mol,  O2  rotational  energy  equals  4.78  kcal/mole,  and 
relative  translational  energy  equals  1.91  kcal/mol. 

(in  the  translational  energy  range  under  consideration).  Small  energy  trans¬ 
fers  dominate.  That  behavior  of  course  accounts  for  the  relatively  slow 
relaxation  of  vibratlonally  excited  diatoms  in  a  rare  gas  atomic  bath.  In 
this  case  we  anticipate  very  rapid  relaxation  of  excited  diatoms  due  to  the 
large  average  energy  transferred  per  collision  and  the  high  frequency  of  very 
large  energy  transfers. 

The  differences  between  the  0-02  relaxation  and  a  system  such  as  Ar-02 
are  clearly  attributable  to  the  much  stronger  interactions  in  the  O-O2  system. 
Experience  with  rare-gas  atom-molecule  collisions  indicates  that  the  transfer 
of  vibration  to  translation  is  not  facile.  In  fact  energy  transfer  often 
seems  to  occur  by  a  vlbration-to-rotation-to-translation  pathway  in  highly 
excited  diatom  and  polyatomic  relaxation  by  rare  gas  atoms.  Similar  behavior 
has  been  noticed  in  the  relaxation  of  other  atom-diatom  systems  in  which 
strong  interactions  are  present.  Osborn  and  Smith^  recently  investigated 


vibrational  energy  transfer  in  a  model  A+BC(v)  system.  They  concluded  that 
even  fairly  moderate  intermolecular  attraction  can  appreciably  enhance  the 
degree  of  vibrational  energy  transfer  and  lead  to  large  energy  transfers. 

This  type  of  interaction  is  clearly  found  on  the  Murrell-Sorbie-Varandas  ozone 
potential  surface. 

We  also  found  that  the  role  of  long-lived  collisions  is  minor  in 
Initially  excited  O2-O  collisions.  At  the  highest  internal  excitation 
energies  virtually  no  collisions  longer  than  1  ps  were  observed.  At  the 
lowest  energies  the  collision  lifetime  Increased  substantially.  However,  the 
observed  energy  transfer  for  longer-lived  collisions  appeared  similar  to  those 
for  more  direct  collisions.  This  again  is  in  substantial  agreement  with  the 
implications  of  the  Osborn  and  Smith  model  study.  We  have  also  compared 
energy  transfer  in  reactive  and  non-reactlve  collisions.  Histograms  of 
colllslonal  energy  transfers  for  reactive  and  non-reactive  encounters  are 
shown  in  Figures  7  and  8,  respectively.  These  histograms  are  for  the  case 
shown  in  Figure  5.  For  larger  energy  transfers,  i.e.,  greater  than  3  kcal/mol 
transferred  from  diatom  internal  energy  to  relative  translational  energy,  the 
distributions  for  reactive  and  non-reactive  encounters  are  qualitatively 
similar.  The  non-reactive  case  is  of  course  dominated  by  collisions  that  pro¬ 
duce  little  or  no  energy  transfer.  A  possible  Interpretation  of  the  similar¬ 
ity  between  the  histograms  to  larger  energy  transfers  is  that  the  magnitude  of 
the  energy  transfer  is  controlled  by  the  strength  of  the  interactions  in  the 
three-body  complex  and  not  by  a  peculiarity  of  the  atom  exchange  process.  A 
few  points  on  the  dependence  of  the  exchange  probability  on  internal  excita¬ 
tion  and  relative  translational  energy  are  appropriate. 

At  zero  initial  diatom  rotational  energy  the  exchange  probability 
increases  monotonically  with  increasing  initial  vibrational  excitation.  Also, 
the  exchange  probability  decreases  with  increasing  translational  energy.  This 
is  shown  clearly  in  Figure  9.  An  Increase  in  initial  rotational  energy  causes 
a  rapid  decrease  in  the  exchange  cross  section.  This  behavior  is  shown  in 
Figure  10  for  initial  diatom  vibration  excitations  of  2.4  and  4.8  kcal/mol. 

The  effegt  of  rotation  is  greatest  at  the  lowest  translational  energies.  The 
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Figure  7.  Relative  Frequency  of  Energy  Transfers  in  Reactive  O2-O  Encounters 
Initial  collision  conditions  same  as  in  Figure  5 
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Figure  8.  Relative  Frequency  of  Energy  Tranfers  in  Non-reactive  O2-O 
Encounters.  Initial  collision  conditions  same  as  in 
Figure  5. 
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calculations  involving  diatoms  initially  excited  at  21.5  kcal/mol  show  similar 
decreases  in  exchange  probability  with  increasing  rotational  er.crg^^.  Vari¬ 
ation  in  exchange  probability  for  initial  total  diatom  excitation  of 
21.50  kcal/mol  with  varying  partitioning  between  rotation  and  vibration  is 
shown  in  Figure  11  for  two  values  of  initial  relative  translational  energy, 
0.48  and  1.91  kcal/mol.  The  variation  in  exchange  probability  with  initial 
rotational  energy  is  again  much  greater  at  the  lower  relative  translational 
energy.  A  possible  physical  explanation  for  this  behavior  is  as  follows.  At 
zero  or  very  low  initial  diatom  rotational  energies,  the  motion  of  the  collid¬ 
ing  pair  is  dominated  by  the  long  range  molecular  attraction  on  the  ozone 
potential  surface  that  orients  the  triatom  system  into  a  lowest  potential 


INITIAL  DIATOM  ROTATIONAL  ENERGY  (kcal/mol) 
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Figure  11.  Oxygen  Atom  Exchange  Probability  as  a  Function  of  O2  Initial 
Rotational  Energy  for  a  Fixed  Initial  O2  Internal  Energy 
Equal  to  21.50  kcal/mol.  Results  are  shown  for  two  values  of 
relative  translational  energy,  E^. 
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energy  configuration  approximating  an  ozone  molecule.  In  such  a  configu¬ 
ration,  exchange  can  occur  without  overcoming  an  energy  barrier.  The  Murrell- 
Sorbie-Varandas  potential  surface  is  highly  directional,  i.e.,  configurations 
deviating  substantially  from  the  lowest  energy  path  to  form  the  triatomic  are 
quite  repulsive.  As  the  Initial  diatom  rotational  energy  is  increased  the 
long  range  attraction  cannot  orient  the  diatom  as  effectively.  Therefore  more 
collisions  occur  in  configurations  where  exchange  is  energetically  difficult. 
Thus  the  reaction  probability  decreases. 

This  effect  should  be  greatest  for  lower  relative  translational  energies 
since  the  ability  of  the  diatom  to  orient  will  decrease  as  the  collision  time 
decreases  at  higher  collision  velocities.  Therefore,  this  picture  is  also 
consistent  with  the  general  decrease  in  exchange  probability  with  transla¬ 
tional  energy.  For  higher  translational  energies  interferences  due  to  diatom 
rotation  should  be  relatively  less  Important  since  the  diatom  has  less  time  to 
orient  itself  into  a  favorable  configuration  for  exchange.  A  less  dramatic 
decrease  with  rotational  energy  is  therefore  expected  at  higher  translational 
energies . 

This  picture  may  also  be  used  to  rationalize  the  general  behavior  of  the 
average  energy  transfer  per  collision  with  relative  translational  and  diatom 
rotational  energy.  Large  energy  transfer  occurs  for  collisions  that  feel  the 
strong  attractive  forces  of  the  triatom  system  near  its  equilibrium  configu¬ 
ration.  Collisions  that  occur  at  configurations  where  the  forces  are  repul¬ 
sive  should  be  qualitatively  similar  to  rare-gas  atom  diatom  collsions  and 
result  in  smaller  energy  transfers.  According  to  this  picture  collisions  at 
low  translational  energies  and  rotational  energies  will  allow  a  larger  frac¬ 
tion  to  probe  the  strong  attractive  forces.  Therefore,  energy  transfer  should 
decrease  with  increasing  translational  energy  and  diatom  rotational  energy. 


5.  CALCULATIONS  OF  STATE-TO-STATE  CROSS  SECTIONS 


The  present  calculations  were  done  for  purely  classical  atoms.  However, 
we  may  use  a  boxing  procedure  to  estimate  quantal  cross  sections.  First  we 
consider  the  case  In  which  molecules  are  excited  21.5  kcal/mol  above  the  well 
minimum.  The  vibrational  levels  of  the  oxygen  molecule  are  listed  In  Table  1. 
These  values  are  derived  from  the  spectroscopic  constants  given  by  Herzberg^ 
assuming  oxygen  can  be  described  adequately  as  a  harmonic  oscillator. 


Table  1.  Vibrational  Energy  Levels  of  O2  (kcal/mol) 


v-0 

2.258 

v-1 

6.773 

v-2 

11.29 

v«3 

15.80 

v»4 

20.32 

v-5 

24.84 

The  Initial  diatom  vibrational  energy  of  21.50  corresponds  to  an  energy 
slightly  above  the  value  for  v*4.  We  assume  that  It  adequately  represents 
v-4.  Calculations  having  Initial  vibrational  energy  of  7.17  kcal/mol  will  be 
used  to  describe  V"1.  To  compute  state  to  state  cross  sections  we  assume  a 
quaslclasslcal  boxing  procedure  as  described  by  Porter  and  Raff.^  The  purpose 
of  the  effort  was  not  primarily  to  compute  cross  sections  but  to  follow  the 
evolution  of  ensemble  of  relaxing  molecules.  This  procedure  although  somewhat 
approximate  will  yield  reasonable  estimates. 

The  cross  section  of  Interest  Is  the  total  cross  section  for  a  transition 
between  two  vibrational  states,  a(v,v'),  without  regard  to  final  rotational 
state.  It  Is  given  by 

a(v;v')  ■  L  a(v,J;  v',J*) 

J' 
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To  compute  the  final  state  vibrational  energy  we  write  the  diatom  rotational 
energy  as 


BJ 


2 


» 


where  J  is  the  diatom  rotational  angular  momentum  and  B  is  the  equilibrium 
rotational  constant.  The  vibratxonal  energy  transfer,  AE^,  is  then  given  by 


AE 


V 


E 


f 


V 


d  d  R  R 


where  E^  and  are  diatom  Initial  and  final  vibrational  energies,  E^  and  E^ 

V  V  d  d 

are  the  diatom  initial  and  final  rotational  energies. 


The  computed  cross  sections  for  transitions  from  the  v-4  oxygen  vibra¬ 
tional  state  as  a  function  of  relative  translational  energy  are  shown  in 
Figure  12.  Values  for  transition  from  v=2  are  shown  in  Figure  13.  The  values 
for  these  cross  sections  are  quite  large.  A  general  trend  for  decreasing 
cross  section  with  increasing  translational  energy  can  be  seen  in  Figures  12 
and  13. 


We  can  also  calculate  the  rate  constant  for  each  transition  using  the 
»  equation 


k  »  vc 

This  quantity  is  probably  of  more  Interest  to  modelers  trying  to  simulate 
upper  atmospheric  behavior.  The  computed  values  are  listed  in  Table  2.  The 
rate  constants  are  listed  as  a  function  of  approximate  temperature  corre¬ 
sponding  to  the  relative  translational  energy. 


839 


Figure  12.  Computed  Oxygen  Vibrational  Deexcitation  Cross  Sections, 
a(i+j),  for  an  O2-O  Collision  as  a  Function  of  Relative 
Translational  Energy 


RELATIVE  TRANSLATIONAL  ENERGY  (kcal/mol )  A-3273 

Figure  13.  Computed  v>2  to  v*0  Oxygen  Vibrational  Deexcitation  Cross 
Section,  a(2-»-0),  for  an  O2-O  Collision  as  a  Function 
of  Relative  Translational  Energy 
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Table  2.  Computed  Vibrational  Deexcitation  Rate  Constants 


60  K 

120  K 

240  K 

480  K 

k(4-»-0) 

4.3x10"1 1 

5.2x10-11 

4.5x10-11 

7.7x10-1  1 

k  {  4-f  1 ) 

4. 1x10-11 

4.3x10-1 1 

6.1x10-11 

6.7x10-11 

k(4>2) 

2.9x10-1 1 

4.0x10-11 

6.1x10-11 

7.7x10-11 

k(4+3) 

2.4x10-11 

3.9x10-11 

5.8x10-1  1 

5.4x10-11 

k( 1+0) 

9.9x10-11 

1.2x10-10 

1.6x10-10 

1.4x10-10 

6.  MULTIPLE  COLLISION  CALCULATIONS 


The  multiple  trajectory  technique  ie  based  on  an  approximate  solution  to 
the  Boltzmann  equation  given  as  follows 


dt 


Lf 


(1) 


where  L  Is  the  Boltzmann  collision  operator  and  £  Is  the  ensemble  distribution 
function.  The  essence  of  the  multiple  collision  technique  Is  1)  to  approxi¬ 
mate  a  solution  to  Eq.  (1)  as 


f(t+T)  -  (1+LT)f(t) 


(2) 


where  t  Is  chosen  to  correspond  to  a  mean  collision  time;  2)  to  Identify  the 
term  (1-»-LT)£(t)  as  the  ensemble  after  one  collision  of  each  member  of  the 
ensemble;  and  3)  to  use  Monte  Carlo  trajectory  calculations  to  compute  the 
collision  dynamics  of  each  member  of  the  ensemble.  The  simulation  Is  accom¬ 
plished  by  picking  a  finite  ensemble  of  molecules  (In  this  case  250}  in  a 
given  Initial  state.  Each  molecule  undergoes  a  sequence  of  collisions.  The 
Initial  conditions  for  the  collision  are  chosen  according  to  the  ensemble 
thermal  characteristics.  The  evolution  Is  constructed  by  using  the  molecule 
post  collision  variables,  i.e.,  energy  and  angular  momentum,  as  initial  condi¬ 
tions  for  the  next  collision.  After  n  collisions  the  ensemble  describes  the 
evolution  after  nr.  The  primary  quantities  of  Interest  to  this  study  are  the 
average  total  diatom  energy,  vibrational  energy  and  rotational  energy. 

We  first  discuss  the  evolution  of  mean  oxygen  diatom  quantities.  The 
oxygen  diatom  relaxation  Is  assumed  to  occur  In  a  bath  of  thermally  equili¬ 
brated  oxygen  atoms  with  a  small  concentration  of  diatoms.  Therefore  diatom- 
diatom  collisions  are  neglected.  It  Is  furthermore  assumed  that  no  transla¬ 
tional  motion  memory  Is  retained  between  collisions  so  that  the  relative 
cranslatlonal  energy  of  distribution  of  colliding  pairs  Is  given  by  a 
Boltzmann  distribution. 
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The  evolution  of  the  average  total  internal  diatom  energy,  diatom  vibra¬ 
tional  and  diatom  rotational  energies  are  shown  In  Figures  14  through  20.  For 
the  cases  shown  In  Figures  14  through  18  the  Initial  rotational  energy  was 
equal  to  zero.  For  the  cases  In  Figures  19  and  20  the  Initial  rotational 
energy  was  set  equal  to  2.39  kcal/mol.  In  all  cases  the  Initial  Internal 
energy  was  equal  to  26.3  kcal/mol.  We  first  discuss  the  cases  with  Initially 
non-rotating  diatoms. 

Calculations  were  performed  for  oxygen  atom  baths  varying  between  300  and 
1000  K.  Qualitatively  and  quantitatively  there  Is  no  dramatic  variation  In 
relaxation  observed  with  bath  temperature.  We  therefore  focus  first  on  the 
30C  K  bath  temperature  case  (Figure  14).  Note  that  we  are  considering  purely 
classical  relaxation  and  that  the  «iqulllbrlum  values  of  the  diatom  rotational 
and  vibrational  energies  are  each  kT. 


NUMBER  OF  COLLISIONS  WITH  OXYGEN  ATOM  BATH  A-3664 


Figure  14.  Variation  of  Oxygen  Diatom  Total  Internal  Energy,  Y1,  Vibrational 
Energy,  Y2,  and  Rotational  Energy,  Y3,  with  Number  of  Collisions 
of  Oxygen  Molecules  with  Oxygen  Atom  Thermal  Bath.  Initial 
conditions  are:  oxygen  Internal  energy,  Ej,  equals  26.3  kcal/mol, 
oxygen  rotational  energy,  E{^,  equals  zero,  and  oxygen  atom  bath 
temperature,  T3,  equals  300  K. 
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Figure  15.  Variation  of  Oxygen  Diatom  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Atom  Thermal  Bath.  Initial  conditions  are:  oxygen 
internal  energy,  Ej,  equals  26.3  kcal/mol,  oxygen  rotational 
energy,  E|^,  equals  0.0,  and  oxygen  atom  bath  temperature,  Tg, 
equals  400  K. 


Figure  16.  Variation  of  Oxygen  Diatom  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Atom  Thermal  Bath.  Initial  conditions  are:  oxygen 
internal  energy,  Ej,  equals  26.3  kcal/mol,  oxygen  rotational 
energy,  Ej^,  equals  0.0,  and  oxygen  atom  bath  temperature,  T3, 
equals  *  500  K. 
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Figure  17.  Variation  of  Oxygen  Diatoa  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Atom  Thermal  Bath.  Initial  conditions  are:  oxygen 
internal  energy,  Ej,  equals  26.3  kcal/mol,  oxygen  rotational 
energy,  E{^,  equals  0.0,  and  oxygen  atom  bath  temperature,  Tg, 
equals  750  K. 


Figure  18.  Variation  of  Oxygen  Diatom  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Atom  Thermal  Bath.  Initial  conditions  are:  oxygen 
internal  energy,  Ej,  equals  26.3  k.cal/mol,  oxygen  rotational 
energy,  Eg,  equals  0.0,  and  oxygen  atom  bath  temperature,  Tg, 
equals  1000  K. 
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NUH8ER  OF  COLLISIONS  MITH  OXYGEN  ATOH  BATH  A- 3661 

Figure  19.  Variation  of  Oxygen  Diatom  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Atom  Thermal  Bath.  Initial  conditions  are:  oxygen 
Internal  energy,  Ej,  equals  26.3  kcal/mol,  oxygen  rotational 
energy,  Ei^,  equals  2.39  kca/mol,  and  oxygen  atom  bath 
temperature,  Tb,  equals  300  K. 


NUMBER  OF  COLLISIONS  HITH  OXYGEN  ATOH  BATH  A-3653 

Figure  20.  Variation  of  Oxygon  Diatom  Total  Internal  Energy,  Vibrational 
Energy  and  Rotational  Energy  with  Number  of  Collisions  with 
Oxygen  Aton  Thermal  Bath.  Initial  conditions  are:  oicygen 
Internal  energy,  Ej.  equals  26.3  hcal/mol,  ojcygen  rotational 
energy,  equals  2.39  hca/mol,  and  oxygen  atom  bath 
temperature,  T3,  equals  400  K. 
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The  ensemble  internal  energy  drops  very  rapidly  for  the  first  several 
collisions.  Approximately  5  kcal/mol  internal  energy  is  lost  in  the  first 
collision.  Very  large  average  energy  transfers  persist  up  to  approximately 
the  fifth  collision.  After  t''e  fifth  collision  the  ensemble  Internal  energy 
is  approximately  12  kcal/mol,  implying  that  over  14  kcal/mol  has  been  trans¬ 
ferred  out  of  the  molecule.  Beyond  approximately  the  fifth  collision  the 
relaxation  slows  quite  dramatically.  Between  collisions  five  and  twenty-eight 
an  additional  8  kcal/mol  are  transferred  from  the  diatom.  At  the  end  of  the 
calculation,  after  the  twenty-eighth  collision,  the  average  diatom  internal 
energy  is  3.9  kcal/mol.  The  equilibrium  value  for  the  internal  energy  of  a 
classical  diatom  in  a  300  K  bath  is  1.2  kcal/mol.  Thus,  after  twenty-eight 
collisions  approximately  90  percent  of  the  excess  internal  diatom  energy  has 
been  lost  to  the  bath.  This  is  indeed  rapid  relaxation  for  a  vibrating 
molecule. 

The  average  diatom  vibrational  energy  is  initially  eqiial  to  the  total 
internal  energy  in  this  case.  It  drops  rapidly  with  collision  number.  Again 
very  large  energy  transfers  occur  in  the  first  few  collisions.  By  the  fifth 
collision  the  vibrational  energy  has  decreased  from  26.3  to  7.4  kcal/mol.  The 
average  vibrational  energy  drops  by  approximately  10  kcal/mol  in  the  first 
collision.  Between  collisions  five  and  twenty-eight  the  vibrational  energy 
decreases  much  more  slowly  from  7.4  to  2.0  kcal/mol.  The  initial  very  rapid 
transfer  of  vibrational  energy  is  associated  with  two  processes:  inter- 
molecular  transfer  of  vibrati  '''  energy  and  intramolecular  transfer  between 
vibrational  and  rotational  ener  ^ 

In  the  case  under  consideration  (Figure  14)  the  initial  rotational  energy 
is  equal  to  zero.  After  one  collision  the  average  diatom  rotational  energy  is 
4.0  kcal/mol.  The  value  of  the  average  rotational  energy  peaks  at  the  second 
collision  at  a  value  of  5.2  kcal/mol.  After  the  second  collision,  the  average 
rotational  energy  decreases.  Note  that  the  rotational  relaxation  does  not 
decouple  from  the  vibrational  energy  relaxation  over  the  time  scales  consid¬ 
ered.  For  this  classical,  strongly  interacting  system  rotational  and  vibra¬ 
tional  energies  relax  at  the  same  rates.  By  the  twentieth  collision  the 
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average  rotational  and  vibrational  energies  are  approximately  equal.  From 
that  collision  on,  both  modes  relax  at  the  same  rate.  This  behavior  is 
observed  in  all  the  cases  considered  over  the  computed  relaxation  times.  A 
plausible  explanation  for  this  behavior  is  the  large  coupling  between  rota¬ 
tional  and  vibrational  motion  on  the  Murrell-Sorbie-Varandas  ozone  potential 
surface.  This  strong  coupling  does  not  allow  separation  of  time  scales  for 
rotational  and  vibrational  relaxation. 

« 

In  the  more  frequently  studied  cases  of  rare  gas  atom-diatom  relaxation, 
the  Interaction  potentials  may  be  separated  into  terms  separately  leading  to 
vibrational  and  rotational  energy  transfers.  The  coupling  between  diatom 
vibration  and  translational  motion  tends  to  be  very  inefficient  in  the  energy 
range  considered  here  and  the  diatom  vibrational  energy  relaxes  very  slowly. 
Rotational  coupling  to  translational  motion  is  generally  quite  efficient 
resulting  in  rapid  rotational  relaxation.  Thus,  at  least  for  low  values  of 
vibrational  excitation  rotational  relaxation  occurs  on  a  faster  time  scale. 

This  behavior  is  generally  not  found  for  high  values  of  vibrational  excitation 
even  in  rare  gas  atom-diatom  systems.^  Therefore  for  low  collision  numbers 
strong  coupling  is  expected.  However,  the  strong  coupling  that  exists  even 
for  vibrational  energies  only  a  few  kilocalories  above  the  well  minimum  would 
seem  due  to  the  strong  interactions  in  the  O-O2  system. 

An  additional  feature  of  Interest  is  the  number  of  collisions  required 
for  the  rotational  energy  to  build  up  in  the  molecule  by  collisions  with  the  , 

bath.  This  is  the  number  of  collisions  for  the  rotational  energy  to  reach  a 
peak  value  before  beginning  to  relax.  As  is  readily  seen  in  Figures  14 
through  18  this  occurs  in  three  to  five  collisions  with  the  bath.  This  number 
of  collisions  corresponds  to  an  intramolecular  energy  transfer  colllsional 
relaxation  time  and  is  very  fast.  The  rapidity  of  this  transfer  is  indicative 
again  of  the  strong  coupling  between  vibration  and  rotation  on  the  Murrell- 
Sorbie-Varandas  ozone  potential  surface. 

It  is  of  interest  to  see  whether  or  not  a  single  relaxation  time  is 
sufficient  to  characterize  the  molecular  relaxation.  This  would  seem  unlikely 
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for  the  non-rotating  initial  oxygen  molecule  ensemble  due  to  the  strong  depen¬ 
dence  of  total  energy  transfer  on  rotational  energy.  As  discussed  earlier  in 
single  collision  calculations  for  highly  excited  molecules,  increasing  rota¬ 
tional  energy  causes  a  sharp  decrease  in  the  average  collislonal  energy  trans¬ 
fer  (Figure  4)  at  low  values  of  relative  translational  energy.  Therefore  for 
the  initially  non-rotating  cases  the  early  time  relaxation  will  be  character¬ 
ized  by  large  energy  transfers  characteristic  of  non-rotating  molecules. 

Within  three  to  five  collisions  with  the  bath,  rotational  energy  is  built  up 
in  the  diatoms  and  relaxation  will  occur  at  the  slower  rate  of  rotating  mole¬ 
cules.  This  can  be  significant  since  as  mentioned  earlier  for  the  300  K 
ncn-rotating  case  since  most  of  the  excess  diatom  energy'  is  luct  within  five 
collisions  with  the  bath. 

The  total  molecular  energy  in  excess  of  its  classical  equilibrium  value 
as  a  function  of  collision  number  is  shown  for  each  case  in  Figures  21 
through  27.  The  plots  are  presented  on  semi-log  paper  so  that  in  a  regime  in 
which  the  relaxation  is  described  by  a  single  exponential  decay  a  linear  curve 
should  result.  First  we  consider  the  cases  in  which  the  molecules  are  non¬ 
rotating  (Figures  21  through  25). 

Two  regimes  regions  are  readily  identified:  a  short  time  (zero  to  five 
collisions)  rapid  relaxation  regime  and  longer  time  (>five  to  ten  collisions) 
slower  relaxation.  As  mentioned  earlier  the  short  time  regime  corresponds  to 
the  time  in  which  rotational  energy  (or  angular  momentum)  is  being  transferred 
into  the  molecule.  As  noted  earlier,  (Figure  1)  collislonal  energy  transfer 
decreases  strongly  with  diatom  rotational  energy.  Therefore  the  molecules 
relax  more  rapidly  during  the  first  few  collisions  before  rotational  energy  is 
built  up.  This  effect  is  most  pronounced  at  the  lower  atom  bath  temperatures 
because  as  seen  earlier  for  non-rotating  diatoms  the  collislonal  energy  trans¬ 
fer  from  the  diatom  decreases  with  increasing  atom-diatom  trat.S4.a  wxonal 
energy.  Therefore  the  initial  molecular  relaxation  rate  (zero  to  five  colli¬ 
sions)  increases  with  decreasing  bath  temperature.  This  is  seen  in  Figures  21 
through  25.  At  longer  times  the  situation  is  reversed.  For  vibrationally- 
rotatlonally  excited  molecules  with  substantial  rotational  energy,  collislonal 
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Figure  21.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  Ei-Eeq, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
initial  conditions  are:  diatom  internal  energy,  E^,  equals 
26.3  kcal/mol,  rotational  energy,  Ef^,  equals  zero,  and  bath 
temperature,  Xb,  equals  300  K. 
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Figure  22.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  Ej-Egq, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
initial  conditions  are:  diatom  Internal  energy,  E^,  equals 
26.3  kcal/mol,  rotational  energy,  E^,  equals  0.0,  and  bath 
temperature,  Tbi  equals  400  K. 
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Figure  23.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  Ei-E^q, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
Initial  conditions  are:  diatom  Internal  energy,  Ej,  equals 
26.3  kcal/mol,  rotational  energy,  Ej^,  equals  0.0,  and  bath 
temperature,  T3,  equals  500  K. 
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Figure  24.  Oxygen  Internal  Energy  Above  the  Equlllbrlvim  Value,  E^-E^q, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
Initial  conditions  are:  diatom  Internal  energy,  E^,  equals 
26.3  kcal/mol,  rotational  energy,  Ej^,  equals  0,0,  and  bath 
temperature,  Tg,  equals  750  K. 
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Figure  25.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  Ej-Egq, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
initial  conditions  are:  diatom  internal  energy,  E^,  equals 
26.3  kcal/mol,  rotational  energy,  equals  0.0,  and  bath 

temperature,  T3,  equals  1000  K. 


energy  transfer  increases  with  increasing  relative  translational  energy. 
Therefore  after  the  initial  transient  in  which  rotational  energy  is  pumped 
into  the  molecule  the  relaxation  rate  Increases  with  bath  temperature.  This 
again  is  seen  in  Figures  21  through  25. 


In  summary  two  relaxation  regimes  are  seen  for  the  relaxation  of  highly 
excited  oxygen  diatoms  in  an  atomic  bath.  The  first  regime  is  a  short  time 
regime  persisting  until  rotational  energy  is  transferred  into  diatom.  During 
this  period  relaxation  is  very  rapid  and  the  rate  increases  with  decreasing 
bath  temperature.  A  second  relaxation  regime  occurs  after  a  rotational  steady 
state  is  built  up  in  the  diatom.  Relaxation  proceeds  at  a  slower  rate  and  the 
rate  of  relaxation  increases  with  bath  temperature. 


852 


The  behavior  of  the  two  cases  in  which  some  of  the  initial  molecular 
energy  was  put  into  rotational  energy  is  shown  in  Figures  26  and  27.  Note 
that  the  initial  value  of  the  rotational  energy  was  below  the  value  of  the 
steady  state  value  and  an  increase  in  average  diatom  rotational  energy  was 
computed  during  the  first  few  collisions.  The  average  rotational  energy 
increases  from  its  initial  value  of  2.3  k.cal/mol  to  about  4.5  kcal/mol  after 
three  collisions  for  both  the  300  and  400  K  atom  baths.  Nevertheless  the 
effect  of  initial  rotational  energy  is  readily  seen  by  comparing  Figure  21 
with  Figure  26  for  the  300  K  bath  case  and  correspondingly  Figure  22  with 
Figure  27  for  the  400  K  bath  case. 
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Figure  26.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  Ej-Egq, 
as  Function  of  Number  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
initial  conditions  are:  diatom  internal  energy,  E^,  equals 
26.3  hcal/mol,  rotational  energy,  Ej^,  equals  2.89  kcal/mol, 
and  bath  temperature,  Tg,  equals  300  K. 
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Figure  27.  Oxygen  Internal  Energy  Above  the  Equilibrium  Value,  E^-E^q, 
as  Function  of  Kumber  of  Collisions  with  Oxygen  Atom 
Thermal  Bath.  Initial  bath  and  oxygen  molecule  ensemble 
Initial  conditions  are:  diatom  Internal  energy,  Ej,  equals 
26.3  kcal/mol,  rotational  energy,  E|^,  equals  2.89  hcal/mol, 
and  bath  temperature,  Tg,  equals  400  K. 


Initial  diatom  rotational  energy  reduces  the  Initial  rapid  relaxation 
rate  In  the  early  time  regime.  The  differences  In  relaxation  vanish  after 
about  ten  collisions  by  which  time  a  rotational  energy  steady  state  Is 
achieved.  The  long  time  relaxation  rates  of  Initially  rotating  and  non¬ 
rotating  ensembles  are  Identical.  However,  nince  much  of  the  ensemble  relaxa¬ 
tion  occurs  during  the  first  few  collisions,  the  non-rotating  cases  show  an 
overall  faster  relaxation  rate.  It  Is  anticipated  that  an  ensemble  containing 
rotational  energy  equal  to  the  steady  state  value  would  relax  according  to  a 
single  exponential  decay  law. 


The  evolution  of  the  average  molecular  energy  distribution  as  a  function 
of  number  of  collisions  with  the  atomic  bath  Is  a  quantity  of  Interest.  The 
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initial  distribution  function  is  a  delta  function  at  26.3  kcal/mol.  At  equi¬ 
librium,  when  fully  relaxed  by  the  thermal  bath,  it  becomes  a  Boltzmann  func¬ 
tion.  The  molecular  energy  distributions  after  one,  five,  ten,  twenty  and 
twenty-five  collisions  are  shown  in  Figures  28  through  32  respectively,  for  an 
initially  non-rotating  diatom  in  a  500  K  atomic  bath. 

After  one  bath  collision  (Figure  28)  the  distribution  of  energies  has 
broadened  substantially  from  the  initial  delta  function  at  26.3  kcal/mol. 
Memory  of  the  initial  distribution  is  seen  from  the  large  peak  for  energies 
between  26  and  28  kcal/mol.  The  occurrence  of  large  single  energy  transfers 
is  readily  seen  by  the  large  number  of  molecules  below  16  kcal/mol.  That  part 
of  the  distribution  is  due  to  single  collision  energy  transfers  in  excess  of 
10  kcal/mol.  After  five  collisions  (Figure  29)  the  molecular  energy  distribu¬ 
tion  is  quite  flat  between  2  and  20  kcal/mol.  Memory  of  the  initial  energy 
distribution  is  not  apparent.  After  ten  collisions  (Figure  30)  the  distribu¬ 
tion  begins  to  take  on  a  Boltzmann-like  form.  The  bulk  of  the  ensemble  mole¬ 
cules  have  relatively  low  internal  energies  with  the  highest  histogram 
peak  for  molecular  energies  between  0  and  2  kcal/mol.  Note,  however,  that 
there  are  still  substantial  numbers  of  molecules  with  high  values  of  total 
energy.  At  twenty  collisions  (Figure  31)  the  high  energy  part  of  the  distri¬ 
bution  is  substantially  reduced.  This  behavior  becomes  more  apparent  at 
twenty-five  collisions  (Figure  32)  by  which  point  80  percent  of  the  molecules 
have  energies  below  6  kcal/mol.  It  is  clear  that  the  rate  of  relaxation  for 
this  O2/O  model  system  is  remarkably  fast. 
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Figure  28.  Distribution  of  O2  Molecular  Energies  for  O2  Initially  Excited 
26.3  kcal/mol  After  One  Collision  with  0-Atom  Bath  at  500  K 
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Figure  29.  Distribution  of  O2  Molecular  Energies  for  O2  Initially  Excited 
26,3  kcal/mol  After  Five  Collisions  with  O-Atom  Bath  at  500  K 
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Figure  30.  Distribution  of  O2  Molecular  Energies  for  O2  Initially  Excited 
26.3  kcal/mol  After  Ten  Collisions  with  0-Atom  Bath  at  500  K 
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Figure  31.  Distribution  of  O2  Molecular  Energies  for  O2  Initially  Excited 
26.3  kcal/mol  After  Twenty  Collisions  with  0-Atom  Bath  at  500  K 
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Figure  32.  Distribution  of  O2  Molecular  Energies  for  O2  Initially  Excited 
26.3  kcal/mol  after  Twenty-Five  Collisions  with  0-Atom  Bath  at 
500  K 
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